AD  A052593 


D D C 

r.[?f 

Vj  APR  13  1978  ' 

r>t 


I ! .T'"TT  , . 
j u ib 1 'J  1 “ 


t i 


l k»i  Vj 

A 

A 


DISTRIBUTION  IT 


Approved  lor  tv.:- 


REVISED  EDITION  PUBLISHED  UNDER 
CIVIL  WORKS  INVESTIGATIONS  PROJECT  CW-I5I 


BY  THE 

U.  S.  ARMY  ENGINEER  DISTRICT 

CORPS  OF  ENGINEERS 
SACRAMENTO,  CALIFORNIA 


JANUARY  1962 


HID  7 


i/ 


-JL 


(This  edition  supersedes  paper  on  same  subject,  dated  July  1952,  as 
distributed  with  OCE  Civil  Works  Engineer  Bulletin  52-24  and  cited 
in  EM  1 1 10-2-1450  on  Hydrologic  Frequency  Estimates) 


(13ni22AW 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Bntand) 


MENTATION  PAGE. 


Z.  OOVT>CCESSIOK  NO 

Dno-fOToi  b92D«f 


«.  title  (mdSubutw  . «■  'fci.T¥WBb  *f»=dhWWB*iob  coVereo 

jb  pTapfrave  abtcco-t  ofpo To^h^H  To  29?vl>-r,r  f*o  bg?pd  , 2ii.r<smr;t29  Y3nr,1-'Pr,'i'r 

“bNo-to  ?utb6  .anortr-ja  pnrpm  r:69T3? 

* ^t^W§Yn^^9f'-r?mrT^9vriB5ft929»ta97  ^Vk-RLSHOAG-oRG. EWort noTb^h- 

2norjfitj2urrf  ,02f'  . nornot  nrij  nr  zr'i-t*  ^r* r r» r a'rb  zuorisv  nr  2or:J;?r*t9l,06*<Brb 


7.  AuTH6Mr«)  c 71 IJ  7*  «eu  to  : •-.nnn  mu  r r ! f.'lM 

. r^vrn 

Leo  R.  Bear#*  efdfiJku?  ampfrinno  nfnrrla  brr, 


».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


12.  REPORT  DATE 

January  1962 


IS.  NUMBER  OF  PAGES 


MONITORING  AGENCY  NAME  A ADDRESS (It  dIHatant  Irom  Controlling  Ottlea) 


1S«.  DECLASSIFICATION/DOWNGRADINO 

schedule 


16.  DISTRIBUTION  STATEMENT  (ol  thla  Raport) 


Distribution  of  this  publication  Is  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  (Am  aba  tract  entered  In  Block  20,  it  different  from  Report) 


19.  KEY  WOROS  (Continue  on  rereree  aide  i/  necttiary  end  Identity  by  block  number) 


Statistical  methods,  hydrologic  aspects,  frequency  analysis,  correlation 
analysis,  regression  analysis 


20.  ABSTRACT  (Continue  on  reveree  aide  ll  neceeaery  end  Identity  by  block  number) 

The  application  of  statistics  in  hydrologic  engineering  is  described  and 
illustrated.  The  subject  matter  covers  the  following  Items: 

(1)  A concise  review  of  the  basic  concepts  of  probability  and  correlation 
analyses  that  are  applicable  In  hydrologic  engineering,  with  a guide  to 
supplemental  reading  for  further  treatment; 

(2)  Presentation  of  detailed  computation  procedures  and  supporting 


DD  , jan  73  1473  EDITION  OF  I NGV  6 S>  IS  .OBSOLETE 


UNCLASSIFIED 

S ECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Ent rrrd) 


) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOEfWh*"  Data  Bntatad) 


|20^  estimates  based  on  analysis  of  hydrologic  records  that  have  been  adjusted 
as  required  to  conform  with  selected  reference  base  conditions^  * > 

(3)  A summary  of  procedures  for  developing  •regional  1 zed*  hydrologic 
frequency  estimates,  based  on  analyses  of  hydrologic  records  available  at 
stream  gaging  stations,  adjusted  to  provide  gene^U^)Vf,l9^d-|rqg^nc^  : 
relations  that  are  considered  most  representative  of  long-period  hydrologic 
characteristics  In  various  drainage  areas  In  the  region.  Also,  illustrations 
and  explanations  of  simple  generalization  procedures  for  use  where  these  are 
adequate  and  advantageous  are  given. 

(4)  A brief  list  of  key  questions  and  simple  problems  suitable  fbr  ^ .v  ooj 
use  in  training  classes. 


":'nf  v'tMinr.r. 

rsr 

''■r'l  ra?p  f Drill 


„-r  : frty  7 r ‘r  .* r r ■ !u"  r'rf|i 


I . 7 | ?v  : r.  • r^'  «-  , D rr;i ! 0'ir-Vf1  . 1,v 

7:  - • f - r; r 


- 7 r t , ; J ? 

; ':vr: 


■ i ■*  n.d'to.'f.  ■ - *>  f r'-1 : ■ ’T  n«'f  jr.j.f  .•  •'*. 

■ 2" ■■*"* f : • 1 ? i'".-- :■>  r?  «ri  .‘••••*•'1+21*;  r •: 

■ i i f •%■'••«<•  n ‘ • • rv'-(  ' r 

• p r*  -t  r ■ ' n n **'  'f  f r*-"’  ' » ■ O ' *"•'  *-yri  . r ! i • : ' • f i • r*  ' r r-  "*  • **  dv-  f r.  -rr. 

, 4 r*r-  ■*-*  <:•  I f ~fn*  -'f*r<i; 


'mi  r‘ ; 'V’f 


'If-.  pop-’' 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  this  P AGZfWhert  Dctn  Er.tere ri 


FOIIKV'OKD 


| 

I Since  publication  of  the  original  paper  of  this  title  in  July 

1952,  the  general  concept  of  runoff  frequency  analysis  contained 
| in  the  paper  has  been  used  and  developed  extensively  by  the  Corps 

of  Engineers  and  other  organizations.  It  is  considered  appropriate 
at  this  time  to  expand  the  paper  to  include  the  new  developments 
| and  associated  computation  techniques  and  examples. 

While  there  are  yet  many  difficulties  encountered  in  applica- 
! Cion,  the  method  of  frequency  analysis  originally  proposed  (loga- 

i rithmic  Pearson  Type  III)  has  been  generally  very  successful  in 

1 the  accurate  and  rapid  determination  of  extreme  flood  frequencies. 

; It  is  considered  adequate  for  all  hydrologic  frequency  applies- 

f tiens,  and  consequently,  the  treatment  herein  is  still  restricted 

l to  the  originally  proposed  method.  In  the  interest  of  simplicity, 

| terms  and  concepts  chat  are  not  strictly  necessary  to  this  method 

j are  not  discussed,  even  though  they  may  be  in  common  use  in 

f hydrologic  statistics. 

[ bo  many  have  contributed  toward  the  material  contained  herein 

i-  that  it  is  impossible  to  acknowledge  even  the  principal  contributors. 

’ However,  most  of  the  newer  developments  originated  under  the  Civil 

j Works  Investigations  program  of  the  Corps  of  Engineers,  particularly 

; the  project  conducted  in  the  Sacramento  District  under  the  direction 

Iof  the  writer  and  under  the  generaL  administration  of  Mr.  F.  Kochis, 

Chief  of  the  Engineering  Division,  and  Mr.  A.  Gomez,  Chief  of  the 
Planning  and  Reports  branch.  Mr.  A.  E.  Cochran  in  the  Office,  Chief 
of  Engineers,  has  guided  the  program  and  has  provided  invaluable  support 
I and  constant  encouragement. 


i 

1 


TABLE  OF  CONTENTS 


Par.  Page 

Foreword 

SECTION  1 - INTRODUCTION 

1-01  NATURE  AND  OBJECTIVES  OF  STATISTICAL  ANALYSIS  1 

1-02  PURPOSE  AND  SCOPE  1 

1-03  REFERENCES  2 

1- 04  ORGANIZATION  OF  THIS  PAPER  2 

SECTION  2 - GENERAL  PROBABILITY  CONCEPTS 
AND  DEFINITIONS 

2- 01  INTRODUCTION  4 

2-02  NATURE  OF  RANDOM  EVENTS  4 

2-03  PROBABILITY  INFERENCE  5 

2-04  CU EMULATIVE  FREQUENCY  CURVES  5 

2-05  TECHNICAL  APPROACH  IN  ESTIMATING  FREQUENCY  CURVES  7 

2- 06  TERMS  AND  SYMBOLS  8 

SECTION  3 - FLOOD  PEAK 
FREQUENCY  - GRAPHICAL  METHOD 

3- 01  USE  AND  RELATIVE  ADVANTAGES  9 

3-02  GENERAL  PROCEDURE  9 

3-03  SELECTION  AND  ARRANGEMENT  OF  DATA  10 

3-04  ADJUSTMENT  TO  UNIFORM  CONDITION  11 

3-05  PLOTTING  FORMULA  12 

3-06  PLOTTING  GRID  13 

3- 07  SAMPLE  COMPUTATION  14 

SECTION  4 - FLOOD  PEAK 
FREQUENCY  - ANALYTICAL  COMPUTATION 

4- 01  USE,  LIMITATIONS  AND  RELATIVE  ADVANTAGES  15 

4-02  EQUATIONS  USED  15 

4-03  ANNUAL-EVENT  CURVES  16 

4-04  PARTIAL-DURATION  CURVES  18 

4-05  USE  OF  HISTORICAL  DATA  19 

4-06  USE  OF  FLOW  ESTIMATES  19 


I 


Par..  Page 

4-07  ABNORMAL  DRY-YEAR  EFFECTS  19 

4-08  USE  OF  SYNTEHTIC  FLOODS  20 

4- 09  SPECIAL  CEOGRAPHIC  CONSIDERATIONS  21 

SECTION  3 - FLOOD  PEAK 
FREQUENCY  - ANALYTICAL  ADJUSTMENT 

5- 01  INTRODUCTION  23 

5-02  ESTIMATING  INDIVIDUAL  EVENTS  23 

5-03  DEGREE  OF  CORRELATION  23 

5-04  ADJUSTMENT  OF  FREQUENCY  STATISTICS  24 

5-05  ADVANTAGE  OF  ADJUSTMENT  24 

5- 06  SUMMARY  OF  PROCEDURE  25 

SECTION  6 - FLOOD  VOLUME  FREQUENCY 

6- 01  NATURE  AND  PURPOSE  26 

6-02  DATA  FOR  COMPREHENSIVE  SERIES  26 

6-03  STATISTICS  FOR  COMPREHENSIVE  SERIES  27 

6-04  FREQUENCY  CURVES  FOR  COMPREHENSIVE  SERIES  28 

6-05  DURATIONS  EXCEEDING  1 YEAR  29 

6- 06  APPLICATIONS  OF  FLOOD  VOLUME-DURATION  FREQUENCIES  31 

SECTION  7 - REGIONAL  FREQUENCY  ANALYSIS 

7- 01  GENERAL  34 

7-02  USE  OF  FREQUENCY  STATISTICS  34 

7-03  SIMPLE  SCHEMES  34 

7-04  DRAINAGE-BASIN  CHARACTERISTICS  35 

7-05  CORRELATION  METHODS  35 

7-06  LINEAR  RELATIONSHIPS  35 

7-07  EXAMPLE  OF  REGIONAL  CORRELATION  36 

7-08  SELECTION  OF  USEFUL  VARIABLES  37 

7-09  USE  OF  MAPS  38 

7-10  SUMMARY  OF  PROCEDURE  38 

7-11  REGIONAL  SKEW  DETERMINATIONS  39 

7-12  SAMPLE  REGIONAL  CRITERIA  40 


SECTION  8 - FREQUENCIES  OF  OTHER 
HYDROLOGIC  FACTORS 

8-01  INTRODUCTION 

8-02  RAINFALL  FREQUENCIES 


41 

41 


II 


IUyfradaMMitulm— mi  —MMawwiai, «■ — nhi  iniiiHlmiMiaJh  u 


Par.  Page 

8-03  LOW  FLOW  FREQUENCIES  41 

8-04  HURRICANE  FREQUENCIES  42 

8-05  SEDIMENT  FREQUENCIES  42 

8- 06  CONINCIDENT  FREQUENCIES  42 

SECTION  9 - CORRELATION  ANALYSIS 

9- 01  NATURE  AND  APPLICATION  44 

9-02  CALCULATION  OF  REGRESSION  EQUATIONS  44 

9-03  THE  CORRELATION  COEFFICIENT  AND  STANDARD  ERROR  46 

9-04  SIMPLE  LINEAR  CORRELATION  EXAMPLE  43 

9-05  FACTORS  RESPONSIBLE  FOR  NON-DETERMINATION  49 

9-06  MULTIPLE  LINEAR  CORRELATION  EXAMPLE  51 

9-07  PARTIAL  CORRELATION  52 

9-OS  VERIFICATION  OF  CORRELATION  RESULTS  53 

9-09  PRACTICAL  GUIDE  LINES  53 


SECTION  10  - STATISTICAL  RELIABILITY  CRITERIA 


10-01 

FUNCTION 

55 

10-02 

RELIABILITY 

OF 

FREQUENCY  STATISTICS 

55 

10-03 

RELIABILITY 

OF 

FREQUENCY  ESTIMATES 

56 

10-04 

RELIABILITY 

OF 

CORRELATION  RESULTS 

57 

SECTION  11  - TERMS  AND  SYMBOLS 


REFERENCES 

APPENDIX  1 - KEY  QUESTIONS  AND  ANSWERS 
APPENDIX  II  - PROBLEMS  AND  SOLUTIONS 


- « - 


Edit  I ti  ITS 


1.  Illustration  of  Frequency  Array  and  Chronologic  Sequence 

2.  Graphical  Frequency  Curve 

3.  Taoulation  of  Peak-Flow  Frequency  Data 

4.  Analytical  Computation  of  l’eak-Flow  Frequency  Curv-. 

5.  Analytical  Frequency  Computation  and  Adjustment 

6.  Errors  of  Estimated  Values 

7.  Error-Limit  Curves 

3.  Partial  -Duration  Frequency  Curve 

9.  Flow-Duration  Curve 

10.  Uiver  Stage  Frequency  Curve 

11.  Analytical  Frequency  Computation  Using  Pre-record  Data 

12.  Analytical  Frequency  Curve  Using  Pre-record  Data 

13.  Analytical  Frequency  Computation  Omitting  Lower  Flows 

14.  Analytical  Frequency  Curve  Omitting  Lower  Flows 

15.  Maximum-Runoff  Volume  Frequency  Data  - Chronoj ogical  Order 

16.  Maximum-Runof f Volume  Frequency  Data  - Order  of  Magnitude 

17.  Maximum- Runoff  Volume  Frequency  Data  - Logarithms 

1«.  Maximum-Runoff  Volume  Frequency  Computations 

ID.  Standard-Deviation  Adjustment 

20.  Maximum-Runoff  Volume  Frequency  Curves 

21.  Volume-Duration  Curves 

22.  Regional  Frequency  Correlation 

23.  Standard  Lrrors  of  Frequency  Statistics 

24.  .-lap  of  C 

25.  Map  of  Sp 

26.  basin  Elevation  Factor,  K 

27.  Index  Frequency  Curves  of  Peak  Runoff 

28.  Frequency  Curve  of  Hurricane  Pressures 

29.  Frequency  Curve  of  Sediment  Load 

30.  Storage  Frequency  Computation  based  on  Coincident  Frequencies 

31.  Cotaputation  of  Simple  Linear  Correlation 

32.  Illustration  of  Simple  Linear  Correlation 

33.  Computation  of  Multiple  Linear  Correlation 

34.  Values  of  t 

33.  Values  of 

36.  Transformed  Plotting  Positions  (2  sheets) 

37.  Table  of  Plotting  Positions  (2  sheets) 

35.  Table  of  k Versus  P.^ 

39.  Table  of  Peatoon  Type  111  Coordinates 

40.  Table  of  P^  versus  P«, 

41.  Table  of  Logarithms 


1-01 


STATISTICAL  METHODS  IN  HYDROLOGY 
SECTION  1 - INTRODUCTION 

1-01.  NATURE  AND  OBJECTIVES  OF  STATISTICAL  ANALYSIS 

Statistical  analysis  as  applied  in  hydrologic  engineering  consists 
of  (a)  estimating  the  future  frequency  or  probability  of  hydrologic 
events  based  on  information  contained  in  hydrologic  records  and  (b) 
correlating  interrelated  hydrologic  variables.  In  probability  analyses, 
statistical  methods  permit  coordination  of  observed  data  to  yield  a more 
accurate  estimate  of  future  frequencies  than  is  indicated  by  the  raw 
data,  and  also  provide  criteria  for  judging  the  reliability  of  the 
frequency  estimates.  In  correlation  analyses,  statistical  methods 
provide  means  for  deriving  the  most  likely  relationship  between  two 
variables,  and  also  provide  criteria  for  judging  the  reliability  of 
forecasts  or  estimates  based  on  the  derived  relationship. 

1-02.  PURPOSE  AND  SCOPE 

a.  This  is  a revision  of  the  paper  of  the  same  title  dated  July 
1952  and  distributed  to  Corps  of  Engineers  offices  by  Engineer  Bulletin 
52-24,  dated  26  August  1952.  Ihe  revision  incorporates  new  material 
developed  under  the  Civil  Works  Investigations  program  of  the  Corps. 

b.  It  is  the  purpose  of  this  publication  to  describe  and  illustrate 
the  application  of  statistics  in  hydrologic  engineering.  The  subject 
matter  covers  the  following  items: 

(1)  A concise  review  of  the  basic  concepts  of  probability  and 
correlation  analyses  that  are  applicable  in  hydrologic  engineering, 
with  a guide  to  supplemental  reading  for  further  treatment. 

(2)  Presentation  of  detailed  computation  procedures  and 
supporting  justifications  and  computation  aids  for  derivation  of 
probability  of  frequency  estimates  based  on  analysis  of  hydrologic 
records  that  have  been  adjusted  ’s  required  to  conform  with  selected 
reference  base  conditions. 

(3)  A summary  of  procedures  for  developing  "regionalized" 
hydrologic  frequency  estimates,  based  on  analyses  of  hydrologic 
records  available  at  stream  gaging  stations,  adjusted  to  provide 
generalized  flood-frequency  relations  that  are  considered  most 
representative  of  long-period  hydrologic  characteristics  in  various 
drainage  areas  in  the  region.  Also,  illustrations  and  explanations 
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of  simple  generalization  procedures  for  use  where  these  are  adequate 
and  advantageous  are  given. 

(4)  A brief  list  of  key  questions  and  simple  problems 
suitable  for  use  in  training  classes.  More  detailed  training  docu- 
ments have  been  prepared  in  connection  with  training  courses  given 
under  Civil  Works  Investigation  Project  CW-151  in  the  Sacramento 
District  of  the  Corps  of  Engineers. 

(5)  Discussions  pertaining  to  certain  aspects  of  statis- 
tical analyses  associated  with  hydrologic  engineering  that  deserve 
special  emphasis. 

c.  For  those  who  are  interested  only  in  frequency  analyses 
of  flood  discharges,  sections  3 and  4 and  exhibits  2 to  4 contain  the 
essential  guide  material.  For  those  interested  only  in  correlation 
analyses,  section  9 contains  the  essential  material.  The  detailed 
procedure  for  computing  the  frequencies  of  flood  peaks  and  volumes 
is  given  step  by  step  on  exhibit  18. 

1-03.  REFERENCES 

There  are  a great  many  textbooks  on  statistics,  probability, 
and  correlation,  and  a multitude  of  technical  papers  on  the  statis- 
tical aspects  of  hydrologic  engineering.  A few  references  con- 
sidered particularly  applicable  and  useful  for  the  purpose  of 
Supplementing  material  contained  herein  are  given  at  the  end  of 
the  text. 

1-04.  ORGANIZATION  OF  THIS  PAPER 

The  technical  presentation  in  this  paper  begins  with  a brief 
discussion  of  statistical  concepts  and  definitions  contained  in 
section  2.  Graphical  methods  of  frequency  analyses  are  covered  in 
section  3,  and  numerical  or  analytical  procedures  are  contained  in 
section  4.  A recommended  procedure  for  adjusting  frequency  esti- 
mates based  on  short  records,  using  data  at  long  record  stations, 
is  contained  in  section  5.  Extension  of  the  methods  discussed  in 
section  4 for  use  in  estimating  frequencies  of  flood  volumes  is 
contained  in  section  6.  Recommended  means  for  coordinating  fre- 
quency estimates  within  a region  for  the  purposes  of  increasing 
the  reliability  of  estimates  and  for  deriving  estimates  for 
ungaged  areas  are  contained  in  section  7.  A brief  description  of 
the  application  of  frequency  procedures  to  hydrologic  factors 
other  than  runoff  is  contained  in  section  8.  Section  9 contains 
a brief  exposition  of  correlation  methods  generally  in  use.  These 
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are  readily  available  In  many  textbooks,  but  they  are  given  here  for 
convenience  and  for  the  purpose  of  standardizing  notation.  Criteria 
for  evaluating  the  reliability  of  frequency  statistics,  frequency 
estimates  and  correlation  results  are  contained  in  section  10.  Terms 
and  symbols  as  used  herein  are  surrmarized  in  section  11.  Following 
this  is  given  a list  of  selected  references.  Common  questions  and 
answers  intended  to  clarify  some  of  the  more  complex  aspects  of 
probability  analysis  are  contained  in  Appendix  I,  and  illustrative 
problems  for  use  as  exercises  in  connection  with  a frequency  course 
are  given  in  Appendix  II. 
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SECTION  2 - GENERAL  PROBABILITY  CONCEPTS  AND  DEFINITIONS 
2-01.  INTRODUCTION 

The  subjects  of  probability  and  statistics  are  becoming  increas- 
ingly applicable  in  engineering  work,  and  it  is  considered  appropriate 
to  provide  background  information  for  orienting  those  engineers  who 
have  not  had  formal  training  in  the  subjects.  This  section  contains  a 
brief  review  of  the  theoretical  basis  for  probability  estimates  based 
on  observed  data.  Details  of  application  will  be  presented  later,  and 
a broader  understanding  of  the  theory  can  be  obtained  from  textbooks 
such  as  reference  1 . 

2-02.  NATURE  OF  RANDOM  EVENTS 

a.  Probability  estimates  made  in  hydrologic  engineering  are  based 
on  records  of  random  events.  To  understand  probability  methods  and 
fully  appreciate  the  degree  of  reliability  of  such  probability  esti- 
mates, one  should  consider  the  nature  and  variation  of  random  samples. 

b_.  Consider  a period  of  2,000  years  during  which  controlling 
hydrologic  conditions  do  not  change.  Annual  maximum  hydrologic  events 
occurring  during  this  period  can  be  divided  into  100  records  of  20  years 
each.  From  knowledge  of  probability,  it  is  expected  that  one  of  these 
records  will  contain  a flood  that  is  exceeded  on  the  average  only  once 
in  2,000  years,  a very  rare  event.  About  13  of  these  records  will 
contain  floods  that  are  exceeded  on  the  average  only  once  in  100  years 
(it  would  be  20,  except  that  some  of  the  records  might  contain  more 
than  one  of  these  large  floods),  and  64  of  the  records  should  contain 
floods  larger  than  that  exceeded  on  the  average  once  in  20  years.  On 
the  other  hand,  about  12  of  the  records  would  not  have  floods  larger 
than  that  exceeded  on  the  average  once  in  10  years. 

c.  When  a hydrologic  engineer  is  studying  a record  of  20  years' 
length,  he  cannot  tell  by  examining  t.he  record  alone  whether  it  is 
one  that  has  a normal  sequence  of  events,  abnormally  rare  events,  or 
an  abnormally  small  number  of  large  events.  If  the  record  contains 
abnormally  large  events,  the  resulting  probability  estimates  for  large 
events  will  be  too  high,  and  vice  versa.  In  order  to  reduce  the 
uncertainties  from  this  source,  it  is  advisable  to  study  all  of  the 
events  in  relation  to  each  other  and  to  introduce  knowledge  obtained 
on  similar  phenomena  at  other  locations. 
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2-03.  PROBABILITY  INFERENCE 

a.  Knowledge  that  a certain  set  of  conditions  can  result  in 
various  sets  of  data  because  of  random  variations  is  used  to  Infer 
that  any  particular  set  of  data  could  have  resulted  from  various  sets 
of  hydrologic  conditions.  The  problem  in  making  probability  estimates 
is  essentially  to  determine  the  set  of  conditions  that  most  likely 
generated  the  sample  of  data  that  has  been  recorded.  This  set  of 
conditions  is  represented  by  a "parent  population"  which  consists  of 
all  of  the  hydrologic  events  that  v/ould  be  generated  if  the  record 
continues  indefinitely  and  controlling  conditions  do  not  change. 

b.  In  probability  analysis,  there  are  two  basic  approaches  to 
estimating  or  Inferring  the  parent  population  from  sample  data.  First, 
data  can  be  arranged  in  the  order  of  magnitude  to  form  a frequency 
array,  illustrated  on  exhibit  1,  and  a graph  of  magnitude  versus 
observed  frequency  plotted,  as  shown  on  exhibit  2.  A smooth  curve 
drawn  through  the  plotted  data  would  represent  an  estimate  of  the 
parent  population  from  which  the  probabilities  of  future  events  can 

be  determined.  The  second  basic  approach  is  to  derive  from  the  data 
general  statistics  representing  the  average  magnitude  of  floods,  the 
variability  from  that  average,  and  any  other  pertinent  statistics 
relating  frequencey  to  magnitude  as  indicated  by  the  data. 

2-04.  CUMULATIVE  FREQUENCY  CURVES 

a.  A cumulative  frequency  curve,  or  simply  frequency  curve,  such 
as  that  illustrated  on  exhibit  2,  relates  the  magnitude  of  an  event 
to  the  frequency  with  which  that  magnitude  is  exceeded  as  events  occur 
at  random.  For  example,  if  25  floods  at  a location  exceed  10,000  c.f.s. 
in  100  years,  on  the  average,  then  the  value  of  10,000  c.f.s.  on  the 
cumulative  frequency  curve  will  correspond  to  an  exceedence  probability 
of  0.25  in  any  1 year  or  an  exceedence  Frequency  of  25  times  per  100 
years,  250  times  per  thousand  years,  etc.,  or  simply  25  percent.  While 
a single  frequency  curve  can  represent  the  frequency  of  peak  discharges 
at  a qiven  location,  the  frequency  of  flood  volumes  would  be  represented 
by  an  entirely  different  curve,  or  by  more  than  one  curve  if  volumes 
for  various  durations  are  concerned.  Frequency  curves  can  also  be 
used  to  represent  the  frequency  of  reservoir  stages,  river  stages, 
precipitation,  and  many  other  phenomena. 

b_.  Frequency  curves  are  most  commonly  used  in  flood  control 
benefits  studies  for  the  purpose  of  evaluating  the  economic  effect 
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of  the  project.  Other  common  uses  of  frequency  curves  Include  the 
determination  of  reservoir  stage  frequency  for  real  estate  acquisi- 
tion and  reservoir-use  purposes,  the  selection  of  rainfall  frequency 
for  storm-drain  design,  and  the  selection  of  runoff  frequencies  for 
interior  drainage,  pumping  plant,  and  local-protection  project  design. 

c.  The  basic  frequency  curve  used  in  hydrologic  engineering  is 
the  frequency  curve  of  annual  maximum  or  annual  minimum  events.  A 
second  curve,  the  partial-duration  curve,  represents  the  frequency 
of  all  events  above  a given  base  value,  regardless  of  whether  two 

or  more  occurred  in  the  same  year.  Either  curve  must  be  supplemented 
by  considerations  of  seasonal  effects  and  other  factors  in  applica- 
tion, as  explained  in  reference  5.  When  both  the  frequency  curve  of 
annual  floods  and  the  partial-duration  curve  are  prepared,  care  must 
be  exercised  to  assure  that  the  two  are  consistent.  Normal  relation- 
ships between  the  two  are  given  in  paragraph  4-04. 

d.  In  almost  all  locations  there  are  seasons  during  which  storms 
or  floods  do  not  occur  or  are  not  severe,  and  other  seasons  when  they 
are  more  severe.  Also,  damages  associated  with  a flood  often  vary 

with  season  of  the  year,  among  other  factors.  In  many  types  of  studies, 
the  seasonal  variation  factor  is  of  primary  importance,  and  it  becomes 
necessary  to  establish  frequency  curves  for  each  month  or  other  sub- 
division of  the  year.  For  example,  one  frequency  curve  might  represent 
the  largest  floods  that  occur  each  January,  a second  one  would  represent 
the  largest  floods  that  occur  each  February,  etc.  In  another  case,  one 
frequency  curve  might  represent  floods  during  the  snowmelt  season,  while 
a second  might  represent  floods  during  the  rain  season.  Occasionally, 
when  seasons  are  studied  separately,  an  annual-event  curve  covering  all 
seasons  is  also  prepared,  and  care  should  be  exercised  to  assure  that 
the  various  seasonal  curves  are  consistent  with  the  annual  curve 
(reference  19). 

e.  In  connection  with  power  studies  for  run-of-river  plants 
particularly,  and  in  some  phases  of  sediment  studies,  the  flow- 
duration  curve  serves  a useful  purpose.  It  simply  represents  the 
percent  of  time  during  which  specified  flow  rates  are  exceeded  at 
a given  location.  Ordinarily,  variations  within  periods  less  than 
one  day  are  not  of  consequence,  and  the  curves  are  therefore  based 

on  observed  mean-daily  flows.  For  the  purpose  served  by  flow-duration- 
curves,  the  extreme  rates  of  flow  are  not  important,  and  consequently 
there  is  no  need  for  refining  the  curve  in  regions  of  high  flow.  The 
procedure  ordinarily  used  in  the  preparation  of  a flow-duration  curve 
consists  of  counting  the  number  of  mean-daily  flows  that  occur  within 
given  ranges  of  magnitude.  Then  the  lower  limit  of  magnitude  in  each 
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range  is  plotted  aqainst  the  percentage  of  days  of  record  that  mean- 
daily  flows  exceed  that  magnitude.  A typical  flow-duration  curve  is 
shown  on  exhibit  9. 

2-05.  TECHNICAL  APPROACH  IN  ESTIMATING  FREQUENCY  CURVES 

a.  There  are  two  basic  approaches  to  estimating  frequency  curves— 
graphical  and  analytical.  Each  of  these  approaches  has  several 
variations  in  current  practice,  but  the  discussion  herein  will  be 
limited  to  selected  methods. 

b.  Graphically,  frequencies  are  evaluated  simply  by  arranging 
observed  values  in  the  order  of  magnitude  and  considering  that  a 
smooth  curve  suggest  d by  that  array  of  values  ir,  representative  of 
future  possibilities.  Each  value  represents  a fraction  of  the  future 
possibilities  and,  when  plotting  the  frequency  curve,  it  is  given  a 
"plotting  position"  that  is  calculated  to  give  it  the  proper  weight 
(see  paragraph  3-05). 

c.  In  the  application  of  analytical  (statistical)  procedures, 
the  concept  of  theoretical  populations  or  distributions  is  employed, 
as  discussed  in  paragraph  2-03.  A distribution  is  a set  of  values 
that  would  occur  under  fixed  conditions  in  an  infinite  amount  of 
time.  Those  that  have  occurred  are  presumed  to  constitute  a random 
sample  and  accordingly  are  used  to  make  particular  inferences  regard- 
ing their  "parent  population"  (i.e.,  the  distribution  from  which  they 
were  derived).  Such  Inferences  are  necessarily  attended  by  considerable 
uncertainty,  because  a given  set  of  observations  could  result  from  any 
of  many  sets  of  physical  conditions  (from  any  one  of  many  distributions). 
However,  by  the  use  of  statistical  processes,  the  most  probable  nature 
of  the  distribution  from  which  the  data  were  derived  can  be  estimated. 
Since  in  all  probability  this  is  not  the  true  parent  population,  the 
relative  chance  that  variations  from  this  "maximum  likelihood" 
distribution  might  be  true  must  be  evaluated.  Each  range  of  possible 
parent  population  is  then  weighted  in  proportion  to  its  likelihood  to 
obtain  a weighted  average  as  demonstrated  in  reference  4.  A probability 
obtained  from  this  weighted  average  is  herein  referred  to  as  the 
expected  probability,  P^.  Computation  procedures  are  given  in  section  4. 

d.  Because  of  the  shortness  of  hydrologic  records,  frequency 
determinations  are  relatively  unreliable  where  based  on  a single  record 
(see  paragraph  10-03).  Also,  it  is  often  necessary  to  estimate 
frequencies  for  locations  where  no  record  exists.  For  these  reasons, 
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regionalized  frequency  studies.  In  which  frequency  characteristics  are 
related  to  drainage-basin  features,  are  ceslrable.  These  are  facilitated 
by  the  use  of  analytical  methods,  as  illustrated  in  section  7. 

2-06.  TERMS  AND  SYMBOLS 

Special  terms  and  symbols  are  defined  where  first  used  herein.  A 
summary  is  given  in  section  11  for  convenience. 
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SECTION  3 - FLOOD  PEAK  FREQUENCY  - GRAPHICAL  METHOD 
3-01.  USE  AND  RELATIVE  ADVANTAGES 

Every  frequency  study  should  be  plotted  graphically,  even  though 
the  results  can  be  obtained  entirely  analytically  as  described  in 
section  4,  in  order  that  observed  data  may  be  visually  compared  with 
the  derived  curve.  The  graphical  method  of  frequency-curve  determin- 
ation can  be  used  for  any  type  of  frequency  study,  but  analytical 
methods  have  certain  advantages  where  they  are  applicable  (see  para- 
graph 4-01).  The  principal  advantages  of  graphical  methods  are  that 
they  are  generally  applicable,  that  the  derived  curve  can  be  easily 
visualized,  and  that  the  observed  data  can  be  readily  compared  with 
the  computed  results.  However,  graphical  methods  of  frequency  analysis 
are  inferior  In  accuracy  to  analytical  methods  where  the  latter  apply, 
and  do  not  provide  means  of  evaluating  the  reliability  of  the  estimates. 
Comparison  of  the  adopted  curve  with  plotted  points  is  not  an  index  of 
reliability  as  in  correlation  analysis,  but  it  is  often  erroneously 
assumed  to  be,  thus  implying  a much  greater  reliability  than  is 
actually  attained.  For  these  reasons,  graphical  methods  should  be 
limited  to  those  cases  where  analytical  methods  do  not  apply  (that  is, 
where  frequency  curves  are  too  irregular  to  compute  analytically)  and 
to  use  as  a visual  aid  or  check  on  analytical  computations. 

3-02.  GENERAL  PROCEDURE 

a.  Graphical  construction  of  a frequency  curve  simply  consists 
of  arranging  the  selected  data  in  the  order  of  magnitude  and  plotting 
the  magnitude  of  each  item  on  the  vertical  scale  against  its  estimated 
exceedence  frequency  (plotting  position  - see  paragraph  3-05)  on  the 
horizontal  scale,  using  a suitable  grid  (see  paragraph  3-06).  A 
smooth  curve  drawn  through  the  points  is  the  desired  frequency  curve. 

jj.  Data  used  in  the  construction  of  frequency  curves  of  peak  flows 
consist  of  the  maximum  flow  for  each  year  of  record  and  all  of  the 
secondary  flows  that  exceed  a selected  base  value.  This  base  value 
must  be  smaller  than  any  floodflow  that  is  of  importance  in  the  analysis, 
and  should  also  be  low  enough  so  that  the  total  number  of  floods  in 
excess  of  the  base  equals  or  exceeds  the  number  of  years  of  record. 
Ordinarily,  the  latter  criterion  controls,  and  the  two  series  of  events 
tabulated  are  equal  in  number.  Exhibit  3 is  a sample  tabulation  of  data 
directly  from  the  record,  of  the  frequency  arrays  in  the  order  of 
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magnitude,  and  of  corresponding  plotting  positions.  In  arranging  the 
data  in  order  of  magnitude,  much  time  can  be  saved  by  taking  the  events 
in  the  order  in  which  they  occurred,  and  placing  them  rapidly  on  a 
blank  sheet  of  paper  in  the  order  of  magnitude,  thus  making  an  irregular 
tabulation.  This  tabulation  is  then  recopied  onto  the  form  shown  on 
exhibit  3 and  cross-checked  with  chronologic  values. 

3-03.  SELECTION  AND  ARRANGEMENT  OF  DATA 

a.  The  primary  consideration  in  selection  of  an  array  of  data  for 
a frequency  study  is  the  use  to  which  the  frequency  estimates  will  be 
put.  If  the  frequency  curve  is  to  be  used  for  estimating  damages  that 
are  related  to  instantaneous  peak  flows  in  a stream,  peak  flows  should 
be  selected  from  the  record.  If  the  damages  are  related  to  maximum 
mean-daily  flows  or  to  maximum  3-day  flows,  these  items  should  be 
selected.  If  the  behavior  of  a reservoir  under  investigation  is  related 
to  the  3-day  or  10-day  rainflood  volume,  or  to  the  seasonal  snowmelt 
volume,  that  pertinent  item  should  be  selected.  Occasionally,  it  is 
necessary  to  select  a related  variable  in  lieu  of  the  one  desired.  For 
example,  where  mean-daily  flow  records  are  more  complete  than  the 
records  of  peak  flows,  it  may  be  more  desirable  to  derive  a frequency 
curve  of  mean-daily  flows  and  then,  from  the  computed  curve,  derive  a 
peak-flow  curve  by  means  of  an  empirical  relation  between  mean-daily 
flows  and  peak  flows.  All  reasonably  independent  values  should  be 
selected,  but  the  annual  maximum  events  should  ordinarily  be  segregated 
when  the  application  of  analytical  procedures  discussed  in  section  4 is 
contemplated. 

b.  Data  selected  for  a frequency  study  must  measure  the  same  aspect 
of  each  event  (such  as  peak  flow,  mean-daily  flow,  or  flood  volume  for  a 
specified  duration),  and  each  event  must  be  controlled  by  a uniform  set 
of  hydrologic  and  operational  factors.  For  example,  it  would  be  improper 
to  combine  items  from  old  records  that  are  reported  as  peak  flows  but  are 
in  fact  only  daily  readings,  with  newer  records  where  the  peak  was 
actually  measured.  Similarly,  care  should  be  exercised  when  there  has 
been  significant  change  in  upstream  storage  regulation  during  the  period 
of  record  so  as  not  inadvertently  to  combine  unlike  events  into  a single 
series.  In  such  a case,  the  entire  record  should  be  adjusted  to  a 
standard  condition, 

c.  Hydrologic  factors  and  relationships  operating  during  a general 
winter  rainflood  are  usually  quite  different  from  those  operating  during 
a spring  snowmelt  flood  or  during  a local  summer  cloudburst  flood, 

Where  two  or  more  types  of  floods  are  distinct  and  do  not  occur  pre- 
dominantly in  mutual  combinations,  they  should  not  be  combined  into 
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a single  series  for  frequency  analysis.  It  Is  usually  more  reliable 
in  such  cases  to  segregate  the  data  in  accordance  with  type  and  to 
combine  only  the  final  curves,  If  necessary.  In  the  Sierra  Nevada 
region  of  California  and  Nevada,  frequency  studies  are  made  separately 
for  rainfloods,  which  occur  principally  during  the  months  of  November 
through  March,  and  for  snowmelt  floods,  which  occur  during  the  months 
of  April  through  July.  Flows  for  each  of  these  two  seasons  are 
segregated  strictly  by  cause  - those  predominantly  caused  by  snow- 
melt and  those  predominantly  caused  by  rain.  In  desert  regions, 
summer  thunderstorms  should  be  excluded  from  frequency  studies  of 
winter  rainfloods  or  spring  snowmelt  floods  and  should  be  considered 
separately. 

d.  Occasionally  a runoff  record  may  be  interrupted  by  a period  of 
one  or  more  years.  If  the  interruption  is  caused  by  destruction  of  the 
gaging  station  by  a large  flood,  failure  to  fill  in  the  record  for  that 
flood  would  have  a biasing  effect,  which  should  be  avoided.  However, 
if  the  cause  of  the  interruption  is  known  to  be  independent  of  flow 
magnitude,  the  entire  period  of  interruption  should  be  eliminated  from 
the  frequency  array,  since  no  bias  would  result.  Knowledge  we  have 
about  floods  observed  at  other  locations  during  periods  of  no  record 
at  the  site  concerned  can  be  utilized  as  discussed  in  section  5.  In 
cases  where  no  runoff  records  are  available  on  the  stream  concerned, 
it  is  usually  best  to  estimate  the  frequency  curve  as  a whole  using 
regional  generalizations  discussed  in  section  7,  instead  of  attempting 
to  estimate  a complete  series  of  individual  floods,  because  ordinary 
methods  of  estimating  individual  floods  tend  to  reduce  the  slope  of 
the  frequency  curve  (the  standard  deviation). 

3-04.  ADJUSTMENT  TO  UNIFORM  CONDITION 

Since  the  frequency  analysis  of  hydrologic  data  is  based  on  the 
assumption  of  random  occurrences,  each  item  of  data  must  have  occurred 
under  similar  hydrologic  conditions  or  must  be  adjusted  to  a standard 
uniform  condition.  If  control  by  reservoirs  or  diversion  for  irri- 
gation, etc.,  has  affected  the  runoff,  some  adjustment  of  the  data 
is  usually  necessary.  Where  it  is  feasible  to  adjust  to  natural 
conditions,  it  is  advisable  to  do  so  in  order  that  the  data  will  more 
nearly  conform  to  theoretical  frequency  functions  that  have  been 
found  to  describe  the  frequency  of  natural  hydrologic  events.  This 
is  accomplished  by  standard  routing  procedures,  and  in  many  cases 
an  approximate  adjustment  is  satisfactory.  Where  the  regulation 
is  complex,  as  in  the  case  of  a large  number  of  upstream  reservoirs, 
it  may  be  advisable  to  adjust  the  data  to  a uniform  condition  with 
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specific  reservoirs  and  diversion  facilities  operating.  For  design 
purposes,  a frequency  curve  of  runoff  under  '‘non-project"  conditions 
that  is  expected  to  prevail  during  the  lifetime  of  the  proposed  project, 
if  the  project  is  not  constructed,  is  required.  A frequency  curve  based 
on  any  specified  uniform  condition  can  be  converted  to  one  for  nonproject 
conditions  using  relationships  developed  by  routing  "balanced"  floods  of 
specified  frequency  ( 1 . e . , floods  having  runoff  for  various  durations  and 
in  various  portions  of  the  drainage  basin  of  equal  exceedence  frequency). 
Techniques  for  doing  this  are  outside  the  scope  of  this  paper. 

3-05.  PLOTTING  FORMULA 

a.  The  reasoning  behind  the  selection  of  an  exact  formula  for 
plotting  the  frequency  of  observed  flood  events  is  extremely  complex. 
Approximate  plotting  positions  can  be  obtained  by  reasoning  that  each 
item  in  a set  of,  say  10,  represents  10  percent  of  the  parent-population 
events  and  should  be  plotted  in  the  middle  of  its  group,  that  is  at  5, 

15,  25  percent,  etc.,  for  successive  events  in  the  order  of  magnitude. 

The  formula  derived  from  this  line  of  reasoning  has  been  used  in  the 
past,  and  is  generally  satisfactory,  considering  the  overall  reliability 
of  the  results.  However,  more  accurate  plotting  positions  have  been 
derived  theoretically  and,  since  their  use  is  very  simple,  it  is 
considered  advantageous  to  use  them.  Plotting  positions  recommended 
are  shown  on  exhibit  37,  and  are  based  on  the  premise  that  if  they  are 
used  repeatedly  in  a great  number  of  random  samples,  they  will  prove  to 
be  too  low  in  half  of  the  cases  and  too  high  in  the  other  half,  compared 
with  the  theoretically  true  values  that  cannot  be  determined  because  of 
random  variations  in  data.  They  are,  therefore,  called  median  plotting 
positions.  There  are  other  systems  of  deriving  plotting  positions  that 
yield  good  results,  but  of  the  formulas  generally  used,  use  of  the 
median  plotting  position  will  most  nearly  duplicate  results  obtained 
by  analytical  methods  of  frequency  analysis  that  do  not  require  plot- 
ting positions. 
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b.  In  ordinary  hydrologic  frequency  work*  exceedence  frequencies 
are  expressed  in  percent  or  in  terms  of  events  per  hundred  years,  as 
shown  on  exhibit  37,  which  give  plotting  positions  for  arrays  up  to 
100  events  in  size.  For  arrays  larger  than  100,  the  plotting  position, 

P,  can  be  obtained  as  were  those  of  exhibit  37  by  use  of  the  following 
equation: 

1 - P1  - (0.5)1/N  (1) 

in  which  P-j  is  the  plotting  position  for  the  largest  event,  and  N is 
the  number  of  years  of  record.  The  plotting  position  for  the  smallest 
event  Is  the  complement  of  this  value,  and  all  other  plotting  positions 
are  interpolated  linearly  between  these  two.  For  partial-duration  curves, 
particularly  where  there  are  more  events  than  years  (N) , plotting  positions 
larger  than  50  percent  are  obtained  by  use  of  the  following  equation: 

P * (2m  - 1)/2K  (la) 

in  which  m is  the  order  number  of  the  event. 

3-06.  PLOTTING  GRID 

If  hydrologic  frequency  data  are  plotted  with  Cartesian  coordinates, 
the  resulting  frequency  relationship  will  curve  rather  abruptly  at  the 
upper  end  and  possibly  at  the  lower  end  also.  Furthermore,  the  extreme 
values  in  which  there  is  the  greatest  interest  would  be  compressed  into 
a very  small  area,  and  extrapolation  of  the  curve  would  be  difficult. 
Accordingly,  it  has  been  found  desirable  to  use  a plotting  grid  on 
which  a frequency  curve  of  hydrologic  data  will  usually  ^proximate  a 
straight  line.  A grid  that  has  been  found  to  be  suitable  for  this 
purpose  is  the  probability  grid.  This  grid  is  designed  so  that  the 
cumulative  frequency  curve  of  a variable  that  is  distributed  in  accord- 
ance with  the  normal  probability  curve  will  plot  as  a straight  line. 

The  grid  is  illustrated  on  exhibit  10.  It  has  been  found  that  items 
such  as  air  temperature  or  river  stage  that  either  do  not  have  a fixed 
lower  limit  of  zero  or  whose  lower  limit  is  far  removed  from  the  range 
of  experience,  will  often  yield  frequency  curves  approximating  a straight 
line  when  plotted  on  this  grid.  Variables  such  as  streamflow  where  a 
lower  limit  of  zero  is  often  approached  in  experience  will  ordinarily 
yield  an  approximately  straight  frequency  curve  only  if  the  logarithms 
are  plotted  on  this  grid.  For  convenience,  the  logarithmic  probability 
grid  illustrated  on  exhibit  2 has  been  devised  so  that  flows  can  be 
plotted  directly  to  yield  an  approximately  straight  line. 
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3-07.  SAMPLE  COMPUTATION 

a.  Exhibit  2 shows  the  plotting  of  a frequency  curve  of  annual 
peak"‘f1ows  corresponding  to  those  tabulated  on  exhibit  3.  The  curve 
should  be  drawn  so  as  to  balance  out  the  plotted  points  to  a reason- 
able degree  and  so  that  there  is  no  abrupt  break  in  the  frequency  curve. 
Unless  computed  as  discussed  in  section  4,  the  frequency  curve  should 

be  drawn  as  a straight  line  on  the  grid  whenever  data  reasonably  indicate 
a straight  line  and  conditions  do  not  exist  that  would  make  a straight 
line  unreasonable. 

b.  The  partial-duration  curve  corresponding  to  the  partial- 
duration  data  on  exhibit  3 has  been  shown  on  exhibit  8.  This  curve 
has  been  drawn  by  generally  balancing  out  the  plotted  points  except 
that  it  was  made  to  conform  with  the  annual -event  curve  in  the  upper 
portion  in  general  accord  with  the  standard  relationship  discussed  in 
paragraph  4-04.  './hen  partial-duration  data  must  include  more  events 
than  there  are  years  of  record  (see  paragraph  3-02)  it  will  be  necessary 
to  use  logarithmic  paper  for  plotting  purposes,  as  on  exhibit  8,  in 
order  to  plot  exceedence  frequencies  greater  than  100  percent.  Otherwise, 
the  curve  can  be  plotted  on  probability  grid,  as  illustrated  on  exhibit  20. 

c.  Exhibit  10  illustrates  the  graphical  construction  of  a river-stage 
frequency  curve.  The  shape  of  this  curve  is  dictated  by  the  plotted  points 
and,  in  some  cases,  by  consideration  of  the  stage  at  which  overbank  flows 
begin.  Whenever  stage  is  a consistent  function  of  flow,  as  is  the  usual 
case,  the  stage-frequency  curve  should  be  obtained  from  the  flow- 
freouency  and  stage-discharge  curves. 
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SECTION  4 - FLOOD  PEAK  FREQUENCY  - ANALYTICAL  COMPUTATION 
4-01.  USE,  LIMITATIONS  AND  RELATIVE  ADVANTAGES 

The  analytical  method  of  computing  a frequency  curve  In  hydrologic 
engineering  is  limited  almost  exclusively  to  curves  of  annual  maximum 
or  annual  minimum  streamflows  for  a specified  duration  (including  peak 
flows)  and  annual  maximum  precipitation  amounts  for  a specified  duration. 
In  general,  the  results  obtained  by  analytical  methods  are  considerably 
more  reliable  than  those  obtained  by  graphical  procedures.  They  have 
the  additional  advantages  that  the  degree  of  reliability  of  frequency 
estimates  can  be  evaluated,  as  discussed  in  paragraph  10-03. 

4-02.  EQUATIONS  USED 

a..  Frequency  curves  are  computed  analytically  by  the  use  of 
moments  of  the  logarithms,  exprc  .d  in  terms  of  the  mean,  M,  (first 
moment),  standard  deviation,  S,  (second  moment)  and  skew  coefficient, 
g,  (third-moment  function).  The  three  corresponding  equations  used 
are  as  follows: 


M = eX/N 


_ EX2  _ EX2  - (EX)2/N 

S - r-T iTi~ 

g = N Ex3 = N2EX3  - 3N  EX  EX2  + 2(eX)3 

(N-l ) (N-2)S3  N(N-1)(N-2)S3 


in  which: 

X = Magnitude  of  an  event  (logarithm) 
x = X - M,  deviation  of  a single  event  from  the  mean 
N = Number  of  events  in  the  record 

b.  The  types  of  cumulative  frequency  curves  fitted  in  hydrologic 
engineering  do  not  require  moments  of  a higher  order  than  these  three, 
and  ordinary  fitting  will  require  only  the  first  two.  There  is  no 
need  to  tabulate  the  individual  deviation  for  each  logarithm,  as  the 
second  parts  of  equations  3 and  4 can  be  used  to  compute  the  standard 
deviation  and  skew  coefficient  much  more  rapidly,  directly  from  the 
original  logarithms.  In  computing  these  quantities  in  this  manner, 
however,  it  is  essential  that  the  intermediate  quantities  in  the 
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computation  be  accurate  to  at  least  four  decimal  places,  which  is  not 
practicable  without  an  automatic  desk  calculator  or  electronic  computer. 

c.  Computation  of  the  maximum  likelihood  frequency  curve  coordinates 
from”the  computed  mean  and  standard  deviation,  using  the  adopted  skew 
coefficient  (usually  zero  for  peak  discharges  and  for  rainfall  frequencies) 
is  accomplished  as  described  below  by  use  of  the  following  equation: 

Log  Q * M + kS  (5) 

4-03.  ANNUAL-EVEI1T  CURVES 

Frequency  curves  of  annual  maximum  or  minimum  events  are  computed  as 
follows  (See  exhibit  4 for  example): 

a.  Mean  Logarithm.  After  tabulation  of  the  data  in  chronological 
order  orTrPthe  order  of  magnitude,  the  logarithm  (exhibit  41)  of  each 
discharge  is  tabulated  to  two  decimal  places.  The  mean  logarithm  is 
obtained  by  dividing  the  sum  of  these  logarithms  by  the  number  of  events 
(equation  2).  Time  can  be  saved  by  obtaining  the  sum  of  the  logarithms 
on  one  register  of  a calculator  at  the  same  time  that  the  sum  of  the 
squares  of  the  logarithms  are  obtained  on  a second  register  for  step  (b). 

b.  Standard  Deviation.  The  standard  deviation  is  computed 
(equation  3)  as  follows: 

(1)  Obtain  the  sum  of  the  squares  of  the  logarithms  in  an 
automatic  calculator.  This  quantity  should  not  be  rounded  off,  but 
all  figures  carried  in  the  computation. 

(2)  The  sum  of  the  logarithms  obtained  in  the  same  machine 
operation,  which  figure  is  also  not  to  be  rounded  off,  is  squared  and 
divided  by  the  number  of  events.  This  is  a single  machine  operation, 
and  the  quotient  should  be  carried  to  as  many  places  as  is  the  sum  of 
the  squares. 

(3)  This  quotient  is  subtracted  from  the  sum  of  the  squares 
to  obtain  a quantity  numerically  equal  to  the  sum  of  the  squares  of 
the  deviations  from  the  mean. 

(4)  Divide  this  quantity  by  one  less  than  the  number  of 
events.  The  square  root  of  the  quotient  is  the  standard  deviation. 
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When  an  automatic  calculator  is  not  available,  steps  (1),  (2)  and 
(3)  can  be  replaced,  as  illustrated  on  exhibit  4,  by  the  following: 


(1.1)  Tabulate  to  two  decimal  places  the  difference 
between  each  logarithm  and  the  mean  logarithm.  This  quantity  is 
called  the  deviation. 

(2.1)  Tabulate  the  square  of  each  deviation  to  three 
decimal  places. 

(3.1)  Add  the  squares  of  the  deviations. 

c.  Skew  Coefficient.  It  is  impractical  to  base  the  skew 
coefficient  to  be  used  in  a frequency  study  on  a single  record  of 
annual  flows  that  is  less  than  100  years  in  length.  Even  if  such 
a long  record  is  available,  it  is  possible  that  a more  accurate 
determination  of  skew  coefficient  can  be  obtained  by  combining 

the  information  from  other  records.  Unless  the  data  show  a radical 
departure  from  usual  values  of  skew,  zero  skew  coefficient  should 
be  used  for  frequency  curves  of  annual  maximum  peak  flows  or  precip- 
itation, and  coefficients  given  in  paragraph  6-03  should  be  used  for 
frequency  curves  of  annual  maximum  flood  volumes.  Any  radical 
departure  of  the  observed  data  from  the  adopted  skew  coefficients 
would  appear  when  the  curve  and  data  are  plotted  graphically. 

Special  regional  determination  of  skew  coefficients  are  discussed 
in  paragraph  7-11. 

d.  Computation  of  Curve.  Each  frequency  curve  can  be  computed 
as  follows  (See  exhibit  7 for  example): 

(1)  For  selected  values  of  P*,,  tabulate  values  of  k obtained 
from  exhibit  39  corresponding  to  the  adopted  skew  coefficient. 

(2)  Multiple  each  of  these  by  the  computed  standard 
deviation,  and  add  each  product  in  turn  to  the  mean  logarithm 
(Equation  5). 
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(3)  Tabulate  values  of  P,,  from  exhibit  40  corresponding  to 

the  selected  values  1 

(4)  Plot  the  anti  logarithm  of  each  of  the  sums  obtained  in 
step  2 against  the  selected  P^  values  obtained  in  step  3. 

4-04.  PARTIAL -DU RAT  I ON  CURVES 

Once  a frequency  curve  of  annual  events  has  been  established,  the 
corresponding  partial-duration  curve  can  be  determined  analytically 
by  use  of  average  criteria  derived  by  Walter  Langbein.  These  criteria 
are  based  on  the  asst  iption  that  there  are  a large  number  of  flood 
events  each  year  and  that  these  events  are  mutually  independent. 

They  should  not  be  used  without  checking  their  applicability  unless 
only  very  approximate  results  are  desired  and  time  does  not  permit  a 
more  accurate  determination.  An  average  relationship  developed 
empirically  from  many  stations  in  a region  would  ordinarily  be 
preferred,  because  experience  indicates  that  the  observed  relation- 
ship is  often  different  from  the  theoretical  relationship  (as  demon- 
strated in  reference  20).  A suimiary  of  the  Langbein  criteria  is 
contained  in  the  following  tabulation,  and  an  example  of  its  use  is 
shown  on  exhibit  20. 


Corresponding  Exceedence  Frequencies  per  Hundred  Years 


Annual -event  curve 
(Ho.  of  years  flow  is 
exceeded  per  hundred  years) 


1.00 

2.00 

5.0 

10.0 

20 

30 

40 

50 

60 

63.2 

70 

80 

90 

95 


Partial-duration  curve 
(No.  of  times  flow  is 
exceeded  per  hundred  years) 


1.00 

2.02 

5.1 

10.5 

22.3 

35.6 
51.0 

69.3 

91.7 
100 
120 
161 
230 
300 
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4-05.  USE  OF  HISTORICAL  DATA 

Where  one  or  more  large  historical  floods  have  been  estimated, 
such  items  can  be  used  as  a guide  in  drawing  the  frequency  curve, 
particularly  in  the  range  of  higher  flows.  Such  adjustment  would 
ordinarily  be  done  graphically,  and  the  historical  flows  would  be 
plotted  in  the  order  of  magnitude  in  which  they  are  known  to  have 
occurred,  using  the  entire  period  of  history  (not  only  that  period 
dating  from  the  earliest  known  flood)  for  the  computation  of 
plotting  positions.  While  it  is  important  that  all  known  informa- 
tion be  used  in  the  construction  of  a frequency  curve,  it  should 
be  recognized  that  historical  estimates  are  not  as  valuable  as 
comparable  recorded  flows,  and  that  the  largest  knov/n  floods  are 
not  always  representative  for  the  period.  If  historical  flows 
are  particularly  outstanding  relative  to  recorded  data,  procedures 
illustrated  on  exhibits  11  and  12  can  be  used  to  compute  a com- 
posite frequency  curve.  This  consists  of  selecting  plotting  posi- 
tions based  on  the  flood  magnitudes  and  periods  of  record  and  of 
historical  knowledge,  converting  these  plotting  positions  to 
linear  distances  as  measured  on  probability  grid  (k  values  on 
exhibit  36),  and  solving  for  a best-fit  slope  by  use  of  equation 
36,  This  techique  is  explained  fully  in  reference  23. 

4-06.  USE  OF  FLOW  ESTIMATES 

As  discussed  in  paragraph  3-03b,  use  of  a large  number  of 
flow  estimates  based  on  a record  at  a nearby  location  might  lead 
to  erroneous  frequency  estimates.  However,  there  are  some  cases 
where  one  or  a few  estimates  are  essential.  In  the  case  where  a 
record  is  not  obtained  because  the  gage  was  washed  out,  it  is 
imperative  that  some  estimate  of  the  value  be  used.  The  fact 
that  the  estimate  may  be  in  error  by  25  percent  is  minor  compared 
to  the  error  introduced  by  omitting  the  value  from  the  record. 
Also,  where  a comprehensive  flood  volume -duration  series  is  being 
studied  (see  section  6),  and  a few  of  the  items  are  missing,  it 
is  ordinarily  advantageTJUS>  "LU-estimate  these  items  rather  than 
to  have  a different  number  of  items  for  each  duration  studied. 

4-07.  ABNORMAL  DRY-YEAR  EFFECTS 

The  shape  of  the  frequency  curve  is  sometimes  seriously  dis- 
torted by  the  dominance  of  minor  runoff  factors  during  dry  years. 
This  is  particularly  true  where  (a)  floods  are  normally  caused  by 
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rainstorms,  yet  high  base  flows  from  underground  sources  or  occasional 
snowmelt  prevent  true  measures  of  rainflood  runoff  during  very  dry 
years,  or  (b)  floods  are  normally  caused  by  rainstorms,  but  excessive 
diversion  or  channel  losses  deplete  flows  during  dry  years  so  that  the 
frequency^curve  drops  sharply  at  the  lower  end.  In  such  cases,  it  may 
be  advantageous  to  fit  only  the  upper  half  of  the  annual  floods  with 
a theoretical  frequency  curve.  A suggested  procedure,  illustrated  on 
exhibits  13  and  14,  simply  converts  plotting  positions  to  a linear 
distance  as  measured  on  probability  grid  (k  values  on  exhibit  36)  and 
solves  for  a bestfit  slope  (standard  deviation)  by  use  of  equation  36. 

4-08.  USE  OF  SYNTHETIC  FLOODS 

a.  Frequency  estimates  at  best  are  not  fully  reliable.  Although 
theoretical  devices  are  employed  to  extrapolate  frequency  curves 
beyond  experienced  values,  such  extrapolation  is  highly  dangerous. 

It  is  sometimes  possible  to  use  synthetic  floods  for  constructing  a 
frequency  curve  or  especially  for  extrapolating  one  more  accurately. 

b_.  Where  runoff  records  are  not  available,  an  attempt  is  some- 
times made  to  compute  floods  that  would  result  from  various  rainfall 
amounts,  and  then  to  construct  a synthetic  flood  frequency  curve, 
using  the  rainfall  frequency  curve  as  a guide.  This  method  is 
considered  satisfactory  for  airport  drainage  design  and  for  urban 
storm  drain  design,  but  is  ordinarily  not  satisfactory  where  infil- 
tration losses  are  a considerable  percentage  of  precipitation.  In 
these  latter  cases,  it  Is  best  to  construct  synthetic  frequency  curves 
from  regional  correlation  studies  as  discussed  in  section  7. 

c.  A large  hypothetical  flood  can  sometimes  be  used  as  a guide 
in  extrapolating  the  frequency  curve.  Where  it  appears  that  major 
storms  may  have  accidentally  missed  the  basin  considered,  a major 
observed  storm  might  be  transposed  to  the  basin,  and  the  flood 
resulting  from  this  storm  on  wet  ground  conditions  could  be  assigned 

a reasonable  frequency  and  used  as  an  "anchor  point"  for  extrapolating 
the  frequency  curve. 

d.  While  the  probable  maximum  flood  is  defined  as  the  largest 
flooJ  that  is  reasonably  possible  at  a location,  it  would  not  be 
said  that  a larger  flood  could  absolutely  not  occur.  The  science  of 
meteorology  and  hydrology  has  not  yet  advanced  to  the  stage  where  an 
absolute  maximum  estimate  cart  be  made.  Consequently,  it  is  not 
considered  necessary  that  a frequency  curve  be  limited  to  values 
smaller  than  the  probable  maximum  flood.  However,  it  is  considered 
that  such  a flood  would  have  an  exceedence  interval  considerably  in 
excess  of  1 ,000  years. 
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4-09.  SPECIAL  GEOGRAPHIC  CONSIDERATIONS 

a.  New  England.  Frequency  studies  of  peak  flows  made  in  the 

New  England  area  have  Indicated  that  the  frequency  curves  have  a strong 
positive  skew.  In  studying  this  problem,  two  peculiar  conditions  were 
observed: 

(1)  The  period  of  the  last  40  years,  which  includes  the  period 
of  record  at  most  of  the  stream  gaging  stations,  has  an  abnormally 
large  number  of  major  floods  in  comparison  with  the  100-year  period 
covered  by  the  longer  records  and  in  comparison  with  the  300-year 
period  of  history. 

(2)  Major  floods  occur  from  at  least  two  independent  causes, 
tropical  hurricane  storms  and  extratropical  cyclones.  Hurricane  floods 
are  comparatively  rare,  but  produce  extreme  flows,  and  therefore  cause 
an  upward  curvature  of  the  frequency  curve  of  annual  maximum  flows. 

Some  improvement  in  frequency  estimates  in  this  region  is  attained  by 
segregating  hurricane  and  non-hurricane  floods  (see  reference  21). 
However,  this  apparently  does  not  solve  entirely  the  problem  of  upward 
curvature  of  the  frequency  curves. 

b.  Desert  Regions.  In  many  desert  areas,  streamflow  from  general 
storms  is  ordinarily  moderate  each  year,  but  an  occasional  intense 
thunderstorm  may  center  over  a basin  and  cause  an  exceptionally  large 
flood.  As  in  the  case  of  the  New  England  studies,  floods  from  these 
two  different  causes  tend  to  produce  a sharp  upward  curvature  of  the 
frequency  curve.  While  a satisfactory  solution  to  this  problem  has 
not  been  attained,  a reasonable  approach  appears  to  be  to  estimate  the 
frequency  of  thunderstorms  on  a regional  basis  and  if  desired,  to  combine 
the  estimated  frequency  of  thunderstorm  flows  with  the  frequency  of  flows 
from  general  storms.  In  some  basins,  an  entire  year  may  pass  with  zero 
runoff.  When  analytical  methods  using  flow  logarithms  are  used,  this 
presents  a particular  difficulty,  since  the  logarithm  of  zero  is  minus 
infinity.  It  may  be  best  to  omit  such  years  from  the  record,  compute  a 
tentative  frequency  curve  based  on  the  remaining  years,  and  then  adjust 
the  exceedence  frequency  by  the  ratio  of  the  number  of  years  of  record 

to  the  number  of  years  with  runoff.  However,  factors  resulting  in  zero 
flow  usually  also  affect  small  flows  to  the  extent  that  the  shape  of 
the  frequency  curve  is  distorted  in  the  range  of  lower  flows.  It  is 
possible  to  compute  only  the  upper  half  of  the  frequency  curves,  in 
such  cases,  using  procedures  described  in  paragraph  4-07. 
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c.  Regions  of  Extreme  Variance.  In  several  regions  in  the  United 
States,  runoff  frequency  curves  are  very  steep.  This  occurs  particu- 
larly in  the  Texas  escarpment  area  and  in  the  Southern  California- 
Arizona  area.  In  extrapolating  these  frequency  curves  beyond  the  range 
of  experienced  floods  by  analytical  means,  unreasonably  high  flood 
estimates  may  result.  In  these  cases,  extra  care  must  be  used  in 
extrapolation,  and  hypothetical  computed  floods  might  well  be  used  as 
a general  guide. 
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SECTION  5 - FLOOD  PEAK  FREQUENCY  - ANALYTICAL  ADJUSTMENT 
5-01.  INTRODUCTION 

In  most  cases  of  frequency  studies  of  runoff  or  precipitation, 
there  are  locations  in  the  region  where  records  have  been  obtained 
over  a long  period.  Additional  record  at  a nearby  station  is  useful 
for  extending  the  record  at  the  site  insofar  as  there  is  correlation 
between  recorded  values  at  the  two  locations. 

5-02.  ESTIMATING  INDIVIDUAL  EVENTS 

It  is  possible  by  correlation  or  other  means  to  estimate  from  the 
base  station  values,  the  individual  events  that  were  not  recorded  at 
the  site.  In  doing  this  by  use  of  regression  methods  discussed  in 
section  9,  however,  the  variance  of  the  estimated  values  is  reduced 
by  the  amount  of  non-determination  between  the  two  stations.  In 
frequency  studies,  therefore,  missing  events  should  not  be  freely 
estimated  by  regression  analysis  (and  probably  not  by  graphical 
relations).  In  order  to  insure  that  such  estimates  do  not  unduly 
reduce  the  computed  standard  deviation,  an  estimate  of  an  annual 
maximum  flow  Q]  at  station  1 based  on  a corresponding  annual  maximum 
flow  Qg  at  a base  station  would  be  made  for  this  purpose  by  the 
following  equation: 


X1  " M1  = (X2  " 'VVS2  (6) 

in  which  X represents  the  logarithm  of  a discharge  Q,  Sg , and  S2  are 
the  standard  deviations  of  the  logarithms  of  annual  maximum  flows  for 
concurrent  periods  at  stations  1 and  2 respectively  and  M-j  and  M2  are 
corresponding  mean  values  of  annual  maximum  logarithms. 

5-03.  DEGREE  OF  CORRELATION 

The  direct  means  of  estimating  the  degree  of  correlation  between 
corresponding  flows  at  two  stations  is  to  arrange  successive  pairs  of 
annual  maximum  flows  in  parallel  columns  for  the  concurrent  period  of 
record.  These  flows  should  not  be  arranged  in  the  order  of  magnitude, 
but  should  be  paired  in  their  chronological  sequence.  The  correlation 
coefficient  R is  computed  as  discussed  in  paragraph  9-03  by  use  of  the 
following  equations: 
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1 - R2  = 0 - R2)  pi- 


r2  s ( Sxy ) c s 

2 2 
l)Clyc 


(exy  - ex  EY/Nr 

[ex2  - (ex)2/n][ey2  - (ey)2/n] 


Symbols  are  defined  In  paragraph  9-02  and  In  section  11.  As  discussed 
in  paragraph  9-03e,  a correlation  coefficient  computed  In  the  above 
manner  may  be  unduly  high  or  low,  depending  on  chance  variation  In  the 
data.  When  many  such  correlations  have  been  computed  within  a region, 
It  may  be  possible  to  modify  these  by  judgment  or  regional  correlation 
procedures  in  such  a way  as  to  result  in  a more  reliable  correlation 
estimate.  It  should  be  remembered,  however,  that  this  correlation 
coefficient  usually  has  a minor  influence  on  the  ultimate  frequency 
determination,  and  extensive  studies  designed  to  improve  its  reli- 
ability might  ordinarily  not  be  warranted. 

5-04.  ADJUSTMENT  OF  FREQUENCY  STATISTICS 

In  cases  where  frequency  curves  are  calculated  analytically,  it 
is  neither  necessary  nor  desirable  to  establish  individual  flows  based 
on  nearby  stations,  but  adjustments  can  be  made  in  the  calculated 
statistics  as  follows: 


St  = (S*  - s2)  R‘ 


(approx. ) 


- M1  s <M2  - V r s; 


in  which  the  primes  indicate  the  long-period  values,  and  those  without 
primes  are  based  on  the  same  Short  period  for  both  stations.  Sub- 
scripts indicate  the  station  number.  An  example  of  these  adjustments 
is  shown  on  exhibit  5. 

5-05.  ADVANTAGE  OF  ADJUSTMENT 


The  reliability  of  an  adjusted  value  may  be  expressed  in  terms  of 
the  equivalent  length  of  record  required  to  establish  an  equally 
reliable  unadjusted  value.  The  equivalent  record  derived  from  a 
nearby  station  is  obtained  as  follows: 
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N2  (approx.)  (11) 

Thus,  on  exhibit  5,  use  of  an  additional  17  years  of  record  at  the 
base  station  is  equivalent  to  adding  about  9 years  at  the  site, 
inasmuch  as  the  coefficient  of  determination  is  0.67. 

5-db.  SUMMARY  OF  PROCEDURE 

The  procedure  for  computing  a frequency  curve  using  data  recorded 
at  the  site  and  at  a nearby  long-record  station  is  summarized  in  the 
first  five  steps  in  paragraph  7-10  and  the  four  steps  in  paragraph  4-03, j. 


c. 
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SECTION  6 - FLOOD  VOLUME  FREQUENCY 
6-01 . NATURE  AMD  PURPOSE 

The  comprehensive  flood  volume-duration  frequency  series  consists 
of  a set  of  frequency  curves  as  follows: 

a.  Maximum  rate  of  flow  for  each  water  year. 

b_.  Maximum  1-day  average  flow  for  each  water  year. 

c.  Maximum  3-day  average  flow  for  each  water  year. 

d.  Maximum  10-day  average  flow  for  each  water  year. 

e.  Maximum  30-day  average  flow  for  each  water  year. 

f.  Maximum  90-day  average  flow  for  each  water  year. 

Average  flow  for  each  water  year. 

Runoff  volumes  are  expressed  as  average  flows  in  order  that  peak  flows 
and  volumes  can  be  readily  compared  and  coordinated.  Whenever  it  is 
necessary  to  consider  flows  separately  for  a portion  of  the  water  year 
such  as  the  rain  season  or  snowmelt  season,  the  same  items  (up  to  the 
30-day  or  90-day  values)  are  selected  from  flows  during  that  season 
only.  A comprehensive  flood  volume-duration  series  is  used  primarily 
for  reservoir  design  and  operation  studies,  and  should  be  developed  in 
the  design  of  reservoirs  having  flood  control  as  a major  function. 

When  reservoir  problems  involve  runoff  durations  greater  than  one  year, 
frequency  studies  might  we 1 1 include  multi-annual  runoff  volumes  and 
consideration  of  seasonal  effects,  as  discussed  in  paragraph  6-06. 

6-02.  DATA  FOR  COMPREHENSIVE  SERIES 

Data  to  be  used  for  a comprehensive  flood  volume-duration  frequency 
study  should  be  selected  from  complete  water-year  records  in  accordance 
with  rules  contained  in  paragraph  3-03.  Unless  overriding  reasons  exist, 
durations  specified  in  paragraph  6-01  should  be  used  in  order  to  assure 
consistency  among  various  studies  for  comparison  purposes.  Peak  flows 
should  be  selected  only  for  those  years  when  recorder  gages  existed  or 
when  peak  flows  were  measured  by  other  means.  Where  a minor  portion  of 
a water-year's  record  is  missing,  the  longer-duration  flood  volumes  for 
that  year  can  often  be  estimated  adequately.  Where  upstream  regulation 
or  diversion  exists,  care  should  be  exercised  to  assure  that  each  period 
selected  is  that  when  flows  would  have  been  maximum  under  the  specified 
(usually  natural)  conditions.  The  dates  and  amounts  of  each  selected 
average  flow  should  be  tabulated  in  chronologic  order  in  c.f.s.  or 
thousand  c.f.s..  A typical  tabulation  is  illustrated  on  exhibit  15. 
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6-03.  STATISTICS  FOR  COMPREHENSIVE  SERIES 


The  analytical  method  used  for  flood  volume-duration  frequency 
computations  is  based  on  fitting  the  Pearson  type  III  function  by  use 
of  moments  of  flow  logarithms.  In  practice,  only  the  first  two  moments, 
computed  by  use  of  equations  2 and  3,  are  based  on  station  data.  As 
discussed  in  paragraph  4-03,  the  skew  coefficient  should  not  be  based 
on  a single  record,  but  should  be  derived  from  regional  studies.  The 
following  coefficients,  based  on  studies  summarized  in  references  12  and 
20,  are  considered  to  be  generally  applicable  for  annual  maximum  flood 
volume  frequency  computations: 


Duration 


Skew  Coefficient 


Instantaneous  0 

1 day  -.04 

3 days  -.12 

10  days  -.23 

30  days  -.32 

90  days  -.37 

1 year  -.40 


A sample  computation  of  frequency  statistics  is  given  on  exhibits  15 
and  19.  This  can  be  accomplished  in  steps  as  follows: 


a,.  Tabulate  annual -maximum  average  flows  for  each  duration  in 
chronological  order  as  shown  on  exhibit  15. 


b.  Tabulate  logarithms  in  chronological  order  as  shown  on  exhibit 
17.  If  a long-record  station  nearby  is  available,  tabulate  corresponding 
logarithms  for  that  base  station.  The  work  thus  far  can  be  checked 
approximately  by  ascertaining  that  the  logarithm  for  each  succeeding 
duration  decreases,  but  by  not  more  than  about  0.5.  Much  inconvenience 
can  be  eliminated  by  using  data  for  the  same  years  for  each  duration 
insofar  as  is  feasible.  If  a year's  record  is  incomplete,  the  missing 
portion  can  usually  be  estimated  satisfactorily. 


£.  By  use  of  a statistical  calculator,  the  sums  of  each  column  of 
logarithms,  the  sums  of  their  squares  and  the  sum  of  their  cross-products 
(doubled)  can  be  obtained  for  each  duration  in  a single  cumulative 
operation,  as  described  in  paragraph  9-02b. 
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Those  should  be  entered  on  exhibit  18.  In  cases  of  peak  flows,  advantage 
can  sometimes  be  gained  by  using  1-day  values  at  the  same  station  as  a 
base  instead  of  peak  values  at  the  base  station,  particularly  if  the 
records  of  peak  flows  are  incomplete.  The  extended  1-day  statistics 
would  then  be  used  as  long-term  values. 

d Compute  the  means,  standard  deviations  and  determination 
*c-:.  i cients  for  each  duration.  If  peak-flow  statistics  are  extended 
y .so  of  1-day  flows  at  the  same  station,  possible  inconsistencies 
car.  Lie  avoided  by  using  a determination  coefficient  of  1.00  instead  of 
• *’.ulat1ng  the  coefficient.  Enter  the  long-term  mean  and  standard 
'.:viation  for  each  duration  at  the  base  station  and  compute  the  long- 
term mean  and  standard  deviation  for  each  duration  for  the  station 
concerned.  These  operations  are  shown  on  exhibit  18. 

e.  Plot  the  extended  standard  deviation  against  the  extended  mean, 
and  adopt  a smooth  relationship,  as  shown  on  exhibit  19.  Tabulate  these 
values  on  exhibit  18. 

f_.  When  riany  peak  flows  are  missing  from  the  record,  it  is  best  to 
add  the  average  difference  between  corresponding  peak  and  1-day 
logarithms  to  the  extended  1-day  mean  in  order  to  obtain  the  peak 
mean  logarithm.  The  peak  standard  deviation  can  then  be  obtained  by 
extrapolation  of  the  curve  similar  to  that  shown  in  exhibit  19. 

g.  Select  a skew  coefficient  for  each  duration  from  the  above 
tabulation  or  from  special  regional  studies  as  discussed  in  paragraph 
7-11.  Tabulate  on  exhibit  18. 

6-04.  FREQUENCY  CURVES  FOR  COMPREHENSIVE  SERIES 

a.,  General  procedure.  Frequency  curves  of  flood  volumes  are 
computed  analytically  using  general  principles  and  methods  of  section 
4.  They  should  also  be  shown  graphically  and  compared  with  the  data 
on  which  they  were  based.  This  is  a general  check  on  the  analytic 
work  and  will  ordinarily  reveal  any  inconsistency  in  data  and  method- 
ology. Data  are  plotted  on  a single  sheet  for  comparison  purposes, 
using  procedures  described  in  section  3. 

_b.  Computation  of  basic  curves.  Frequency  curves  are  obtained 
from’  the  frequency  statistics  and  compared  with  observed  frequencies 
for  each  of  the  seven  basic  durations  as  follows: 


- 28  - 


(1)  Tabulate  the  average  flows  for  each  duration  in  order 
of  magnitude  and  obtain  the  plotting  position  for  each  event  fron 
exhibit  37.  This  operation  is  shown  on  exhibit  16. 

(2)  Compute  flows  corresponding  to  related  values  of  P«, 
for  each  frequency  curve  from  the  adjusted  means  and  smoothed 
standard  deviations,  using  equation  5 and  coefficients  obtained  fron 
exhibit  39,  as  shown  on  exhibit  18. 

(3)  Tabulate  values  of  Pfj  from  exhibit  40  for  each  selected 
value  of  Pm  , using  the  average  value  of  M ' for  all  durations. 

(4)  Plot  the  points  obtained  in  step  (1)  and  the  curves  from 
coordinates  obtained  in  steps  (2)  and  (3)  as  shown  on  exhibit  20. 

c.  Interpolation  between  fixed 'durations.  The  runoff  volune  for 
any  specified  frequency  can  be  determined’  for  any  duration  between 
24  hours  and  1 year  by  drawing  a curve  on  logarithmic  paper  as 
illustrated  on  exhibit  21,  relating  volume  to  duration  for  that 
specified  frequency,  using  maximum  24-hour  criteria  derived  in 
reference  15  and  summarized  in  terms  of  average  flov.'s  as  follows: 

log  ^2 4 -hr  = 0,77  109  ^1-day  + 0,23  1og  Vak  (18) 

When  runoff  volumes  for  durations  shorter  than  24  hours  are  very 
important,  special  frequency  studies  should  be  made.  These  could  be 
done  in  the  same  manner  as  for  the  longer  durations,  using  skew 
coefficients  interpolated  in  some  reasonable  manner  between  those 
used  for  peak  and  1-day  flows.  An  approximate  determination  of  short- 
duration  frequencies  is  illustrated  in  reference  15. 

6-05.  DU  RAT  I0IJS  EXCEEDING  1 YEAR 

a.*  Introduction.  In  the  design  of  reservoirs  for  conservation 
purposes  (and  occasional ly  for  flood  control  purposes),  the  volumes  of 
runoff  that  can  be  expected  to  occur  during  the  lifetime  of  the 
structure  within  durations  exceeding  one  year  can  be  of  primary  concern. 
In  general,  the  design  of  such  reservoirs  has  been  based  on  the  lowest 
(or  highest)  volume  of  runoff  observed  for  the  critical  duration  during 
the  period  of  record,  which  generally  encompasses  40  to  100  years. 
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Frequently,  however,  It  is  felt  that  the  observed  minimum  or  maximum 
is  much  more  extreme  than  should  normally  be  anticipated  in  a similar 
future  period.  In  order  to  determine  expected  volumes  with  a greater 
degree  of  reliability,  frequency  studies  of  long-duration  volumes  can 
be  made.  However,  direct  frequency  analysis  is  ordinarily  not 
practicable,  because  a relatively  small  number  of  independent  long- 
duration  volumes  is  contained  in  a single  record.  (Or.ly  20  independ- 
ent 5-year  volumes  are  contained  in  a 100-year  record,  for  example.) 
In  order  to  overcome  this  limitation,  a study  has  been  made  under  the 
Civil  Works  Investigation  program  of  the  Corps  of  Engineers  to  relate 
long-duration  volumes  to  annual  volumes.  Using  criteria  developed  in 
this  study,  long-duration  volumes  can  be  derived  from  a frequency 
curve  of  annual  volumes  and  the  correlation  coefficient  between 
successive  annual  flows.  Where  this  correlation  coefficient  is 
considered  to  be  zero,  the  criteria  should  be  fairly  dependable,  but 
where  there  appears  to  be  substantial  correlation  between  successive 
annual  runoff  values,  the  reliability  decreases,  principally  because 
of  the  uncertainty  as  to  the  true  correlation  coefficient.  A 
description  of  the  study  and  the  derived  criteria  is  contained  in 
reference  18. 

b.  Annual  runoff.  Criteria  for  determining  multi -annual  runoff 
are  based  on  the  logarithmic  mean  and  standard  deviation  of  annual 
runoff.  These  should  be  computed  as  described  in  paragraph  6-04. 

c.  Persistence  effects.  In  many  river  basins,  surface  or  sub- 
surface storage  effects  cause  the  flows  in  one  year  to  reflect 
conditions  in  the  preceding  year  to  some  extent.  This  will  result 
in  a positive  correlation  between  successive  years1  runoff.  A 
measure  of  this  persistence  effect  is  the  correlation  coefficient, 
or  preferably  its  square,  the  determination  coefficient,  between 
successive  years'  runoff  logarithms.  This  is  determined  by  pairing 
each  year's  runoff  logarithm  (except  the  first)  with  that  of  the 
preceding  year,  and  computing  the  determination  coefficient  by  use 
of  equations  7 and  8.  Because  of  the  important  effect  of  the  deter- 
mination coefficient  on  long-duration  volume  estimates,  its  degree 
of  unreliability  as  discussed  in  paragraph  10-04  should  be  given 
special  consideration. 

d_.  liulti-annual  runoff.  A frequency  curve  of  total  runoff 
volume  expected'  to  occur  during  a period  consisting  of  an  integral 
number  of  water  years  can  be  obtained  directly  from  the  mean  and 
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standard  deviation  of  annual  (water-year)  runoff  logarithms  and  the 
determination  coefficients  between  successive  annual  runoff  logarithms, 
by  use  of  chart  14  of  reference  10.  In  using  such  a frequency  curve, 
it  should  be  remembered  that  there  are  (N  - T + 1)  T-year  periods  in 
an  1 t-year  record.  For  example,  only  46  different  5-year  volumes  can 
occur  in  a 50-year  period. 

e.  Seasonal  effects.  It  must  be  recognized  that  multi-annual 
runoTf  based  on  water-year  volumes  is  not  representative  of  critical 
conditions.  A drier  or  wetter  period  of  the  same  duration  can  usually 
be  found  by  starting  the  period  some  days  or  months  earlier  or  later. 
Also,  an  integral  number  of  years  does  not  represent  a critical  duration, 
because  adding  one  more  dry  season  (or  wet  season)  will  ordinarily  worsen 
the  condition.  For  these  reasons,  criteria  for  determining  maximum  or 
minimum  runoff  volumes  for  any  duration  between  1 and  20  years  and  for 
any  frequency  were  derived,  as  described  in  reference  1G,  and  are 
summarized  on  chart  15  of  that  report. 

£.  Illustrative  example.  A sample  computation  of  runoff  volume 
frequencies  for  durations  longer  than  1 year  is  illustrated  on  chart  16 
of  reference  18. 

6-06.  APPLICATIONS  OF  FLOOD  VOLUME-DURATION  FREQUENCIES 

a_.  Volume-duration  curves.  The  use  of  flood  volume-duration 
frequencies  in  solving  reservoir  planning,  design,  and  operation 
problems  usually  involves  the  construction  of  volume-duration  curves 
for  specified  frequencies.  These  are  drawn  first  on  logarithmic 
paper  for  interpolation  purposes,  as  discussed  in  paragraph  6-04c  and 
illustrated  on  exhibit  21,  and  are  then  replotted  on  arithmetic  grid 
as  shown  on  the  same  exhibit.  A volume-duration  curve  for  durations 
longer  than  1 year  is  illustrated  on  chart  17  of  reference  18.  A 
straight  line  on  this  grid  represents  a constant  rate  of  flow  (so  many 
acre-feet  per  day).  The  straight  lines  on  exhibit  21  represent  a 
uniform  flow  of  2,000  c.f.s.,  and  placement  on  the  100-ycar  volume- 
duration  curve  demonstrates  that  a reservoir  capacity  of  36,000  acre- 
feet  is  required  to  control  the  indicated  runoff  volumes  to  a project 
release  of  2,000  c.f.s..  The  curve  also  indicates  that  durations  of 
4 to  7 days  are  critical  for  this  project  release  and  flood  control 
space. 
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b_.  Simple  reservoir  problems.  In  the  case  of  a single  flood 
control  reservoir  located  Immediately  upstream  of  the  only  important 
damage  center  concerned,  the  volume  frequency  problems  are  relatively 
simple.  A series  of  volume-duration  curves  similar  to  that  shown  on 
exhibit  21  corresponding  to  selected  frequencies  should  first  be  drawn. 
The  project  release  rate  should  be  determined,  giving  due  consideration 
to  possible  channel  deterioration,  encroachment  into  the  flood  plain, 
and  operational  contingencies.  Lines  representing  this  flow  rate  are 
then  drawn  tangent  to  each  volume-duration  curve,  and  the  intercept  in 
each  case  determines  the  reservoir  space  used  to  control  the  flood  of 
that  selected  frequency.  The  point  of  tangency  represents  the  critical 
duration  of  runoff.  This  procedure  can  be  used  not  only  as  an 
approximate  aid  in  selecting  a reservoir  capacity,  but  as  an  aid  in 
drawing  filling-frequency  curves. 

c . Complex  reservoir  problems.  Where  reservoir  operation 
schedules  are  variable  or  where  many  reservoirs  are  operated  jointly, 
it  may  be  necessary  to  route  historical  flows  month  by  month  or  day  by 
day  in  order  to  demonstrate  the  adequacy  of  a design  or  operation 
procedure,  however,  the  techniques  described  in  the  preceding  para- 
graph may  bo  applicable  approximately,  and  may  shod  considerable 
light  on  the  problem.  In  applying  such  techniques,  the  following 
quidos  should  be  used; 

(1)  Volume-duration  curves  are  needed  for  unregulated  flows 
at  each  important  damage  center. 

(2)  The  straight  line  corresponding  to  the  average  nondamaging 
flow,  allowing  for  operational  contingencies,  when  drawn  tangent  to  the 
volume-duration  curve  corresponding  to  the  selected  design  frequency 
will  indicate  the  storage  required  in  the  system  if  it  is  located  and 
operated  so  as  to  be  fully  effective. 

(3)  The  same  straight  line  will  indicate  the  range  of  critical 
uurotions  for  design  and  operation  studies.  A system  f reservoirs 
should  be  "tuned"  to  this  range  of  durations,  insofar  ;s  is  feasible, 
because  a reservoir  that  fills  and  empties  in  5 days  may  be  of  no  value 
if  the  critical  duration  at  a downstream  damage  center  is  15  days. 
Likewise,  a reservoir  that  is  only  half  full  in  15  days  would  not  have 
provided  its  best  control  at  the  damage  center. 
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d.  Representative  hydroqraphs.  In  solving  complex  reservoir 
problems' , re p re sentat.i ve  I lydrog r a p h s at  all  locations  can  be  patterned 
after  one  or  more  past  floods.  The  ordinates  of  these  hydrographs  can 
be  adjusted  so  that  their  volumes  for  the  critical  durations  will  equal 
corresponding  magnitudes  at  each  location  for  the  selected  frequency. 

A design  or  operation  scheme  based  on  regulation  of  such  a set  of 
hydrographs  would  be  reasonably  well  balanced. 
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SECTION  7 - REGIONAL  FREQUENCY  ANALYSIS 


7-01 . GENERAL 

Runoff  frequency  estimates  for  any  specific  location  based  on  data 
recorded  at  that  location  and,  where  practicable,  adjusted  by  use  of 
longer-record  stations  nearby  as  described  in  section  5,  can  ordinarily 
be  improved  by  a study  of  frequency  characteristics  throughout  the 
hydrologic  region.  Such  a regional  study  will  also  help  to  assure 
consistency  of  estimates  for  different  locations  and  will  provide  means 
for  estimating  frequencies  at  locations  where  data  are  not  available. 
Where  time  is  limited  and  approximate  estimates  are  satisfactory,  some 
very  simple  schemes  are  helpful.  Where  time  permits,  more  elaborate 
schemes  are  often  justified. 

7-02.  USE  OF  FREQUENCY  STATISTICS 

A regional  frequency  correlation  study  is  based  on  the  two  principal 
frequency  statistics  - the  mean  and  standard  deviation  of  annual  maximum 
flow  logarithms.  Prior  to  relating  these  frequency  statistics  to 
drainage-basin  characteristics,  it  is  essential  that  the  best  possible 
estimate  of  each  frequency  statistic  be  made.  This  is  done  by  adjusting 
short-record  values  by  the  use  of  longer  records  at  nearby  locations- 
When  many  stations  are  involved,  it  is  best  first  to  select  long-record 
base  stations  for  each  portion  of  the  region.  It  might  be  desirable  to 
adjust  the  base  station  statistics  by  use  of  the  one  or  two  longest- 
record  stations  in  the  region,  and  then  adjust  the  short-record  station 
values  by  use  of  the  nearest  or  most  appropriate  base  station.  Methods 
of  adjusting  statistics  are  discussed  in  section  5. 

7-03.  SIMPLE  SCHEMES 

For  preliminary  studies  where  a high  degree  of  accuracy  is  not 
required,  regional  frequency  analyses  might  consist  simply  of  plotting 
the  standard  deviation  against  drainage  area  size  or,  for  various 
locations  on  the  same  river,  against  river  mile  distance,  if  preferred. 
Similarly,  the  mean  logarithm  representing  general  magnitude  can  be 
plotted  against  drainage  area  size.  Another  simple  scheme  is  to  plot 
the  standard  deviation  or  mean  logarithm  of  discharge  per  square  mile 
on  a map,  and  drawing  lines  of  equal  standard  deviation  or  mean 
logarithm.  Such  an  analysis  can  be  used  to  modify  the  estimates 
slightly  to  improve  consistency  and  to  select  statistics  for  ungaged 
areas  in  the  region.  This  type  of  analysis  is  illustrated  in 
reference  20. 
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7-04.  DRAINAGE-BASIN  CHARACTERISTICS 

A regional  analysis  involves  the  determination  or  the  main  factors 
responsible  for  differences  in  precipitation  or  runoff  regimes  between 
different  locations.  This  is  done  by  correlating  Important  factors 
with  the  long-record  mean  and  with  the  long-record  standard  deviation 
of  the  frequency  curve  for  each  station.  (The  long-record  values  are 
those  based  on  extension  of  the  records  as  discussed  in  section  5.) 
Statistics  based  on  rainfall  measurements  may  be  correlated  in  moun- 
tainous terrain  with  the  following  factors: 

a.  Elevation  of  station 

b.  General  slope  of  surrounding  terrain 
c^.  Orientation  of  that  slope 

d.  Elevation  of  windward  barrier 

e.  Exposure  of  gage 

T.  Distance  to  leeward  controlling  ridge 

Statistics  based  on  runoff  measurements  may  be  correlated  with  the 
following  factors: 

a.  Drainage  area  (contributing) 

F.  Slope  of  drainage  area  or  of  main  channel 

c.  Surface  storage  (lakes  and  swamps) 

7.  Mean  annual  rainfall 

e.  Number  of  rainy  days  per  year 

f.  Infiltration  characteristics 
£.  Stream  length 

7-05.  CORRELATION  METHODS 

Correlation  methods  and  their  application  are  discussed  in 
section  9. 

7-06.  LINEAR  RELATIONSHIPS 

In  order  to  obtain  satisfactory  results  using  multiple  linear 
correlation  techniques,  all  variables  must  be  expressed  so  that 
the  relation  between  the  dependent  and  any  independent  variable 
can  be  expected  to  be  linear,  and  so  that  the  interaction  between 
two  independent  variables  is  reasonable.  An  illustration  of  the 
first  condition  is  the  relation  between  rainfall  and  runoff.  If 
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the  runoff  coefficient  is  sensibly  constant,  as  in  the  case  of  urban 
or  airport  drainage,  then  runoff  can  be  expected  to  bear  a linear 
relation  to  rainfall.  However,  in  many  cases  initial  losses  and 
infiltration  losses  cause  a marked  curvature  in  the  relationship. 
Ordinarily,  it  will  be  found  that  the  logarithm  of  runoff  is  very 
nearly  a linear  function  of  rainfall,  regardless  of  loss  rates,  and 
In  such  cases,  linear  correlation  of  logarithms  would  be  most  suitable. 
An  illustration  of  the  second  condition  is  the  relation  between  rain- 
fall, 0,  drainage  area,  A,  and  runoff,  Q.  If  the  relation  used  for 
correlation  is: 

Q s aD  + bA  + c (13) 

then  it  can  be  seen  that  one  inch  change  in  precipitation  would  add 
the  same  amount  of  flow,  regardless  of  the  size  of  drainage  area. 

This  is  not  reasonable,  but  again  a transformation  to  logarithms 
would  yield  a reasonable  relation: 

Log  Q = d log  D + e log  A + log  f (14) 

or  transformed: 


Q = fDdAe  (15) 

Thus,  if  logarithms  of  certain  variables  are  used,  doubling  one 
independent  quantity  will  multiply  the  dependent  variable  by  a 
fixed  ratio,  regardless  of  what  fixed  value  the  other  independent 
variables  have.  This  particular  relationship  is  reasonable  and 
can  be  easily  visualized  after  a little  study.  There  is  no  simple 
rule  for  deciding  when  to  use  the  logarithmic  transformation.  It 
is  only  feasible,  however,  when  the  variable  has  a fixed  lower 
limit  of  zero. 

7-07.  EXAMPLE  OF  REGIONAL  CORRELATION 

An  illustrative  example  of  a regional  correlation  analysis  of 
standard  deviation  with  drainage  area  and  number  of  rainy  days 
per  year  is  given  on  exhibit  22.  Since  many  important  variables 
are  neglected,  the  analysis  is  not  of  the  scope  necessary  for  a 
complete  study,  but  is  useful  for  illustrating  various  techniques 
and  problems  involved  in  such  a study.  In  the  example,  Xj  is  one 
plus  the  logarithm  of  the  adjusted  standard  deviation  (one  is  added 
to  eliminate  negative  values),  Xj  is  the  logarithm  of  the  drainage 
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area  size,  and  X3  Is  the  logarithm  of  the  average  number  of  rainy  days 
per  year  for  the  drainage  area.  The  regression  equation  Is  derived 
as  shown,  and  the  calculted  coefficient  of  determination  Is  0.31, 
which  means  that  31  percent  of  the  variance  of  X,  Is  explained  by 
the  regression  equation. 

7-08.  SELECTION  OF  USEFUL  VARIABLES 

In  the  regression  equation  derived  on  exhibit  22,  the  coefficient 
of  X2  Is  very  small,  which  Indicates  that  this  factor  has  very  little 
effect.  To  determine  the  usefulness  of  this  factor,  It  is  necessary 
to  make  an  additional  analysis  using  all  variables  except  this  one. 

In  this  case,  the  problem  would  resolve  into  a simple  correlation 
analysis  using  X,  of  exhibit  22  as  X,  in  equations  19  to  21. 

Then:  J 1 

b2  = -1.1749/2.3484  = -0.50  (16) 

and 

a = 0.3578  - (-.50)1.925  = 1.32  (17) 

Hence 

Xi  * 1.32  - O.SOXj  (18) 

— ? 

A solution  for  R (equations  31  and  33)  would  yield  0.32.  Thus,  a 
better  correlation  is  obtained  neglecting  drainage  area  as  a factor. 

If  additional  factors  were  considered  in  the  analysis,  the  effect  of 
drainage  area  should  be  reconsidered,  as  it  is  possible  that  its 
effect  is  obscured  in  the  example  by  neglecting  some  other  important 
variable.  The  final  test  of  importance  of  a particular  factor  is  a 
comparison  of  the  correlation  coefficient  using  all  factors  and  then 
omitting  only  the  factor  whose  influence  is  being  tested.  Even  in 
the  case  of  a slight  increase  in  correlation  obtained  by  adding  a 
variable,  consideration  of  the  increased  unreliability  of  R,  as 
discussed  in  paragraph  10-04,  might  indicate  that  such  factor  should 
be  eliminated  in  cases  of  small  samples. 
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7-09.  USE  OF  MAPS 

Many  hydrologic  factors  cannot  be  expressed  numerically.  Examples 
are  soil  characteristics,  vegetal  cover,  and  geology,  For  this  reason, 
numerical  regional  analysis  will  explain  only  a portion  of  the  regional 
variation  of  runoff  frequencies.  The  remaining  unexplained  variance 
is  contained  in  the  regression  constant,  which  can  be  considered  to 
vary  from  station  to  station.  These  regression  constants  (residuals) 
can  be  computed  by  inserting  the  drainage  basin  characteristics  and 
frequency  statistic  for  each  station  in  the  regression  equation  and 
solving  for  the  regression  constant.  These  constants  can  be  plotted 
on  a regional  map,  and  lines  of  equal  values  drawn  (perhaps  using  soils 
or  vegetation  maps  as  a guide).  Use  of  such  a map  for  selecting  a 
regression  constant  should  be  much  better  than  using  the  single 
constant  for  the  entire  region  derived  from  equation  28.  In  smoothing 
lines  on  such  a map,  consideration  should  be  given  to  the  reliability 
of  computed  statistics.  Exhibit  23  shows  standard  errors  of  estimating 
means  and  standard  deviations.  As  an  example,  If  a computed  standard 
deviation  based  on  30  years  of  record  is  0.300,  there  is  about  one 
chance  in  20  that  the  mean  is  in  error  by  more  than  0.110  (twice  the 
standard  error)  or  that  the  standard  deviation  is  in  error  by  a factor 
of  1.3  (antilog  of  0.114). 

7-10.  SUMMARY  UF  PROCEDURE 

A regional  analysis  of  precipitation  or  floodflows  is  accomplished 
in  the  following  steps: 

a^.  Select  long-record  base  stations  within  the  region  as  required 
for  extension  of  records  at  each  of  the  short-record  stations. 

b_.  Tabulate  maximum  events  of  each  station,  corresponding 
logarithms,  and  logarithms  of  base-station  values  for  the  corres- 
ponding years.  Logarithms  should  be  rounded  to  2 decimal  places. 

c.  Calculate  M and  S (equations  2 and  3)  for  each  base  station. 

d.  Calculate  M and  S for  each  other  station  and  for  the 
corresponding  values  of  the  base  station,  and  calculate  the  cor- 
relation coefficient  (equations  31  and  33).  Summation  of  logarithms 
and  their  squares  for  both  stations  and  their  cross-products  can 

1 
3 
3 
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be  obtained  In  a single  cumulative  operation  of  an  automatic  calculator, 
as  discussed  in  paragraph  9-02B. 

e.  Adjust  all  values  of  M and  S by  use  of  the  base  station, 
(equations  9 and  10).  (If  any  base  station  is  first  adjusted  by  use 
of  a longer-record  base  station,  the  longer-record  statistics  should 
be  used  for  all  subsequent  adjustments.) 

f.  Select  meterological  and  drainage  basin  parameters  that  are 
expected  to  correlate  linearly  with  M and  log  S,  and  tabulate 
estimated  values  of  these  for  each  area.  (The  physical  significance 
of  log  S is  not  important,  as  the  transformation  simply  eliminates  a 
lower  limit  of  zero  from  the  regression  variables.) 

£.  Calculate  the  regression  equations  relating  fi  and  log  S in 
turn  to  these  statistics,  using  procedures  explained  in  section  9, 
and  compute  the  corresponding  determination  coefficients. 

h,.  Eliminate  variables  in  turn  that  contribute  the  least  to  the 
determination  coefficient,  recomputing  the  determination  coefficient 
each  time,  and  select  the  regression  equation  having  the  highest 
determination  coefficient,  or  one  with  fewer  variables  if  the  deter- 
mination coefficient  is  about  as  high. 

i_.  Compute  the  regression  constants  (residuals)  for  each  station, 
plot  on  a suitable  map,  and  draw  isopleths  of  tru  regression  constant 
for  the  regression  equations  of  M and  S (two  maps),  considering  that 
the  regression  constant  for  a station  represents  a basin-mean  value. 

j.  A frequency  curve  can  be  computed  from  constants  obtained  for 
any  basin  on  the  map,  using  the  computed  regression  equations  to 
obtain  M and  S,  and  using  procedures  discussed  in  paragraph  4-03  for 
computing  a frequency  curve  therefrom. 

7-11.  REGIONAL  SKEW  DETERMINATIONS 

Skew  coefficients  for  use  in  hydrologic  studies  should  be  based 
on  regional  studies,  since  values  based  on  individual  records  in 
the  order  of  100  years  or  less  of  length  are  highly  unreliable. 

This  can  be  done  by  computing  skew  coefficients  for  available 
records  and  using  the  average,  weighted  in  accordance  with  record 
length.  Unless  an  electronic  computer  is  available,  such  procedure 
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Is  ‘laborious.  An  alternative  procedure  suitable  for  desk  calculation 
is  to  compute  an  average  skew  coefficient  as  illustrated  on  chart  4 
of  reference  14. 


7-12.  SAMPLE  REGIONAL  CRITERIA 


The  regional  frequency  analysis  described  in  reference  20  is 
considered  to  be  a moderately  elaborate  type  of  analysis.  The  criteria 
derived  are  suitable  for  selecting  frequency  statistics  for  deriving 
frequency  curves  of  runoff  peaks  and  volumes  for  durations  up  to  30 
days.  Standard  deviations  for  frequency  curves  for  the  various 
durations  of  runoff  are  obtained  directly  from  maps  constructed  for 
each  duration.  The  mean  logarithms  of  runoff  for  each  duration, 
however,  have  been  related  to  drainage  area  size,  normal  annual 
precipitation,  elevation,  and  by  means  by  maps,  to  geographical 
location.  To  illustrate  the  relative  simplicity  of  this  scheme, 
criteria  for  determining  frequency  curves  of  peak  flows  are  Included 
herein  as  exhibits  24  to  27.  Data  and  computations  required  for 
synthetic  frequency  curves  for  the  drainage  basin  used  in  illustrat- 
ing graphical  and  analytical  methods  on  exhibits  2 to  4 are  as  follows: 

Drainage  basin  characteristics  (location  shown  on  exhibit  24): 


Drainage  area 

Normal  annual  precipitation 
Average  elevation 
Average  latitude 


Frequency  constants: 

Standard  deviation 
(from  exhibit  25) 
Cp  (from  exhibit  24) 

K (from  exhibit  26) 
Computation  of  Q 


134  sq.  mi . 

47  in. 

2900  ft.  m.s.l . 
40°  - 03‘ 


.31 

42 

.71 


(equation  11,  reference  20) 


Q = .001  C A’85  P2K 
P P 

= .001  (42)  (64)  (47)4  (.71) 
s 4220  c.f.s. 


Flow  of  specified  frequency  (say  once  per  100  years 

q.Oi  = (4»220)  = 24,000  c.f.s.  (from  exhibit  27) 
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SECTION  8 - FREQUENCIES  OF  OTHER  HYDROLOGIC  FACTORS 
8-01.  INTRODUCTION 

The  frequency  methods  described  in  sections  3 and  4 have  appli- 
cation to  runoff,  and  can  also  be  used  in  estimating  frequencies  of 
various  other  hydrologic  factors.  Some  of  the  more  common  applications 
are  described  in  the  following  paragraphs. 

8-02.  RAINFALL  FREQUENCIES 

Procedures  for  the  computation  of  frequency  curves  of  station 
precipitation,  both  graphically  and  analytically,  are  generally 
identical  to  those  for  streamflow  analysis.  In  precipitation  studies, 
however,  instantaneous  peak  intensities  are  ordinarily  not  analyzed, 
since  they  are  virtually  impossible  to  measure  and  of  little  application. 
Precipitation  amounts  for  specified  durations  are  commonly  analyzed, 
mostly  for  durations  of  less  than  three  or  four  days.  The  few  studies 
made  thus  far  have  indicated  that  the  logarithmic  normal  function 
(with  zero  skew  coefficient)  is  fitted  fairly  well  with  annual  maximum 
station  precipitation  data,  regardless  of  the  duration  used.  Station 
precipitation  alone  is  not  adequate  for  most  hydrologic  studies,  and 
some  means  of  evaluating  the  frequency  of  simultaneous  or  near 
simultaneous  precipitation  over  the  area  is  necessary. 

8-03.  LOW  FLOW  FREQUENCIES 

a.  The  design  of  hydroelectric  powerplants  and  the  design  of 
reservoirs  for  supplementing  low  river  flows  for  water  quality  and 
other  purposes  requires  the  evaluation  of  the  frequencies  of  low 
flows  for  various  durations.  The  method  of  frequency  analysis 
previously  discussed  is  usually  applicable,  except  that  minimum 
instead  of  maximum  runoff  for  each  period  is  selected  from  the  basic 
data.  In  studying  low  flows,  it  will  be  found  that  the  effects  of 
basin  development  are  relatively  great.  For  example,  a relatively 
moderate  diversion  can  be  neglected  when  studying  floodflows,  but 
might  greatly  modify  or  even  eliminate  low  flows.  Accordingly,  one 
of  the  most  important  aspects  of  low  flows  concerns  the  evaluation 
of  past  and  future  effects  of  basin  developments. 

b.  Civil  Works  Investigations  Project  No.  154  has  been  estab- 
lished by  the  Corps  of  Engineers  to  study  low  flows  and  their 
frequency,  and  has  met  with  considerable  success  in  applying  analy- 
tical procedures  described  herein.  These  studies  have  not  progressed 
sufficiently  to  provide  firm  criteria,  however. 
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c.  In  regions  of  water  scarcity  and  where  a high  degree  of 
development  has  been  attained,  basin  development  projects  that  entail 
carryover  of  water  for  several  years  are  often  planned.  A study  of 
frequencies  of  flows  having  durations  exceeding  one  year  is  contained 
In  reference  18  and  may  be  of  considerable  use  in  connection  with 
such  projects. 

8-04.  HURRICANE  FREQUENCIES 

In  studies  of  hurricane  wind  velocities  by  tne  Corps  of  Engineers 
and  U.S.  Weather  Bureau,  the  central  pressure  index  (estimated  minimum 
sea  level  pressure  for  individual  hurricane)  has  been  used  in  conjunction 
with  pressure  vs.  wind  relationships  to  determine  wind  frequencies.  The 
index  frequency  for  each  of  three  large  geographic  zones  was  determined 
as  illustrated  on  exhibit  28,  and  subsequently  divided  into  frequencies 
for  various  subdivisions  of  each  zone.  The  minimum  hurricane  pressure 
apparently  plots  close  to  a straight  line  on  arithmetic  probability 
paper. 

8-05.  SEDIMENT  FREQUENCIES 

Another  illustration  of  use  of  frequency  techniques  discussed  herein 
is  shown  on  exhibit  29,  where  the  frequency  of  annual  sediment  load  of 
the  Colorado  River  is  shown  to  approximate  a linear  relationship  on 
log  probability  paper. 

8-06.  COINCIDENT  FREQUENCIES 

a_.  In  many  cases  of  hydrologic  design,  it  is  necessary  to  consider 
only "those  events  which  occur  coincidentally  with  other  events.  For 
example,  a pumping  station  is  usually  required  to  pump  water  only  when 
interior  drainage  occurs  at  a time  that  the  main  river  stage  is  above 
the  gravity  outlet.  In  constructing  a frequency  curve  of  interior 
drainage  flows  that  occur  only  at  such  times,  data  selected  for  direct 
use  should  be  limited  to  that  recorded  during  high  river  stages.  In 
some  cases,  such  data  might  be  adequate,  but  it  is  usually  possible  in 
cases  where  the  two  types  of  events  do  not  correlate  to  make  indirect 
use  of  noncoincident  data  in  order  to  establish  a more  reliable 
frequency  curve  of  coincident  events.  The  general  procedure  used  is 
as  follows: 

(1)  Select  the  more  stable  of  the  two  variables  whose 
coincidental  frequency  is  to  be  determined.  This  will  be  designated 
as  variable  B,  and  the  other  as  variable  A. 
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(2)  Determine  the  time  limits  of  the  seasons  during  which 
high  stages  of  variable  B are  about  uniformly  likely.  If  all  important 
stages  of  variable  B can  be  limited  to  one  season,  the  analysis  will 

be  simplified.  Otherwise  the  following  steps  must  be  duplicated  for 
each  season. 

(3)  Compute  a frequency  curve  of  stages  (or  flows,  rainfall 
amounts,  etc.)  for  variable  A,  using  data  obtained  only  during  the 
selected  seat'v.i. 

(4)  Construct  a duration  curve  of  stages  (or  flows,  etc.) 
for  variable  B,  using  data  obtained  only  during  the  selected  season. 

(5)  The  exceedence  frequency  of  any  selected  magnitude  of 
variable  A that  is  coincidental  with  any  specified  range  of  stage  for 
variable  B is  equal  to  the  product  of  the  exceedence  frequency 
indicated  by  the  curve  derived  in  (3)  and  the  proportion  of  time 
flows  at  B are  within  the  selected  range  of  stage,  as  indicated  by 
the  curve  derived  in  (4). 

Exhibit  30  illustrates  a computation  of  reservoir  stage  or 
storage  frequency  curve  from  consideration  of  coincidental  frequencies. 
In  this  hypothetical  case  of  a 300,000  acre-foot  reservoir,  the  top 
100,000  acre-feet  is  reserved  for  flood  control,  the  next  50,000  for 
seasonal  irrigation  requirements,  and  the  remainder  for  carryover 
irrigation  storage  and  power  head.  Because  flood  control  space  would 
rarely  be  used,  routings  of  recorded  floods  do  not  adequately  define 
the  frequency  of  storage  in  the  flood  control  space.  The  example  shows 
a flood  frequency  curve  in  terms  of  the  project  design  (50-year)  flood, 
and  a storage  duration  curve  determined  from  monthly  routings  of 
recorded  runoff.  This  latter  curve  represents  only  those  months  when 
major  storms  are  likely  to  occur,  and  does  not  reflect  reservoir  rises 
that  would  occur  during  floods.  Therefore,  it  represents  conditions 
that  can  exist  at  the  beginning  of  a flood.  The  duration  curve  was 
divided  into  four  ranges,  and  average  storages  determined  for  each 
range.  Routings  (not  shown)  of  various  percentages  of  the  project 
design  flood  with  these  four  initial  storages  were  made,  and  four 
stage-frequency  curves  drawn  as  shown.  A composite  frequency  curve 
was  then  drawn  as  illustrated  in  the  inset  table. 
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SECTION  9 - CORRELATION  ANALYSIS 
9-01.  NATURE  AND  APPLICATION 

a. .  Correlation  is  the  process  of  determining  the  manner  in  which 
the  changes  in  one  or  more  independent  variables  affect  another 
(dependent)  variable.  The  dependent  variable  is  the  value  sought  and 
is  to  be  related  to  various  independent  variables  which  will  be  known 
in  advance,  and  which  will  be  physically  related  to  the  dependent 
variable.  For  example,  the  volume  of  spring  runoff  on  a river 
(dependent  variable)  might  be  correlated  with  the  depth  of  snow  cover 
in  the  area  (independent  variable).  Recorded  values  of  such  variables 
over  a period  of  years  might  be  plotted  as  a graph  and  the  apparent 
relation  sketched  in  by  eye.  However,  correlation  methods  will 
generally  permit  a more  dependable  determination  of  the  relation  and 
have  the  additional  advantage  of  providing  means  for  evaluating  the 
dependability  of  the  relation  or  of  estimates  based  on  the  relation. 

b. .  The  function  relating  the  variables  is  termed  the  "regression 
equation,"  and  the  proportion  of  the  "variance"  of  the  dependent 
variable  that  is  explained  by  the  regression  equation  is  termed  the 
"coefficient  of  determination,"  which  is  the  square  of  the  "correla- 
tion coefficient."  Regression  equations  can  be  linear  or  curvi- 
linear, but  linear  regression  suffices  for  most  applications,  and 
curvilinear  regression  is  therefore  not  discussed  herein. 

9-02.  CALCULATION  OF  REGRESSION  EQUATIONS 

a.  In  a simple  correlation  (one  in  which  there  is  only  one 
independent  variable),  the  linear  regression  equation  is  written 

X1  s a + b2X2  (19) 

in  which  Xi  is  the  dependent  variable,  X2  is  the  independent  variable, 
a is  the  regression  constant,  and  b2  is  the  regression  coefficient. 

The  coefficient  b2  is  evaluated  from  the  tabulated  data  by  use  of  the 
equation 

b2  - t(x1x2)/t(x2)2  (20) 

in  which  x-)  is  the  deviation  of  a single  value  X]  from  the  mean  M]  of 
its  series,  and  X2  is  similarly  defined.  The  regres  ion  constant  is 
obtained  from  the  tabulated  data  by  use  of  the  equation 
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b.  All  summations  required  for  a simple  linear  correlation  can 
be  obtained  In  a single  cumulative  operation  of  a single-keyboard  10- 
bank  automatic  calculator,  using  equations  29  and  30,  as  follows: 

(1)  Round  all  values  of  X}  to  the  same  decimal  place  so  that 
the  median  value  has  two  significant  figures.  Logarithms  should  be 
rounded  to  two  decimal  places.  Repeat  for  values  of  X2.  If  there  are 
negative  values  of  any  variable,  add  a constant  to  all  values  of  that 
variable  and  subsequently  subtract  that  constant  from  the  mean. 

(2)  Enter  the  first  value  of  X}  on  the  left  of  the  keyboard 
using  banks  2,  3 and  if  necessary,  bank  1.  Enter  the  corresponding 
value  of  X2  on  the  right  of  the  keyboard  using  banks  9,  10  and  if 
nceessary,  bank  8.  Square  the  quantity  on  the  keyboard. 

(3)  Lock  both  the  multiplier  and  product  dials  and  repeat 
the  process  with  each  pair  of  values.  When  all  values  have  been 
squared,  the  sum  of  the  X]  and  X£  values  will  appear  on  the  multiplier 
dial,  the  sum  of  their  squares  will  appear  at  each  end  of  the  product 
dial,  and  twice  the  sum  of  their  cross-products  will  appear  in  the 
middle  seven  digits  of  the  product  dial.  Care  should  be  exercised 
that  no  cumulative  product  exceeds  the  machine  capacity  (7  digits, 
usually). 

c_.  In  a multiple  correlation  (one  in  which  there  is  more  than 
one  independent  variable)  the  linear  regression  equation  is  written 

Xi  = a + b2X2  + b3X3  . . . + bnXn  (22) 

In  the  case  of  two  independent  variables,  the  b coefficients  are 
evaluated  from  the  tabulated  data  by  solution  of  the  following 
simultaneous  equations: 


2 

t(x2)  b2  + t(x2x3)b3  = Z(x1x2) 

(23) 

E(x2x3)b2  + £(x3)2b3  s E(XjX3) 

(24) 
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f In  the  case  of  3 independent  variables,  the  b coefficients  can  be 

| evaluated  from  the  tabulated  data  by  solution  of  the  following 

I simultaneous  equations: 

^(Xg)  62  £(x2x3)b3  + £ (x2x4 ) ^4  3 (25) 

■ £ ( X2X3 ) ^2  * ^(x3^  ^3  ^ e(x3x4^4  3 £(xix3)  (26) 

^(X2X4^2  ^ ^(X3X4^3  + ^(x4^  ^4  s ^(xix4)  (27) 

For  cases  of  more  than  three  independent  variables,  the  appropriate 
set  of  simultaneous  equations  can  be  easily  constructed  after  study- 
ing the  patterns  of  the  above  two  sets  of  equations.  In  such  cases, 
solution  of  the  equations  becomes  tedious,  and  considerable  time 

■ can  be  saved  by  use  of  the  Crout  method  outlined  in  reference  2. 

Also,  programs  are  available  for  solution  of  simple  or  multiple 

r linear  regression  problems  on  practically  any  type  of  electronic 

| computer. 

k 

d.  For  multiple  regression  equations,  the  regression  constant 
should  be  determined  as  follows: 

* 

| a = M 1 3 b2M2  - b3M3  ...»  bnMn  (28) 

F o 

: e.  In  equation  20  and  equations  23  - 27,  the  quantities  t(x) 

and  T(Xlx2)  are  obtainable  rapidly  by  use  of  the  equations 

E(x)2  = E(X)2  - (£X)2/N  (29) 

; e(Xix2)  = - SX]  ex2/n  (30) 

9-03.  THE  CORRELATION  COEFFICIENT  AND  STANDARD  ERROR 

a_.  The  correlation  coefficient  is  the  square  root  of  the 

I coefficient  of  determination,  which  is  the  proportion  of  the  vari- 
ance of  the  dependent  variable  that  is  explained  by  the  regression 

equation.  A correlation  coefficient  of  1.00  would  correspond  to  a 
coefficient  of  determination  of  1.00,  which  is  the  highest  theo- 
retically possible  and  indicates  that  whenever  the  values  of  the 
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independent  variables  are  known  exactly,  the  corresponding  value  of  the 
dependent  variable  can  be  calculated  exactly.  A correlation  coefficient 
of  0.5  would  correspond  to  a coefficient  of  determination  of  0.25,  which 
would  indicate  that  25  percent  of  the  variance  is  accounted  for  and  75 
percent  unaccounted  for.  The  remaining  variance  (error  variance)  would 
bo  75  percent  of  the  original  variance  and  the  remaining  standard  error 
would  bo  the  square  root  of  this  or  07  percent  of  the  original  standard 
deviation  of  the  dependent  variable.  Thus,  with  a correlation 
coefficient  of  0.5,  the  average  error  of  estimate  would  be  37  percent 
of  the  average  errors  of  estimate  based  simply  on  the  mean  observed 
value  of  the  dependent  variable  without  a regression  analysis. 

b.  The  correlation  coefficient  (R)  Is  determined  by  use  of  the 
following  equations: 

R2  = 1 - (1  - R2) (U  - 1 )/df . (31) 


b2  ^x-jXg)  + b3  e(XjX3)  . . . + bn  r,(x^xi  ) 


zU,)' 


In  the  case  of  simple  correlation,  equation  32  resolves  to 


( I 


(33) 


c.  The  number  of  degrees  of  freedom  (df),  is  obtained  by  sub- 
tracting the  number  of  variables  (dependent  and  independent)  from  the 
number  of  events  tabulated  for  each  variable. 

d.  The  standard  error  (S«)  of  a set  of  estimates  is  the  root- 
mean-square  error  of  those  estimates.  On  the  average,  about  one  out 
of  three  estimates  will  have  errors  greater  than  the  standard  error 
and  about  one  out  of  20  will  have  errors  greater  than  twice  the 
standard  error.  The  error  variance  is  the  square  of  the  standard 
error.  The  standard  error  or  error  variance  of  estimates  based  on  a 
regression  equation  is  calculated  from  the  data  used  to  derive  the 
equation  by  use  of  either  of  the  following  equations: 


i 

1 


j 

fj 


- 47  - 





9-03 


<.2  _ £ ( X-. ) - ^2^ ^ J'  1 X2 ) **  b3^xlx3^  • • • • " (*^xg^  (34) 

e — — — — — - — — 


df 

Sg  s (1  - K2)E(x1)2/(N  - 1)  (35) 

Inasmuch  as  there  Is  some  degree  of  error  involved  In  estimating 
the  regression  coefficients,  the  actual  standard  error  of  an  esti- 
mate based  on  one  or  more  extreme  values  of  the  Independent  vari- 
ables is  somewhat  larger  than  Is  Indicated  by  the  above  equations, 
but  this  fact  is  usually  neglected. 

e,  In  addition  to  considering  the  amount  of  variance  that  is 
indicated  by  the  correlation  coefficient  and  standard  error  to  be 
solved  by  the  regression  equation,  it  is  important  to  consider  the 
reliability  of  these  indications.  There  is  some  chance  that  any 
correlation  is  accident^;,  but  the  higher  the  correlation  and  the 
larger  the  sample  upon  which  It  is  based,  the  less  is  the  chance 
that  it  would  occur  by  accident.  Also,  the  reliability  of  a regres- 
sion equation  decreases  rapidly  as  the  number  of  Independent  vari- 
ables increases,  and  extreme  care  must  be  exercised  in  the  use  of 
multiple  correlation  in  cases  based  on  small  samples. 

9-04.  SIMPLE  LINEAR  CORRELATION  EXAMPLE 

a.  An  example  of  a simple  linear  correlation  analysis  is 
illustrated  on  exhibits  31  and  32.  The  study  from  which  this  example 
was  taken  Involved  the  determination  of  the  areal  distribution  of 
short-duration  precipitation  in  * mountainous  region.  Inasmuch  as 
short-duration  measurements  were  available  at  a relatively  small 
number  of  locations,  it  was  decided  to  Investigate  the  relationship 
of  short-duration  to  long-duration  precipitation  measurements,  which 
were  available  at  many  locations. 

b^.  Inasmuch  as  long-duration  precipitation  is  made  up  of  the 
sum  of  short-duration  precipitation  amounts,  there  is  no  question 
as  to  the  existence  of  a physical  relationship,  and  It  is  therefore 
obvious  that  the  first  requirement  of  a correlation  analysis  (log- 
ical physical  relationship)  is  satisfied.  Values  of  maximum  recorded 
12-hour  precipitation  and  of  moan  annual  precipitation  were  tabulated 
as  shown  on  exhibit  31  and  plotted  as  shown  on  exhibit  32.  It  was 
determined  that  the  relation  on  logarithmic  paper  would  logically 
approximate  a straight  line.  Accordingly,  the  logarithms  of  the 
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values  were  tabulated,  and  a linear  correlation  study  of  them 
made,  as  Illustrated  on  exhibit  31,  using  equations  given  In 
paragraph  9-02.  As  the  Item  to  be  calculated  is  the  short-dura- 
tion precipitation,  the  logarithm  of  that  item  Is  selected  as  the 
dependent  variable  (Xj). 

c.  The  regression  equation  is  plotted  as  curve  A on 
exhibit  32.  This  curve  represents  the  best  estimate  of  what  the 
maximum  12-hour  precipitation  would  be  at  a location  where  the 
mean  annual  precipitation  is  known  and  the  maximum  12-hour  pre- 
cipitation is  not. 

d.  In  addition  to  the  curve  of  best  fit,  approximate  relia- 
bility-limit curves  are  established  at  a distance  of  2 standard 
errors  from  curve  A.  As  logarithms  are  used  in  the  regression 
analysis,  the  effect  of  adding  (or  subtracting)  twice  the  standard 
error  to  the  estimate  is  equivalent  to  multiplying  (or  dividing) 
the  precipitation  values  by  the  antilogarithm  of  twice  the  standard 
error.  In  this  case,  the  standard  error  is  0.081,  and  the  anti- 
logarithm  of  twice  this  quantity  is  1.45.  Hence,  values  of  12- 
hour  precipitation  represented  by  the  limit  curves  are  those  of 
curve  A multiplied  and  divided  respectively  by  1.45.  In  about 

95  percent  of  all  cases,  the  true  value  of  the  dependent  variable 
will  lie  between  these  limit  curves. 

9-05.  FACTORS  RESPONSIBLE  FOR  NON-DETERMINATION 

a.  Factors  responsible  for  correlations  being  less  than  1 
(perfect  correlation)  consist  of  pertinent  factors  not  considered 
in  the  analysis  and  of  errors  in  the  measurement  of  those  factors 
considered.  If  the  effect  of  measurement  errors  Is  appreciable, 

It  is  possible  in  some  cases  to  evaluate  the  standard  error  of 
measurement  of  each  variable  (see  paragraph  9-03d)  and  to  adjust 
the  correlation  results  for  such  effects. 

b.  If  an  appreciable  i-onion  of  the  variance  of  X,  (depend- 
ent variable)  is  attributable  to  measurement  errors,  theft  the 
regression  equation  would  be  more  reliable  than  is  Indicated  by 
the  standard  error  of  estimate  computed  from  equations  34  or  35. 
This  Is  because  the  departure  of  some  of  the  points  from  the 
regression  line  on  exhibit  32  is  artificially  Increased  by  measure- 
ment errors  and  therefore  exaggerates  the  unreliability  of  the 
regression  function.  In  such  a case,  the  curve  is  generally  closer 
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to  the  true  values  than  to  the  erroneous  observed  values.  Where 
there  Is  large  measurement  error  of  the  dependent  variable*  the 
error  of  regression  estimates  should  be  obtained  by  subtracting 
the  measurement  error  variance  from  the  error  variance  obtained 
from  equation  34  or  35.  If  well  over  half  of  the  variance  of  the 
points  from  the  best-fit  line  Is  attributable  to  measurement 
error  In  the  dependent  variable,  then  the  regression  line  would 
actually  yield  a better  estimate  of  a value  than  the  original  meas- 
urement. 


c.  If  appreciable  errors  exist  In  the  values  of  the  independ- 
ent variable,  the  regression  coefficient  and  constant  will  be 
affected,  and  fallacious  estimates  will  result.  Hence,  It  is  Impor- 
tant that  values  of  the  independent  variable  be  rather  accurately 
determined. 

d.  In  the  example  used  in  paragraph  9-C4,  there  may  well  be 
factors  responsible  for  high  short-time  Intensities  that  do  not 
contribute  appreciably  to  annual  precipitation.  Consequently,  some 
locations  with  extremely  high  mean  annual  precipitation  may  have 
maximum  short-time  Intensities  that  are  rot  correspondingly  high, 
and  vice  versa.  Therefore,  the  station  having  the  highest  mean 
annual  precipitation  would  not  automatically  have  the  highest  short- 
time  intensity,  but  would  In  general  have  something  less  than  this. 

On  the  other  hand,  if  mean  annual  precipitation  were  made  the 
dependent  variable,  the  station  having  the  highest  short-time  inten- 
sity would  be  expected  to  have  something  less  than  the  highest  value 
of  mean  annual  precipitation.  Thus,  by  Interchanging  the  variables, 
a change  In  the  regression  line  1$  effected.  Curve  B of  exhibit  32 
is  the  regression  curve  obtained  by  interchanging  the  variables  X, 
and  X~.  As  tnere  Is  a considerable  difference  In  the  two  regres- 
sion curves,  it  is  important  to  use  the  variable  whose  value  Is  to 
be  calculated  from  the  regression  equation  as  the  dependent  vari- 
able In  those  cases  where  some  Important  factors  have  not  been  con- 
sidered in  tne  analysis.  If  It  Is  obvious  chat  all  of  the  pertinent 
variables  are  Included  in  the  analysis,  then  the  variance  of  the 
points  about  the  regression  line  Is  due  entirely  to  measurement 
errors,  and  the  resulting  difference  In  slope  of  the  regression 
lines  is  entirely  artificial.  In  cases  where  all  pertinent  variables 
are  considered  and  the  great  preponderance  of  the  measurement  error 
Is  in  one  variable,  that  variable  should  be  used  as  the  dependent 
variable,  as  its  errors  will  then  not  affect  the  slope  of  the  regres- 
sion line.  In  other  cases  where  all  pertinent  variables  are 
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considered,  an  average  slope  should  be  used.  An  average  slope 
can  be  obtained  by  use  of  the  following  equation: 


EX^ 

1/2 

sd 

eX^  - (EX^)2  / N 

*7 

EX22  - (EX22  / N 

1/2 


(36) 


9-06.  MULTIPLE  LINEAR  CORRELATION  EXAMPLE 

a.  An  example  of  a multiple  linear  correlation  is  illus- 
trate? on  exhibit  33.  In  this  case,  the  volume  of  spring  runoff 
Is  correlated  with  the  water  equivalent  of  the  snow  cover  meas- 
ured on  April  1,  the  winter  low-water  flow  (index  of  ground  water) 
and  the  precipitation  falling  on  the  area  during  April.  Here 
again,  It  was  determined  that  logarithms  of  the  values  would  be 
used  in  the  regression  equation.  Although  the  loss  of  4 degrees  of 
freedom  of  12  available,  as  in  this  case,  is  not  ordinarily 
desirable,  the  correlation  attained  (0.96)  Is  particularly  high, 
and  the  equation  Is  consequently  fairly  reliable.  Note  that  the 
column  arrangement  of  the  cross-product  sums  identical  to  their  arrange- 
ment in  the  simultaneous  equations. 

In  determining  whether  logarithms  should  be  used  for 
the  dependent  variable  as  above,  questions  such  as  the  following 
should  be  considered:  "Would  an  increase  In  snow  cover  contrib- 
ute a greater  increment  to  runoff  under  conditions  of  high  ground 
water  (wet  ground  conditions)  than  under  conditions  of  low  ground 
water?"  If  the  answer  Is  yes,  then  a logarithmic  dependent 
variable  (by  which  the  effects  are  multiplied  together)  would  be 
superior  to  an  arithmetic  dependent  variable  (by  which  the  effects 
are  added  together).  Logarithms  should  be  used  for  the  independ- 
ent variables  when  they  would  increase  the  linearity  of  the 
relationship.  Whenever  logarithms  are  used,  the  logarithms  should 
be  taken  of  values  that  have  a natural  lower  limit  of  zero  and  a 
natural  upper  limit  that  is  large  compared  to  the  values  used 
In  the  study. 

c.  It  will  be  recognized  that  multiple  correlation  performs 
a function  that  is  difficult  to  perform  graphically.  Reliability 
of  the  results,  however,  Is  highly  dependent  on  the  availability 
of  a large  sampling  of  all  Important  factors  that  influence  the 
dependent  variable.  In  this  case,  the  standard  error  of  an  estimate 
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as  shown  on  exhibit  33  is  approximately  0.037,  which,  when  added 
to  a ’•ngarithm  of  a value,  Is  equivalent  to  multiplying  that  value 
by  1.09.  Thus,  the  standard  error  is  about  9 percent,  and  the  1 - In- 
20  ( rror  is  roughly  18  percent.  As  discussed  in  paragraphs  9-03e 
and  10-r  „ however,  the  calculated  correlation  coefficient  may  be 
accidentally  high.  It  can  be  demonstrated  that  the  calculated 
correlation  of  the  parameters  as  low  as  0.89  (one  chance  in  20). 

With  a correlation  coefficient  of  0.89,  and  therefore  a coefficient 
of  determination  (lr)  of  0.79,  the  standard  error  (from  equation  35) 
would  be  0.061.  The  antlloqarithm  is  1.15,  so  the  true  standard 
error  might  easily  be  almost  double  that  estimated. 

9-07 . PARTIAL  CORRELATION 

The  value  gained  by  using  any  single  variable  (such  as  April 
precipitation)  In  a regression  equation  can  be  measured  by  making 
a second  correlation  study  using  all  of  the  variables  of  the  regres- 
sion equation  except  that  one.  The  loss  in  correlation  is  expressed 
in  terms  of  the  partial  correlation  coefficient,  which  is  a measure 
of  the  decrease  in  error  attributable  to  adding  one  variable  to  the 
correlation.  The  square  of  the  partial  correlation  coefficient  is 
obtained  as  follows: 
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in  which  the  first  subscript  ahead  of  the  decimal  indicates  the 
dependent  variable  and  the  second  indicates  the  variable  whose 
partial  correlation  coefficient  is  being  computed,  and  the  sub- 
scripts after  the  decimal  indicate  the  independent  variables. 

This  procedure  is  fairly  laborious  except  where  electronic  com- 
puters are  used,  and  approximation  of  the  partial  correlation  can 
be  made  by  use  of  beta  coefficients,  which  are  very  easy  to  obtain 
by  use  of  the  following  equation  after  the  regression  equation  has 
been  calculated: 
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(38) 
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The  beta  coefficients  of  the  variables  are  proportional  to  the 
Influence  of  each  variable  on  the  result.  While  the  partial 
correlation  coefficient  measures  the  Increase  in  correlation  that 
Is  obtained  by  addlton  of  one  more  Independent  variables  to  the 
correlation  study,  the  beta  coefficient  Is  a measure  of  the  pro- 
portional Influence  of  a given  Independent  variable  on  the 
dependent  variable.  These  two  coefficients  are  related  closely 
only  when  there  Is  no  interdependence  among  the  various  "Inde- 
pendent" variables.  However,  some  "Independent"  variables 
naturally  correlate  with  each  other,  and  when  one  Is  removed 
from  the  equation,  the  other  will  take  over  some  of  Its  weight 
in  the  equation.  For  this  reason,  it  must  be  kept  In  mind  that 
beta  coefficients  indicate  partial  correlation  only  approximately. 

9-08.  VERIFICATION  OF  CORRELATION  RESULTS 

Acquisition  of  basic  data  after  a correlation  study  has  been 
completed  will  provide  an  opportunity  for  making  a check  of  the 
correlation  results.  This  Is  done  simply  by  comparing  the  values 
of  the  dependent  variable  observed,  with  corresponding  values  cal- 
culated from  the  regression  equation.  The  differences  are  the 
errors  of  estimate,  and  their  root-mean-square  Is  an  estimate  of 
the  standard  error  of  the  regression-equation  estimates  (paragraph 
9-03).  This  standard  error  can  be  compared  to  that  already  estab- 
lished In  equation  34  or  35.  If  the  difference  is  not  significant, 
there  is  no  reason  to  suspect  the  regression  equation  of  being 
Invalid,  but  If  the  difference  is  large,  the  regression  equation 
and  standard  error  should  be  recalculated  using  the  additional  data 
acquired. 

9-09.  PRACTICAL  GUIDE  LINES 

The  most  important  thing  to  remember  in  making  correlation 
studies  Is  that  accidental  correlations  occur  frequently,  par- 
ticularly when  the  number  of  observations  Is  small.  For  this 
reason,  variables  should  be  correlated  only  when  there  Is  reason 
to  believe  that  there  is  a physical  relationship.  It  is  helpful 
to  make  preliminary  examination  of  relationships  between  two  or 
more  variables  by  graphical  plotting.  This  is  particularly  helpful 
for  determining  whether  a relationship  is  linear  and  In  selecting 
a transformation  for  converting  curvilinear  relationships  to  linear 
relationships.  It  should  also  be  remembered  that  the  chance  of 
accidentally  high  correlation  Increases  with  the  number  of  correla- 
tions tried.  If  a variable  being  studied  is  tested  against  a 
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dozen  other  variables  at  random*  there  Is  a good  chance  that  one 
of  these  will  produce  a good  correlation,  even  though  there  may 
be  no  physical  relation  between  the  two.  In  general,  the  results 
of  correlation  analyses  should  be  examined  to  assure  that  the 
derived  relationship  Is  reasonable.  For  example.  If  streamflow 
Is  correlated  with  precipitation  and  drainage  area  size,  and  the 
regression  equation  relates  streamflow  to  some  power  of  the  drainage 
area  greater  than  one,  a maximum  exponent  value  of  one  should  be 
used,  because  the  flow  per  square  mile  cannot  increase  with  drainage 
area  when  other  factors  remain  constant. 
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SECTION  10  - STATISTICAL  RELIABILITY  CRITERIA 
10-01.  FUNCTION 

One  of  the  principal  advantages  of  the  use  of  statistical 
procedures  Is  that  they  provide  means  for  evaluating  the  relia- 
bility of  the  estimates.  This  permits  a broader  understanding 
of  the  subject  and  provides  criteria  for  decision  making.  The 
common  statistical  Index  of  reliability  Is  the  standard  error 
of  estimate,  which  Is  defined  as  the  root-mean-square  error. 

In  general,  It  Is  considered  that  the  standard  error  is  exceeded 
on  the  positive  side  one  time  out  of  six  estimates,  and  equally 
frequently  on  the  negative  side,  for  a total  of  one  time  In 
three  estimates.  An  error  twice  as  large  as  the  standard  error 
of  estimate  Is  considered  to  be  exceeded  one  time  In  40  in  either 
direction,  for  a total  of  one  time  In  20.  These  are  only  approxi- 
mate frequencies  of  errors,  and  exact  statements  as  to  error 
probability  must  be  based  on  examination  of  the  frequency  curve 
of  errors  or  the  distribution  of  the  errors.  Both  the  standard 
error  of  estimate  and  the  distribution  of  errors  will  be  dis- 
cussed In  this  section. 

10-02.  RELIABILITY  OF  FREQUENCY  STATISTICS 

a.  The  standard  errors  of  estimate  of  the  mean,  standard 
deviation,  and  skew  coefficient,  which  are  the  principal  statistics 
used  In  frequency  analysis,  are  given  in  the  following  equations: 


s„  - / SS/N 

(39) 

Ss  « / S“V2N 

(40) 

/6N(N-"I  )/(N-2)  (N+l  )(N+3) 

(41) 

These  have  been  used  to  considerable  advantage  as  discussed  in 
paragraph  7-09  in  drawing  maps  of  mean  and  standard  deviation  for 
a regional  frequency  study. 

b.  The  distribution  of  errors  of  estimating  tht  mean  Is 
a function  of  the  t distribution,  exhibit  34,  and  is  given  by  the 
following  equation: 
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In  any  specific  problem,  values  of  N,  M,  and  S are  obtained  from 
the  recorded  data.  To  find  a value  of  M-y  (difference  between 
computed  and  true  mean)  that  Is  exceeded  with  a specified  proba- 
bility, ascertain  from  exhibit  34  the  value  of  t corresponding  to 
that  probability,  equate  this  to  C and  solve  for  M-y  by  replacing 
the  Inequality  sign  by  an  equal  sign  In  the  left  bracket  of  equa- 
tion 42.  The  t distribution  Is  symmetrical  and  therefore  the 
probability  that  M-u  exceeds  a specified  value  is  equal  to  the 
probability  that  y-M  exceeds  that  same  value,  and  therefore  the 
probability  that  the  absolute  value  of  M-y  is  exceeded  regardless 
of  direction  Is  twice  as  great  as  Indicated  in  exhibit  34.  (Most 
published  tables  of  t show  twice  the  probabilities  Indicated  in 
exhibit  34  and  therefore  represent  both  tails  of  the  t distri- 
bution.) 


c.  The  distribution  of  errors  In  estimating  the  standard 
deviation  is  a function  of  the  chi-square  distribution,  exhibit 
35,  and  is  given  by  the  following  equation: 


Prob 


Prob 


(43) 


Here  again  N and  S are  obtainable  from  recorded  data.  To  find  a 
value  of  a (true  standard  deviation)  that  Is  exceeded  with  a 2 
specified  probability,  ascertain  from  exhibit  35  the  value  of  x 
(chi-square)  that  is  exceeded  with  that  probability,  and  compute 
a by  replacing  the  inequality  sign  in  the  left  bracket  of  equation 
43  by  an  equal  sign.  (Note  that  equation  43  is  also  valid  If  both 
Inequality  signs  are  reversed  simultaneously.) 

10-03.  RELIABILITY  OF  FREQUENCY  ESTIMATES 

The  reliability  of  analytical  frequency  determinations  can 
best  be  illustrated  by  establishing  error-limit  curves.  The 
error  of  the  estimated  flow  for  any  given  frequency  is  a function 
of  the  errors  in  estimating  the  mean  and  standard  deviation,  assum- 
ing that  the  skew  coefficient  is  known.  Criteria  for  construction 
of  error-limit  curves  are  based  on  the  distribution  of  the  "non- 
central t".  Selected  values  transformed  for  convenient  use  are 
given  on  exhibit  6.  By  use  of  this  exhibit  and  equation  5,  error- 
limit  curves  shown  on  exhibit  7 were  calculated  as  illustrated  on 
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that  exhibit.  While  the  expected  frequency  is  that  shown  by  the 
middle  curve,  there  is  one  chance  in  20  that  the  true  value  for 
any  given  frequency  is  greater  than  that  indicated  by  the  .05 
curve  and  one  chance  in  20  that  it  is  smaller  than  the  value 
Indicated  by  the  .95  curve.  There  are,  therefore,  nine  chances 
in  ten  that  the  true  value  lies  between  the  .05  and  .95  curves. 

10-04.  RELIABILITY  OF  CORRELATION  RESULTS 

The  reliability  of  correlation  results  has  been  discussed 
to  some  extent  in  section  9.  It  has  been  demonstrated  that  an 
estimate  based  on  a regression  equation  has  a standard  error  of 
estimate  expressed  by  equation  35.  The  advantage  gained  in 
making  regression  analysis  is  measured  by  the  determination 
coefficient  or  correlation  coefficient  computed  by  use  of  equa- 
tion 31.  There  are  further  reliability  criteria  that  are  useful 
in  making  decisions  in  correlation  studies.  To  determine  whether 
it  is  worthwhile  to  Introduce  a particular  variable  in  an  analysis, 
compute  the  partial  correlation  coefficient  as  explained  in  para- 
graph 9-07.  Even  though  the  partial  correlation  coefficient  is 
positive,  however,  it  is  possible  that  it  is  accidentally  high 
and  therefore  not  dependable.  The  chance  that  a partial  correla- 
tion coefficient  is  accidental  (that  the  true  coefficient  is  zero) 
can  be  computed  by  use  of  the  following  equation: 


1 


r 

■ 

Prob 

[jyH 

= Prob 

t > c 

(44) 

This  equation  is  also  valid  when  applied  to  a simple  correlation 
coefficient,  where  df  ■ N-2, 
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SECTION  11  - TERMS  AND  SYMBOLS 


TERMS 

Analytical  method  - Method  of  fitting  frequency  curves  by  moments 
(par  2**05c) 

Annual  event  - The  largest  (or  smallest)  event  In  the  year 
Cumulative  frequency  curve  - Relation  of  magnitude  to  percentage 
of  events  exceeding  that  magnitude  (par  2-04) 

Deviation  - Difference  between  an  Individual  magnitude  and  the 
average  for  the  frequency  array  (par  4-02a) 

Distribution  - Function  describing  the  relative  frequency  with  which 
events  of  various  magnitudes  occur  (par  2-05c) 

Duration  curve  - Curve  Indicating  the  percentage  of  time  that 

various  rates  of  flow  are  exceeded  at  a specified  river  station 
(par  2-04e) 

Error  variance  - Square  of  the  standard  error  (par  9~03d) 

Exceedence  interval  - The  average  interval  of  time  between  values 
that  exceed  a specified  magnitude;  reciprocal  of  the  exceedence 
frequency 

Exceedence  frequency  - The  percentage  of  values  that  exceed  a 
specified  magnitude  (par  2-04a) 

Exceedence  probability  - Probability  that  an  event  selected  at 
random  will  exceed  a specified  magnitude  (par  2-04a) 

Expected  probability  - Statistical  average  of  estimated  future 
probabilities  (par  2-05c) 

Frequency  array  - List  of  events  arranged  in  the  order  of  magnitude 
(par  2-05b) 

Frequency  curve  - Graphical  representation  of  a frequency  distribu- 
tion, usually  with  the  ab^cl  sa  as  magnitude  and  the  ordinate 
as  relative  frequency,  but  aiso  used  Interchangeably  with 
cumulative  frequency  curve. 

Geometric  mean  - Antllogarithm  of  the  average  logarithm 
Logarithmic  normal  grid  - Grid  on  which  a cumulative  frequency  curve 
of  a variable  whose  logarithms  are  normally  distributed  will 
plot  as  a straight  line,  (par  3-06) 

Logarithmic  probability  grid  - Same  as  logarithmic  normal  grid 
Maximum  likelihood  - A designation  applied  to  a statistical  estimate 
that  is  more  likely  than  any  other  estimate  to  be  true  on  the 
basis  of  data  employed  in  its  estimation  (par  4-02c) 

Normal  distribution  - Ideal  frequency  distribution  approached  by 
many  observed  distributions  which  are  symmetrical  and  in  which 
the  small  deviations  greatly  outnumber  the  large  ones,  the 
possible  deviations  being  virtually  unlimited  in  magnitude  by 
physical  conditions. 


TERMS 


Parent  population  - Sum  total  of  events  that  will  occur  In  the 
future  If  pertinent  existing  controls  continue  Indefinitely 
(par  2-03) 

Partial -duration  curve  - Cumulative  frequency  curve  of  all  events 
above  a base  value  (par  2-04c) 

Pearson  Type  III  function  - Family  of  asymmetrical , 

ideal  frequency  distributions  of  which  the  normal  distribu- 
tion Is  a special  case  (ex  39) 

Plotting  position  - Exceedence  probability  of  a magnitude,  estimated 
from  its  position  In  a frequency  array  (par  3-05) 

Probability  grid  - Grid  on  which  a cumulative  frequency  curve  of 
a normal  distributed  variable  will  plot  as  a straight  line 
(par  3-06) 

Random  - Condition  under  which  events  occur  if  magnitudes  of  suc- 
cessive events  are  not  correlated. 

Recurrence  interval  - Exceedence  interval,  In  the  case  of  floods 
or  storms 

Return  period  - Recurrence  interval 

Standard  error  - Root-mean-square  error  (par  9-03) 

Standard  deviation  - Root-mean-square  deviation  (par  4-02) 

Skew  coefficient  - Function  of  the  third  moment  of  magnitudes 
about  their  mean,  a measure  of  asymmetry  (par  4-02) 

Variance  - Mean  square  deviation;  square  of  the  standard  deviation 
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SYMBOLS 

a Reqression-equation  constarf  (par  9-02a) 

A Drainage  area  size 

b Regression-equation  coefficient  (par  9-02a) 

0 Beta  coefficient  (par  9-07) 

C Any  constant 
D Depth  of  rainfall 
df  Degrees  of  freedom  (par  9-03) 
g Skew  coefficient  (par  4-03c) 

k Coefficient  of  S (eq  5) 

m Order  number  of  event  (par  3-05b) 

M Computed  mean 
u True  mean 

N Number  of  events  or  number  of  years  of  record 
P Probability,  specifically  exceedence  frequency  per  hundred  years 
P«  Maximum-likelihood  exceedence  probability  (par  4-03d) 

Q Flow  or  runoff 

R Unadjusted  correlation  coefficient  (par  9-03) 

R Correlation  coefficient,  square  root  of  determination  coef- 
ficient (par  9-03) 

r Partial -correlation  coefficient  (par  9-07) 

S Computed  standard  deviation  (par  4-03b) 
o True  standard  deviation 
$e  Standard  error  (par  9-03) 

T Duration  of  runoff 
t A theoretical  function  (ex  34) 
x Deviation  of  X from  (par  4-02a) 

X Any  variable,  the  dependent  variable  when  subscript  Is  1,  but 

usually  an  independent  variable;  frequently  the  logarithm  of  Q 

y Deviation  of  Y from  My  (par  4-02a) 

Y Any  variable,  usually  the  dependent  variable 
2 

X Chi-square,  a theoretical  variable  (ex  35) 
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APPENDIX  I - KEY  QUESTIONS  AND  ANSWERS 

The  following  questions  include  general  questions  that  might 
normally  occur  to  a person  after  completing  this  paper  and  are 
intended  to  Insure  that  the  main  points  in  the  paper  have  been 
properly  understood. 

1.  What  are  the  principal  advantages  of  frequency  analysis 
over  a simple  assumption  that  a record  will  repeat  Itself? 

First,  a properly  conducted  frequency  study  is  more 
accurately  representative  of  future  probabilities,  and  secondly, 
a frequency  study  can  be  used  to  represent  periods  of  varying 
length  regardless  of  the  period  of  record  and  various  geographic 
locations  regardless  of  the  location  of  the  record. 

2.  How  can  we  tell  that  one  method  of  hydrologic  frequency 
analysis  yields  more  accurate  results  than  another,  since  it  would 
take  hundreds  or  thousands  of  years  to  prove  or  disprove  an  esti- 
mate? 


Reference  24  demonstrates  a method  of  dividing 
long  records  into  parts  so  that  estimates  based  on  one  part  are 
judged  by  the  remaining  (unused)  portion  of  the  record.  When 
this  is  done  for  alternative  methods  using  a great  many  stations, 
one  is  able  to  select  as  more  dependable  the  method  that  yields 
the  best  results  on  the  average. 

3.  Why  not  design  spillways  for  major  dams  to  pass  the 
1000-year  flood  only? 

First,  the  1000-year  flood  cannot  be  estimated 
dependably--there  have  already  been  cases  where  the  1000-year 
flood  estimate  has  subsequently  been  exceeded.  Secondly,  one 
structure  out  of  10  or  20  would  fail  during  their  economic  life- 
time and  many  more  during  the  expected  physical  lifetime. 

4.  If  there  were  1,000  independent  10-year  records  in  the 
country,  how  many  should  you  expect  to  contain  one  or  more  floods 
larger  than  the  true  100-year  magnitude? 

About  90 

5.  What  ic  an  objection  to  using  the  terms  50-year  flood, 
100-year  flood,  etc.? 

To  some  it  implies  that  a 50-year  flood  will  not 
be  exceeded  in  50  years  whereas  it  can  be  exceeded  next  year.  To 
others  it  implies  that  it  is  a high  degree  of  protection,  whereas 
a structure  designed  to  be  safe  against  a 50-year  flood  will  have 
a 64  percent  probability  of  failure  in  a 50-year  period,  or  an 
18  percent  chance  of  failure  in  a 10-year  period. 
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6.  Why  Is  the  largest  flood  In  a 10-year  period  not  equal 
to  the  10-year  flood,  on  the  average? 

Because  the  10-year  flood  is  defined  as  one  that 
is  exceeded  on  the  average  once  in  10  years,  not  equalled.  There- 
fore the  10-year  flood  should  be  smaller  than  the  largest  flood 
in  a 10-year  period,  on  the  average. 

7.  What  are  the  main  factors  affecting  the  accuracy  of 
ordinary  frequency  determinations? 

The  length  of  record  and  variance  of  runoff.  Com- 
pared to  the  effects  of  random  occurrence,  the  accuracy  of  measuring 
events  that  do  occur  is  a minor  consideration  in  ordinary  cases. 

8.  Why  are  extreme- frequency  estimates  far  less  reliable 
than  ordinary- frequency  estimates? 

Because  factors  that  are  possible  but  of  rare  occur- 
rence such  as  a rare  type  of  storm,  changes  in  river  course, 
landslide-caused  floods,  etc.,  may  not  be  properly  reflected  •> 
the  record. 

9.  If  we  were  willing  to  accept  an  accuracy  such  that  only 
one  case  out  of  20  would  exceed  our  estimate  plus  allowable  error, 
how  much  would  our  error  allowance  be  for  estimating  the  2-year  and 
100-year  floods  from  a 40-year  record? 

About  20  percent  and  40  percent,  respectively,  where 
the  stream  is  moderately  variable.  About  1/3  of  these  values  for 
stable  streams  like  the  Mississippi  and  far  greater  for  erratic 
streams. 


10.  What  does  a simple  correlation  coefficient  or  partial 
correlation  coefficient  of  0.5  imply? 

That  the  Independent  variable  being  tested  reduces 
the  error  of  estimate  by  a fraction  represented  by  one  minus  the 
square  root  of  one  minus  the  square  of  the  correlation  coefficient, 
or  about  13  percent,  in  this  case. 
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APPENDIX  II  - PROBLEMS  AND  SOLUTIONS 

The  following  problems  can  be  used  by  students  as  exercises 
necessary  to  a thorough  understanding  of  the  principles  and  pro- 
cedures contained  in  this  paper.  Solutions  should  be  read  only 
after  the  problem  has  been  worked,  or  to  the  extent  necessary  to 
clarify  any  misunderstanding. 

1.  What  Is  the  chance  that  exactly  three  50-year  floods 
will  occur  In  a 100-year  period?  What  Is  the  chance  that  three 
or  more  will  occur? 

This  question  needs  interpretation.  Since  the 
50-year  flood  Is  an  exact  figure,  the  chance  that  It  will  occur 
exactly  1$  negligible,  or  zero.  We  will,  however,  interpret 
this  question  to  mean  exceeded.  The  general  formula  for  the 
exact  number  of  chance  events,  1,  out  of  N trials  is: 


P - 


N! 

TTTfpryr 


^o-p) 


N-i 


In  which  P Is  the  probability  of  obtaining  In  N trials  exactly 
1 events  having  a probability  p of  occurring  In  a single  trial. 
Substituting  100  for  N,  3 f or  i , and  .02  for  p,  we  obtain 
P=. 183,  or  about  one  chance  in  five  or  six. 

To  obtain  the  answer  to  the  second  question,  we 
could  substitute,  3,  4,  5,  etc.  up  to  100  in  the  formula  and  add 
the  probabilities  of  these  mutually  exclusive  (if  one  occurs, 
the  others  can't)  events,  or  we  can  shorten  the  work  by  substi- 
tuting 2,  1 and  0 In  the  formula  and  subtracting  the  sum  of  proba- 
bilities from  1.00  (certainty),  since  it  is  certain  that  either 
(a)  3 or  more  or  (b)  2 or  less  must  occur.  Substituting  2,  1 
and  C in  the  formula,  we  obtain  .274,  .271,  .271,  and  ,133  in  turn  for 
P.  One  minus  their  sum  is  .322,  or  about  one  chance  in  three. 

It  is  necessary  to  know  for  this  solution  that  zero  factorial 
equals  1.00. 

2.  What  is  the  chance  that  a 1000-year  flood  will  be 
exceeded  in  the  50-year  economic  lifetime  of  a project? 

We  interpret  this  to  include  the  chance  that  the 
flood  will  be  exceeded  more  than  once, as  well  as  exactly  once. 

We  should,  therefore,  solve  for  the  chance  that  it  will  not  occur 
(i=o)  and  subtract  that  probability  for  1.00,  as  explained  in 
problem  1.  Note  that  where  1 = o,  the  equation  resolves  to 


P = O-p) 


N 


II  - 1 


APP.  II 


Substituting  .001  for  p and  50  for  N,  we  obtain  .952  for  P and, 
therefore,  .048  for  the  answer,  about  one  chance  In  20. 

3.  Using  the  peak  flows  tabulated  in  problem  5,  compute  the 
mean,  median,  mode,  variance,  standard  deviation,  skew  coefficient 
and  standard  errors  of  the  mean,  standard  deviation  and  skew  coef- 
ficient for  the  flows  and  not  for  the  logarithm.  What  Is  the  chance 
that  the  true  mean  Is  greater  than  6,000  c.f.s.?  That  the  true 
standard  deviation  Is  greater  than  the  5,000  c.f.s.? 

The  mean,  obtained  from  equation  2,  Is  8.04.  The 
median  (middle  magnitude)  Is  taken  as  the  average  of  the  two  middle 
values  of  5.47  and  7.48,  since  there  are  an  even  number  of  values, 
and  equals  6.48.  The  mode  (most  frequent  value)  should  be  estimated 
from  a plotting  as  follows: 


The  variance,  obtained  from  equation  3,  Is  41.0.  Its  square  root 
(S)  Is  the  standard  deviation,  or  6.40.  The  skew  coefficient, 
obtained  from  equation  4,  Is  0.96.  The  standard  error  of  the 
mean,  obtained  from  equation  39,  Is  1.43.  The  standard  error  of 
the  standard  devlaltlon  Is  1.01,  from  equation  40.  The  standard 
error  of  the  skew  coefficient,  obtained  from  equation  41,  Is  0.26. 
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According  to  equation  42,  the  chance  that  the  true  mean  exceeds 

6.00  Is  equal  to  the  chance  that  a value  of  t equal  to  2.04  (4.36)/ 
6.40  * 1.39  Is  exceeded  In  the  positive  direction.  This  Is  0.91, 
as  determined  from  exhibit  34.  Likewise,  according  to  equation 
43,  the  chance  that  the  true  standard  deviation  Is  greater  than 

5.000  c.f.s.  Is  equal  to  the  chance  that  a value  of  equal  to 
19(6.40)z/5.002  3 31.1  Is  not  exceeded.  This  is  0.96,  as  deter- 
mined from  exhibit  35. 

4.  Using  the  table  of  random  numbers  on  the  following 
page,  generate  twenty  5-event  random  samples  from  a normal  parent 
population  whose  mean  Is  5.00  and  standard  deviation  1.00.  Com- 
pute the  true  exceedence  probability  of  the  mean  plus  2.18  standard 
deviations  for  each  sample.  Average  the  20  probabilities. 

This  problem  Is  designed  to  teach  speed  In  analytical 
computations  and  to  demonstrate  the  nature  of  random  samples  and 
of  statistical  inference.  First,  draw  a straight-line  frequency 
curve  on  arithmetic  probability  paper  with  the  given  mean  and 
standard  deviation  as  demonstrated  on  page  1 1-5.  This  represents 
the  parent  population.  Using  the  Monte  Carlo  Method,  in  order  to 
obtain  a random  magnitude,  enter  the  curve  with  an  exceedence 
frequency  equaling  the  random  number  whose  first  two  digits 
(Including  zero)  are  percentages  and  remaining  digits  decimals 
of  a percent.  A random  number  should  be  selected  using  a pre- 
determined pattern  such  as  starting  with  the  upper  right  corner 
and  working  down  each  column.  Four  digits  are  sufficient  except 
that  a pre-determlned  procedure  is  required  for  obtaining  addi- 
tional digits  if  the  first  two  or  more  are  zeros  or  nines. 

For  each  group  of  5 random  magnitudes,  compute  the 
mean  and  standard  deviation  using  equations  2 and  3.  Multiply 
the  computed  standard  deviation  by  2.18  and  add  to  the  computed 
mean  to  obtain  an  estimate  of  the  20-year  event.  Enter  the  true 
frequency  curve  with  this  magnitude  to  obtain  its  true  exceedence 
frequency.  See  If  the  20  true  exceedence  frequencies  determined 
from  20  different  5-event  random  samples  average  about  5 percent, 
as  exhibit  38  indicates  that  they  should. 

5.  Using  data  on  page  II-6,  compute  volume  frequency  curves 
and  plot  the  curves  and  data  as  Illustrated  on  exhibits  15  to  20. 
Check  to  make  sure  that  the  plotted  points  and  computed  curves  are 
reasonably  consistent. 
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Average  flows  in  thousand  c.f.s. 


Year 

Peak 

1936-37 

2.18 

1.94 

1.79 

1937-38 

16.5 

9.70 

7.40 

-39 

1.55 

1.04 

0.63 

-40 

14.5 

10.8 

7.80 

-41 

13.2 

7.86 

6.20 

-42 

8.12 

4.65 

3.43 

-43 

18.6 

10.2 

5.93 

-44 

1.16 

0.82 

0.64 

-45 

2.31 

1.48 

1.23 

-46 

8.99 

5.76 

5.30 

-47 

5.21 

3.20 

1.66 

-48 

5.47 

3.39 

2.66 

-49 

1.66 

1.48 

1.33 

1949-50 

3.04 

2.13 

1,51 

-51 

7.48 

4.39 

3.40 

-52 

8.86 

5.36 

3.63 

-53 

10.8 

7.75 

4.17 

-54 

5.39 

4.21 

2.71 

-55 

2.19 

1.50 

0.72 

-56 

23.5 

16.3 

11.1 

10- day 

30-day 

90-day 

1 year 

1.13 

0.71 

0.64 

0.16 

2.84 

1.58 

1.42 

0.62 

0.60 

0.42 

0.22 

.089 

3.40 

1.77 

1.21 

0.39 

3.05 

1.64 

1.44 

0.60 

2.17 

1.43 

0.99 

0.40 

2.42 

1.19 

0.78 

0.28 

0.44 

0.31 

0.25 

.092 

0.98 

0.59 

0.39 

0.17 

2.67 

1.41 

0.67 

0.28 

0.73 

0.46 

0.33 

0.11 

1.38 

0.67 

0.47 

0.19 

0.91 

0.79 

0.59 

0.18 

1.13 

0.76 

0.52 

0.18 

1.96 

1.40 

0.76 

0.30 

2.00 

1.36 

1.06 

0.44 

2.95 

1.75 

0.84 

0.35 

1.56 

1.23 

0.89 

0.33 

0.56 

0.39 

0.25 

0.13 

5.47 

2.89 

1.52 

0.60 

6.  Using  the  data  on  the  following  page,  compute  the  multiple  linear 
regression  equation,  determination  coefficient,  standard  error  of 
estimate,  and  the  partial  correlation  coefficient  of  the  variable 
having  the  smallest  beta  coefficient.  Using  the  regression  equation, 
compute  the  residual  (regression  constant)  for  each  station,  plot  in 
the  drainage  area  on  the  following  map,  and  draw  lines  of  equal  resid- 
ual . 

Extensions  required  for  computing  the  regression 
equation  are  computed  as  illustrated  on  exhibit  18.  The  equation 
coefficients  and  constant  are  computed  cost  easily  by  use  of  the 
Crout's  Method  (paragraph  9-02c).  The  equation  should  be: 
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Station 

mu&er 

toon  log  of 
pooh  flow  in  o.f.s. 

Log  Of 

drainage  tree 
in  oq.  d1  . 

Log  of 

normal  annual 
prealp.  in  inahee 

ELev.  above 
enovllna 

In  hundred  ft 

14 

5.08 

3.97 

1.71 

54 

15 

3.72 

2.27 

1.62 

55 

1 6 

3.71 

2.30 

1 .66 

54 

it 

4.41 

2.76 

1.58 

17 

16 

2.23 

1.30 

1.58 

49 

4l 

5.00 

3.57 

1.72 

46 

42 

3.33 

1.89 

1.60 

49 

43 

2.46 

1.62 

1.79 

71 

44 

4.18 

2.23 

1.94 

42 

45 

3.64 

2.52 

1.61 

60 

55 

3.19 

£.21 

1.74 

72 

56 

2.83 

1.83 

1.64 

46 

5T 

3-55 

2.26 

1.67 

59 

58 

4.27 

3.13 

1.60 

61 

59 

4.6l 

3.28 

1.72 

55 

60 

4.17 

2.79 

1.80 

57 

61 

3.43 

2.25 

1.77 

63 

62 

4.20 

2.87 

I.81 

55 

69 

3.98 

2.95 

1.67 

58 

70 

3-87 

2.80 

1.64 

54 

71 

4.27 

2.63 

1.65 

34 

72 

3.78 

2.56 

1.63 

51 

73 

3.45 

1-97 

1.50 

27 

74 

3.61 

2.15 

1.51 

36 

75 

3-72 

2.17 

1.78 

46 

76 

3-19 

1.48 

I.60 

22 

77 

3-81 

2.30 

1.77 

54 

73 

3.58 

1.77 

1 .64 

20 

79 

4.51 

3.08 

1.78 

54 

80 

3.61 

1.85 

1.71 

27 

81 

3.43 

2.00 

l'.Sl 

17 

124 

2.89 

1.31 

1.38 

12 

125 

2.75 

1.68 

1.28 

11 

126 

2.61 

1.36 

1.78 

81 

127 

4.21 

2.63 

1.80 

68 

128 

2.44 

2.09 

1.62 

68 

129 

3-73 

1.81 

1.86 

44 

130 

3-25 

2.58 

1.75 

73 

131 

3-79 

2.36 

1.75 

42 

132 

4.26 

2.98 

1.65 

22 

133 

3.61 

2.09 

1.59 

29 

134 

3-64 

2.13 

1.67 

37 

135 

3-49 

1.83 

1.72 

28 

136 

4.70 

3.56 

1.68 

61 

137 

3.21 

1.54 

1.79 

43 

138 

3.89 

2.39 

1.T9 

67 

139 

4.22 

2.68 

1.81 

55 

140 

3.31 

1.71 

1.70 

77 

141 

4.0  6 

2.54 

1.76 

48 

142 

2.46 

1.36 

1.69 

65 

143 

3.56 

2.00 

1.54 

24 

144 

3.76 

2.30 

1.61 

35 

145 

3.96 

2.05 

1.73 

18 

146 

3.64 

1.91 

I.61 

16 

147 

3.85 

1.94 

1.74 

16 

148 

5.16 

3.49 

1.76 

36 

149 

4.62 

2.73 

1.65 

21 

150 

4.23 

2.30 

1.86 

42 

151 

4.70 

3.45 

1.76 

52 

152 

4.85 

2.79 

1-96 

37 

180 

4.23 

2. 65 

1-TT 

51 

l8l 

2.8l 

1.60 

1.81 

72 

182 

3-69 

2.47 

1-59 

45 

323 

2.12 

1.80 

1.41 

87 

324 

2.48 

2.26 

1.49 

92 

325 

2.13 

1.30 

1.54 

89 

326 

1.73 

1.15 

1.57 

90 

327 

2.32 

1.82 

1.51 

90 

330 

2.49 

1.52 

1.58 

75 

331 

3-03 

2.29 

1.08 

72 

332 

2.93 

2.40 

1.15 

76 

333 

2.84 

2.44 

I.32 

73 
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X1  = 0.945  X2  + 1.37  X3  - 0.014  X4  - 0.147 

The  determination  coefficient,  computed  by  use  of  equations  31  and 
32  is  0.925.  The  standard  error  of  estimate,  from  equation  35,  is 
0.21.  Beta  coefficients,  from  equation  38,  are  0.75,  0.29,  and 
-0.38  in  turn.  In  order  to  obtain  the  partial  correlation  coef- 
ficient for  the  second  variable,  it  is  necessary  to  repeat  Grout's 
Method  using  all  but  the  second  variable.  The  new  determination 
coefficient  is  0.846,  and  the  partial  correlation  coefficient, 
from  equation  37,  is  0.51.  The  residual  for  each  station  is  com- 
puted by  transposing  the  regression  equation  as  follows  and  sub- 
stituting observed  values  of  X2,  X3  and  X4  for  each  station  in  turn: 

C = X1  - 0.945  X2  - 1.37  X3  + 0.014  X4 

7.  Adjust  the  statistics  for  1-day  flows  in  problem  5 to 
take  advantage  of  the  longer  record  on  exhibit  15  and  the  further 
adjustment  of  statistics  for  that  station  on  exhibit  18. 

We  note  on  exhibit  15  that  there  are  an  additional 
10  years  of  1-day  flow  record  on  Mill  Creek  that  are  not  contained 
in  problem  5.  Also,  the  Mill  Creek  statistics  were  adjusted  by 
use  of  Feather  River  flows  recorded  over  a 47-year  period  (exhibit 
18)  to  obtain  the  equivalent  of  an  additional  13  years,  or  43 
years  in  all. 

Sy  pairino  the  20  years  of  concurrent  record,  the 
determination  coefficient  is  found  by  use  of  equations  7 and  8 to 
be  0.70,  and  other  statistics  are  simultaneously  computed  as 
illustrated  for  Mill  Creek  on  exhibit  18.  These  turn  out  to  be: 


Mill  Creek 

mean 

Ion 

20  vears 
3,474 

43  years 
3.TTTT 

Problem  5 

1. 

II 

3.578 

Mill  Creek 

std 

dev 

.317 

.287 

Problem  5 

ii 

ii 

.372 

By  use  of  equations  9 to  11,  adjusted  statistics  and  the  length  of 
equivalent  record  are  found  to  be: 

Mean  log  3.524 

Standard  deviation  .347 

Equivalent  record  36 

8.  Using  exhibits  24  to  27,  compute  a regionalized  estimate 
of  logarithmic  mean  and  standard  deviation  of  peak  flows  for  Mill 
Creek.  Compare  these  with  the  computed  values  on  exhibit  18  and 
select  a value  for  each  that  you  would  recommend  for  adoption. 
Required  data  are  as  follows: 


i 

I 

! 
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Drainage  area  134  sq.  mi. 

Basin  normal  annual  precip.  47  in. 

Average  basin  elevation  2900  ft.  m.s.l. 

This  example  is  already  worked  out  in  paragraph  7-12. 
The  standard  deviation  in  hundredths  Is  a basin-mean  value  estimated 
from  exhibit  25  and  is  0.31.  The  basin-mean  value  of  C determined 
from  exhibit  24  for  entry  into  equation  11  of  referencep21 , also 
given  in  paragraph  7-12,  is  42.  A value  of  K for  use  in  that  equa- 
tion, obtained  from  exhibit.  26,  is  0.71.  The  computed  Q is  4220 
c.f.s.  and  the  corresponding  logarithm  is  3.625.  This  compares  to 
3.661  cn  exhibit  13,  and  the  regional  standard  deviation  of  0.31 
compares  to  0.281  on  exhibit  18.  Considering  that  the  equivalent 
length  of  record  on  exhibit  18  is  43  years,  the  standard  error  of 
the  calculated  mean  from  exhibit  23  is  .152  (0.281)  or  0.43  and 
the  standard  error  of  the  calculated  standard  deviation  is  antilog 
.047,  or  a factor  of  1.11.  Since  the  regional  values  are  less  than 
one  standard  error  from  the  computed  values,  and  since  computed 
values  can  easily  be  off  by  one  standard  error,  it  is  best  to  adopt 
the  regional  values  for  the  sake  of  consistency  with  adjoining 
basins. 
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EXHIBIT  I 


ILLUSTRAJltfE  EXAHPU 

AMALYTlCAL  FREQUENCY  CCKPUTATIOH 
US  ISC  PRE -RECORD  DATA 


Location!  U i 1 Larao  ttc  H ct  Albany,  Orator. 
I.rlotl  of  recordi  1893-1958  £69  years) 


LAJ  UM  tiUlll  " * * - ..  ••  - 7"  ‘ ’ V 

Period  of  racordi  1895-1998  £69  years)  , 

Years  of  historical  flood  estinatesi  1861,  1887 , 1890  (5  largest  Known) 
Period  covered  hy  history  nnd  recordi  1858-1938)000  years) 


(Se«  por.  4-05) 


k - 6,00 

H-69 

CF 


5.86 

1.98 

5.91 

1.97 

5.77 

1.97 

3.73 

1.96 

3.69 

1.96 

5.65 

1.96 

3.60 

1.96 

5.56 

1.96 

3.52 

1.91 

3.6? 

1.91 

3.63 

1.90 

5.38 

1.89 

3.36 

1.88 

3.29 

1.88 

5.26 

1.68 

3.16 

1.88 

5.15 

1.86 

5.07 

1.85 

5.01 

1.85 

2.95 

1.06 

2.87 

1.78 

2.80 

1.77 

2.71 

1.76 

2.62 

1.73 

2.52 

1.72 

2.59 

1.72 

2.23 

1.69 

2.00 

1.67 

1.57 

1.61 

68 

66 

277-67 

139-73 

6.080 

2.053 

EX2 

<EX)2/H" 
dlff . 


1206.2039  290.1317 
1132.2000  267.1 266 
76.0039  ?.0271 


S = b - (3*0271/76. 00 39/*  - .202  (Eq.  36)  i 

M = a . 2.093  - .202(6.080-6.000)  = 2.039  (Eq.  21) 


Computation  of  Frequency  Curve  (H  = 65) 


(u> 

(12) 

(13) 

(16) 

(15) 

(16) 

(17) 

(18) 

P, 

0.25 

L.O 

10 

50 

90 

99 

99.75 

k (Sx.  38) 

2.96 

2.61 

1.31 

0.00 

-1.31 

.2.61 

-2.96 

Log  Q (Bn.  5) 

2-633 

2.526 

2.306 

2.039 

1.776 

1.552 

1.665 

Q.  thouG  cfa 

630 

336 

201 

109 

59.6 

33.6 

27-9 

Events  nur.bered  4 to  J could  not  have  occupied  higher  position*  if  the 
entire  100-year  record  were  available,  (but  might  occupy  lower  positions). 
Accordingly!  k /aloes  were  selected  as  shown  above. 

4.00  was  aided  to  ail  k values  in  order  t-tart  all  numbers  are  positive 
and  have  3 digits  for  sir..pliclty  of  machine  operation. 


ILLUSTRATIVE  EXAMPLE 

TABULATION  OF  PEAK  FLOW  FREQUENCY  DATA 
MILL  CREEK  NEAR  LOS  MOLUTOS,  CALIFORNIA 
(See  sec.  3) 
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ANALYTICAL  COMPUTATION  Of  PtaK-’TjOW  FREQUENCY  CURVE 

HILL  CREEK  NEAR  LOS  MOLINOS,  CALIFORNIA 

(See  par.  4-Oj) 

(1) 

(2) 

(3) 

(4) 

(5) 

Water 

Flow 

Log 

Dev. 

Year 

(c.f .6. ) 

(X) 

(*) 

X2 

1928-29 

1,520 

3.18 

-.49 

.240 

30 

6,000 

3-78 

.11 

.012 

31 

1,500 

3-18 

-.49 

.240 

32 

5,440 

3-74 

•07 

.005 

33 

1,080 

3*03 

-.64 

.410 

34 

2,630 

3-42 

-*25 

.062 

1934-35 

4,010 

3.60 

-.07 

.005 

36 

4,380 

3-64 

-.03 

.001 

37 

3,310 

3-52 

-.15 

.022 

38 

23,000 

4.36 

.69 

.476 

39 

1,260 

3.10 

-•57 

.325 

1939-40 

11,400 

4.06 

■39 

.152 

41 

12,200 

4.09 

.42 

.176 

42 

11,000 

4.04 

•37 

.137 

43 

6,970 

3-84 

.17 

.029 

44 

3,220 

3.51 

-.16 

.026 

1944-45 

3,230 

3*51 

-.16 

.026 

46 

6,180 

3*79 

.12 

.014 

47 

4,070 

3.61 

-.06 

.004 

48 

7,320 

3-86 

•19 

.036 

49 

3,870 

3.59 

-.08 

.006 

1949-50 

4,430 

3.65 

-.02 

.000 

51 

3,870 

3-59 

-.08 

.006 

52 

5,280 

3-72 

■05 

.002 

53 

7,710 

3*89 

.22 

.048 

54 

4,910 

3-69 

.02 

.000 

1954-55 

2,480 

3-39 

-.28 

.078 

56 

9,180 

3-96 

• 29 

.084 

57 

6,140 

3-79 

.12 

.014 

58 

6,880 

3.84 

• 17 

.029 

(6)  N 

30 

(7)  EX 

109-97 

-.13 

2.665 

(8)  M 

3.666 

(9)  2X2 

405.7813 

s2  - 2.6679/29  ■=  .091996 

(10)  (Exr/N 

403.U.34 

(11)  £x2 

2.6679 

s - .303 

(12) 

(13) 

(14) 

(15) 

(16)  (17)  (18)  (19) 

(20)  Pn 

0.25 

1.0 

10 

50  90  99 

99.75 

(21)  k (Ex|38) 

3-09 

2.50 

1.33 

0 -1.33  -2.50  -3.09 

(22)  Log  Q (Eq 

5)  4.602 

4.424 

4.069 

3.666  3.263  2.908  2.730 

(23)  Q,  cfs 

40,000 

26,600 

11,700 

4,630  1,830  609 

537  j 

Note:  Columns  4 and  5 are  not  required  when  desk  calculator  is  available, 
but  are  shown  to  illustrate  procedure  usable  without  desk  calculator.  The 
difference  between  2.665  and  2.668  is  due  to  rounding  values  in  column  4. 


EXHIBIT  4 


ANALYTICAL  FREQUENCY  COMPUTATION  AND  ADJUSTMENT 


(See  sec.  S) 


I Computation  of  Statletlce 


II  Ajjuetwnt  of  Statletlce 


MlX  flow(Q) 

! 

( 

Log  4 

su. 

1 St*.  2 

Sta.  1 

SU.  2 

Sta.  2 

cf« 

cf» 

(X) 

(T) 

(2) 

(3) 

(1) 

(5) 

(6) 

1911- 12 

1912- U 

1913- 1“* 
1911-15 

1915- 16 

1916- 17 

1917- 10 
1910-19 

1919- 20 

1920- 21 

1921- 22 

1922- 23 

1923- 21 
1921-25 

1925- 26 

1926- 27 

1927- 20 
1920-29 

1929- 30 

1930- 31 

1931- 32 
1932  33 
1933-31 
1931-35 

1935- 36 

1936- 37 

1937- 30 
1930-39 
I939-IO 

1910- 11 

1911- 12 

1912- 13 

I9U3.UU 

1911-15 

1915- 16 

1916- 17 

1917- 18 
1916-19 
1919-50 

1950- 51 

1951- 52 

1952- 53 

1953- 51 
1951-55 
193,-56 

1956- 57 

1957- 58 


1,520 

5,500 

3-18 

3.71 

6,000 

25.500 

3-78 

U.ll 

1,500 

5,570 

3.18 

3.75 

5,110 

9,900 

3-71 

1.00 

1,000 

5,100 

3.03 

3.71 

2,630 

11,100 

3.I2 

1.05 

1,010 

25,500 

3.60 

1.11 

1,300 

30,200 

3-61 

1.50 

3.310 

7,920 

3-52 

3.90 

23.000 

93,000 

1.36 

1.97 

1,260 

3,230 

3.10 

3.51 

11,100 

60,200 

1.06 

1.70 

12,200 

30,300 

1.09 

1.18 

11,000 

35.100 

1.01 

1-55 

6.970 

51,300 

3.61 

1.73 

3,220 

8,160 

3-51 

3-93 

3,230 

20,600 

3-51 

1.46 

6,100 

22,000 

5.79 

I.3I 

1,070 

17,000 

3.61 

1.25 

7,320 

16,600 

3.86 

1.22 

3.B70 

6,110 

3.59 

3.79 

1,130 

17,900 

3-65 

1.25 

3,070 

50,200 

3 59 

1.70 

5,200 

21,000 

3-72 

1,32 

7,710 

10,000 

3-09 

1.60 

1.910 

22,900 

3.69 

1.36 

2,100 

5,900 

3 39 

3.77 

9,100 

101,000 

3.96 

5.02 

6,110 

32.700 

3.79 

1.51 

6,680 

39,300 

3.01 

4.59 

EX 

109.97 

120.66 

200.61 

M 

3.666 

1.209 

. 1.269 

(7) 

(0) 

(9) 

EX2 

105.7013 

556.5276 

062.3761 

(EX)2./* 

W.U31 

55) .9511 

855.5191 

Ex2 

2.6679 

1.5762 

>.8570 

E^/N  - 1 

.0920 

.1570 

.1273 

S 

•303 

■ 397 

■357 

171.5925 

171.6900 


( 7) 

S1  * 

.303 

(Eq.  3) 

( 8) 

S2 

.397 

(Eq.  3) 

( 9) 

02  * 

.357 

(Eq.  3) 

(10) 

*1  * 

3.666 

(Eq.  2) 

(11) 

*2  * 

1.209 

(Eq.  2) 

(12) 

"4  • 

1.269 

(Eq.  2) 

(13) 

H2  - 

( 2.0915  )2 

_ - -665  (Eq.  6) 

( 2.6679  > ( 1-5762  ) 


(11)  ii2  . 1 - ( 1 - .685  ) 2b  ’ ,67  (E'1'  7) 

(15)  - .303  ♦ (.357“. 397)  (.605)  (-303/-3 97) 

■ .202  (Eq.  9) 

(16)  - 3.666  ♦ (1.269  - 1.209)  ( .685 )^( .303/ .397) 

- 3.653  (&»•  10) 

(17)  g - .00  (Paragraph  l-03c) 


(18)  - - 

1 1 

_J9 . 

-§(•*> 

39.6  yMr* 

(Eq.  u) 

111  Computation  of 
t*»ed  on  *2 » 

frequency 
Sx  «nd  £ 

curve 

k 

(SXa  38) 

Log  Q, 

(tV^) 

fi 

(19) 

(20) 

(21) 

(22) 

• 25 

3.01 

1.501 

31,700 

1.0 

2.15 

1.3U 

22,100 

5.0 

1/fO 

1.132 

13,600 

10.0 

1.32 

1.025 

10,600 

30.0 

0.51 

3.005 

6,380 

50.0 

0.00 

3.653 

1,500 

70.0 

-0.51 

3.501 

3,170 

90.0 

-1.36 

3.281 

1,910 

95*0 

-1.70 

3.171 

1,190 

99-0 

-2.15 

2.962 

Sa6 

99.75 

-3.01 

2.801 

637 

EXHIBIT  5 


ERRORS  OF  ESTIMATED  VALUES 


As  Coefficients  of 

(See  par 

Standard  Deviation 

'.  10-03) 

Level  of 

CM  nni  fM  - 

Years 

of 

Exceedence  Frequency  in  Percent 

cance# 

(N) 

0.1 

1 

10 

50 

90 

99 

99-9 

.05 

5 

4.41 

3-4i 

2.12 

.95 

.76 

1.00 

1 22 

10 

2.11 

1.65 

1.07 

.58 

• 57 

.76 

.94 

15 

1.52 

1.19 

•79 

.46 

.48 

• 65 

.80 

20 

1.23 

.97 

.64 

• 39 

.42 

• 58 

.71 

30 

•93 

.74 

• 50 

• 31 

• 35 

.49 

.60 

4o 

• 77 

.61 

.42 

.27 

.31 

• 43 

.53 

50 

.67 

• 54 

.36 

.24 

.28 

•39 

• 49 

70 

• 55 

.44 

• 30 

.20 

.24 

• 3* 

.42 

100 

.45 

.36 

.25 

.17 

.21 

• 29 

• 37 

• 25 

5 

1.4l 

1.09 

.68 

• 33 

• 31 

.41 

.49 

10 

.77 

.60 

• 39 

.22 

.24 

• 32 

•39 

15 

• 57 

• 45 

.29 

.18 

.20 

.27 

•34 

20 

•47 

•37 

.25 

.15 

.18 

.24 

.30 

30 

.36 

.29 

• 19 

.12 

• 15 

.20 

• 25 

40 

‘30 

.24 

.16 

.11 

• 13 

.18 

.22 

50 

.27 

.21 

.14 

.10 

.12 

.16 

.20 

70 

.22 

•17 

.12 

.08 

.10 

.14 

.18 

100 

.18 

.14 

.10 

• 07 

•09 

.12 

• 15 

• 75 

5 

- .49 

- .41 

- .31 

-.33 

- .68 

-I.09 

-l.4l 

10 

- .39 

- -32 

- .24 

-.22 

- .39 

- .60 

- .77 

15 

- .34 

- .27 

- .20 

-.18 

- -29 

- .45 

- .57 

20 

- .30 

- .24 

- .18 

-.15 

- .25 

- -37 

- .47 

30 

- -25 

- .20 

- .15 

-.12 

- -19 

- .29 

- .36 

40 

- .22 

- .18 

- .13 

-.11 

- .16 

- .24 

- .30 

50 

- .20 

- .16 

- .12 

-.10 

- .14 

- .21 

- .27 

70 

- .18 

- .14 

- .10 

-.08 

- .12 

- .17 

- .22 

1 

100 

- .15 

- .12 

- .09 

-.07 

- .10 

- .14 

- .18 

• 95 

5 

-1.22 

-1.00 

- .76 

-.95 

-2.12 

-3.41 

-4.41 

10 

- .94 

- .76 

- -57 

-.58 

-1.07 

-1.65 

-2.11 

15 

- .80 

- .65 

- .48 

-.46 

- <79 

-1.19 

-1.52 

20 

- .71 

- .58 

- .42 

-39 

- .64 

- .97 

-1.23 

30 

- .60 

- .49 

- .35 

-.31 

- .50 

- .74 

- .93 

4o 

- -53 

- .43 

- .31 

-.27 

- .42 

- .61 

- .77 

50 

- .49 

- -39 

- .28 

-.24 

- .36 

- .5^ 

- .67 

70 

- .42 

- .3* 

- .24 

-.20 

- .30 

- .44 

- .55 

100 

- -37 

- .29 

- .21 

-.17 

- .25 

- .36 

- >5 

* Chance  of  true 
given  error. 

value  being  greater 

than  sum  of  : 

normal -curve  value 

and 

EXHIBIT  6 


mrr-sffivm  eosswasioii 


(See  per.  10-03) 

(Based  on  equivalent  of  39-year  record) 

(1) 

(2) 

(3) 

<i) 

(5) 

(6) 

(7) 

(8) 

(9) 

P„ 

0. 1 

1.0 

10 

50 

90 

99 

99-9 

(I?) 

k (Ex  39) 

3.09 

2.33 

1.28 

0 

-1.28 

-2.33 

-3.09 

(11) 

Log  Q (Eq  5) 

L.52I 

1,310 

L.OlL 

3.653 

3.293 

2.996 

2.781 

(12) 

Q.  cfe 

33, 1 LOO 

20, LOO 

10,300 

L.5O0 

1,960 

991 

60, 

(13) 

P„,  (Ex  Lo) 

0,20 

1.31 

10.6 

50 

89. L 

98.67 

99.80 

(Plot  Q vs.  PH) 

uv) 

.05  error  la  8 unite  (ESt  6) 

.79 

.62 

.L3 

.27 

.31 

.LI* 

.9* 

(15) 

.05  error,  log 

.223 

.175 

.121 

.076 

.008 

.121 

.152 

.05  Limit-curve  value,  log 

1,71*1 

L.L85 

L.135 

3.729 

3.381 

3.120 

2.933 

(17) 

i05  Limit-curve  value,  cfs 

55,800 

30,500 

13,600 

5,360 

2, LOO 

1,320 

858 

(Plot  ve.  PJ 

(10) 

■95  error  In  S units  (Ex  6) 

-.51* 

-.LL 

-.31 

-.27 

-,L3 

-.62 

-.79 

(19) 

•95  error,  log 

-.152 

-.12L 

-.038 

-.076 

— a 12  L 

-.175 

-.223 

(20) 

.95  limit-curve  value,  log 

1.372 

L.186 

3.926 

3.577 

3.172 

2.861 

2.558 

(21) 

•95  limit-curve  value,  cfo 

23,500 

15,300 

8.L80 

3,780 

1,190 

662 

361 

(Plot  vs.  Pj 

NOTE: 

See  Exhibit  5 for  supporting  computations 

HIIIIIIIIE 

IHUIlii 


BUI 

miii 

■inn 


.05  Limif  curve 


Frequency  curve 
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ILLUSTRATIVE  EXAMPLE 
DURATION  FREQUENCY  CURVE 


iimiinpnpmiiii  impiiipp 


ILLUSTRATIVE  EXAMPLE 

HIAUTICAL  PRMUWCY  CCMPUTAllO* 

USUtO  PRE-RECORD  DATA 

Locationi  WiUaaaUe  R at  Albany,  Oregon 

Period  of  record!  13y 3-1958  (65  years) 

Ysart  of  hiatorioal  flood  e»U4»tesi  1861,  1687.  1890  (5  largest  known) 
Period  corered  by  history  and  recordi  1856“l95fl)(100  years) 

(See  par.  4-08) 

Even 
■ *0. 
(1) 

k 

W-100 

TFp 

. 6.00 
r.'tt 

TS1 

Smaller 

~wr 

log 

0 

(5) 

Event 

No.  *.100 

Wi TFT — 

K * 4.0C 
5^38 

fir 

Sa&ll.r 

51 

log 

“Tiol 

1 

u*r'8 

6.58 

2.5) 

40 

3.84 

J.04 

1.98 

2 

6.17 

6.17 

2.46 

41 

3.81 

3.91 

1.97 

3 

5 * 

2*4? 

42 

3.77 

3.77 

1.97 

4 

>.ni 

6.4  5 

5."1 

f * 

43 

3.7) 

3.7) 

1.96 

5 

5*70 

6*00 

5.'*^ 

?.5* 

44 

3.69 

5.69 

1.96 

6 

5,60 

5.77 

5.60 

2.36 

45 

3.65 

5.65 

1.96 

7 

5.51 

5.61 

5*51 

2.3) 

46 

5.60 

3.60 

1.94 

0 

5.44 

5.48 

5.44 

2.52 

47 

5.56 

3.56 

1.94 

9 

5.57 

5.58 

5.57 

2.31 

40 

3.52 

1.91 

10 

6.51 

5.?9 

6.29 

2.51 

49 

3.47 

3.47 

1.91 

u 

5.25 

5.20 

5.70 

2.29 

50 

3.43 

3.4) 

1.90 

12 

5.20 

5.15 

5.13 

2.29 

51 

3.38 

3.58 

1.89 

15 

5.15 

5.05 

5.05 

2.20 

52 

3.34 

3.54 

1.88 

14 

5.10 

4.99 

4.99 

2.26 

53 

5.29 

3.29 

i.ea 

15 

6.06 

4.95 

4.9) 

2.25 

54 

3.24 

5.24 

i«0d 

16 

4.9? 

4.87 

2.25 

55 

5.18 

3.18 

1,08 

17 

4.02 

4.82 

2.23 

56 

3.1) 

3.13 

i,e6 

ie 

4.76 

4.76 

2.22 

57 

3.07 

3.07 

1.85 

19 

4.71 

4.71 

2.22 

56 

5.01 

5.01 

1.85 

20 

4.66 

4.66 

2.18 

59 

2-95 

2.95 

1.84 

21 

4.62 

4.62 

2.16 

60 

2.07 

2.07 

1.70 

22 

4.57 

4.57 

2.14 

61 

2.80 

2.00 

1.77 

25 

6.53 

4.53 

2.14 

62 

2.71 

2.71 

1.76 

24 

4.48 

4.40 

2.14 

63 

2.62 

2.62 

1.73 

25 

4.44 

4.44 

2.13 

64 

2.52 

2.52 

1.72 

26 

4.40 

4.40 

2.15 

65 

".39 

2.39 

1.72 

27 

6.35 

6.35 

2.11 

66 

2.2) 

2.2) 

1.69 

29 

4.31 

4.31 

2.11 

67 

2.00 

2.00 

1.67 

29 

4.27 

4.27 

2.11 

66 

1.57 

1.57 

1.61 

50 

4.2) 

4.23 

2.10 

68 

68 

51 

4.1V 

4.19 

2.10 

K" 

52 

4.16 

4.16 

2.10 

EX 

277.47 

139.73 

35 

4.12 

4.12 

2.09 

M* 

4.000 

2.055 

54 

35 

4.08 

4.04 

4.08 

4.04 

2.06 

2.07 

EX2 

1206.2039 

290.1517 

36 

4.00 

4.00 

2,06 

(EX)2/*" 

1132.2000 

207.1246 

37 

5.96 

5.96 

2.06 

diff. 

74.0039 

3.0271 

38 

39 

3.92 

).an 

5.92 

5.80 

2.04 

1.99 

S = b . (3.0271/74.0039)^  - .202  (Eq.  36) 

M = a = 2.055 

- .202(4. 

000-4.000)  - 

2.039  (Eq.  21) 

C copulation  of 

Frequency  Curve  (N  * 65) 

(11) 

(12) 

(13) 

(14) 

(15)  (16) 

U7) 

(18) 

P« 

0.25 

r.o 

30 

50  90 

99 

99.75 

k (Ex. 

36) 

2.94 

2.41 

1.31 

0.00  -1.31 

-2.4l 

-2.94 

Log  Q 

if-'i'  5) 

2.633 

2.526 

2.304 

2.039  1.774 

1.552 

1.445 

Q,  thous  cfs 

630 

336 

201 

109  59.4 

35-6 

27-9 

ROTES: 

Events  numbered  4 to  9 could  not  have  occupied  higher  positions  If  the 
entire  100-yeor  record  were  available,  (but  might  occupy  lover  position*). 

Accordingly,  k values  were  selected  as  shown  above. 

4.00  was  added  to  all  k values  ^n  order  that  all  numbers  are  positive 
and  have  3 digits  for  simplicity  of  ^achinif  operation. 

iflii 


ILLUSTRATIVE  EXAMPLE 


(See  par.  4-07) 


ANALYTICAL  FREQUENCY  COMPUTATION  OMITTING  LOWER  FLOWS 


Chronological  order 

Order  of  magnitude  | 

Water 

Year 

Date 

Peak 

flow 

(cfa) 

No. 

Plotting 

position 

<*) 

Peak 

flow 

(cfs) 

Log 

of 

peak 

k 

(1) 

(2) 

(3) 

<»0 

(5) 

mm 

<7) 

(8) 

1921 

16  May 

1250 

1 

1.9 

1890 

3.28 

2.20 

1922 

6 May 

1380 

2 

4.7 

1780 

3.25 

1.74 

1923 

10  May 

1450 

3 

7.4 

1480 

3.17 

1.48 

1924 

4 May 

618 

4 

10.2 

1450 

3.16 

1.30 

1925 

8 May 

523 

5 

12.9 

1430 

3.16 

1.15 

1926 

21  Apr 

508 

6 

15.7 

1380 

3.14 

1.03 

1927 

30  Apr 

1220 

7 

18.4 

1300 

3.11 

0.91 

1928 

28  Apr 

1180 

8 

21.1 

1280 

3-11 

0.81 

1929 

14  May 

1060 

9 

23.9 

1250 

3.10 

0.72 

1930 

25  Apr 

412 

10 

26.6 

1220 

3.09 

0.63 

1931 

5 May 

170 

11 

29.4 

1180 

3.07 

0.55 

1932 

14  May 

1480 

12 

32.1 

1090 

3.04 

0.47 

1933 

21  May 

876 

13 

34.9 

1060 

3-02 

0.39 

1934 

21  Jul 

113 

14 

37-6 

1020 

3-01 

0.32 

1935 

10  May 

516 

15 

40.4 

1000 

3.00 

0.25 

1936 

4 May 

1780 

16 

43.1 

995 

3.00 

0.17 

1937 

8 May 

1090 

17 

45.9 

985 

2.99 

0.11 

1938 

22  Apr 

760 

18 

48.6 

876 

2.94 

0.04 

1939 

30  Apr 

397 

19 

51.3 

788 

1940 

21  Apr 

282 

20 

54.1 

788 

1941 

2 May 

353 

21 

56.9 

760 

1942 

13  Apr 

597 

22 

59.6 

678 

1943 

23  Apr 

995 

23 

62.3 

618 

1944 

14  May 

6ll 

24 

65.1 

611 

1945 

4 May 

985 

25 

67.9 

611 

1946 

18  Apr 

1430 

26 

70.6 

597 

1947 

3 May 

788 

27 

73*4 

523 

1948 

17  May 

1280 

28 

76.1 

516 

1949 

24  Apr 

1020 

29 

78.9 

508 

1950 

18  May 

1300 

30 

81.6 

412 

1951 

12  May 

1000 

31 

84.3 

409 

1952 

3 May 

1890 

32 

87.I 

397 

1953 

29  May 

611 

33 

89.8 

353 

1954 

25  Apr 

409 

34 

92.6 

282 

1955 

6 May 

788 

35 

95.3 

170 

1956 

24  Dec 

678 

36 

98.1 

113 

N" 

18 

18 

Computation  of  curve  (N=36) 


(9) 

PH 

.0025 
.01 
.1 
• 5 


(10) 

(11) 

(12) 

k 

log  Q 

Q 

3-04 

3.437 

2,740 

2.47 

3049 

2,230 

1.32 

3.172 

1,490 

0 

2.969 

931 

S = b 
M = a 


EX  55.64 

M"  3.091 

EX2  172.1336 

(ex)2/m"  171.9694 
Ex2  .1442 

(.144/6.101)^ 
3.091  - .154( .793) 


14.2? 

0.793 

17.4139 

11.3129 

6.1010 

= 0.154 
= 2.969 


EXHIBIT  15 


Peak  flow  in  c 
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StotlomugGS  ll-^RlO  1 & 
EXHIBIT  15 
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MAXIMUM-RUNOFF  VOLUME  FREQUENCY  DATA  - LOGARITHMS 


MAX t MUM- RUNOFF  VOLUME  FREQUENCY  COMPUTATIONS 
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COMPUTATION  PROCEDURE 
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EXHIBIT  19 


STANDARD  DEVIATION  ADJUSTMENT 


Average  flow  in 


Exceodonce  frequency  per  hundred  years 


l 


LHC  END 

■ Annual -event  curves  adjusted 
with  base  station. 

.Partial-duration  curve  based 
on  Langbein  criteria. 

■adopted  partial-duration  curve 
based  on  plotted  pointB* 
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NOTES i 1.  Annual -l -ent  curves 
adjusted  using  record  of  Feather 
River  at  Bidvell  Ba?, 

2.  See  Exhibit  18  for 
supporting  computations . 


Exceedence  interval  in  youi 


ILLUSTRATIVE  EXAMPLE 


MAX-RUNOFF  VOLUME  FREQUENCY  CURVES 
MILL  CREEK  NR  LOS  MOLINOS 


(See  par.  6-04) 


volume  in  thousond 


ILLUSTRATIVE  EXAMPLE 
REGIONAL  FREQUENCY  CORRELATION 


(See  seoa.  7 A 9) 


X^  ■ 1 1 log. 

, S 

X2  - log. 

D.A. 

X j =*  log. 

No.  rainy  days  per  year 

: j 

Sta.  No. 

X2 

X3 

X1 

Sta.  No. 

X2 

X3 

X1 

i 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(?) 

(8) 

1 

1.6l 

2.11 

0.29 

33 

1.94 

1.87 

0.20 

] 

2 

2.89 

2.12 

0.18 

34 

2.73 

1.36 

0.58 

■ 7 

3 

4.38 

2.11 

0.17 

35 

3.63 

1.81 

0.64 

■ 3 

4 

3.20 

2.04 

0.44 

36 

1.91 

1.58 

0.37 

1 

5 

3.92 

2.07 

O.38 

37 

2.26 

1.48 

0.27 

, | 

6 

1.6l 

2.04 

0.37 

38 

2.97 

1.89 

0.54 

7 

3.21 

2.09 

0.30 

39 

0.70 

1.32 

O.63 

8 

3.65 

1.99 

0.35 

40 

0.30 

1.54 

0.78 

i 

9 

3.23 

2.15 

0.16 

4l 

3-38 

1.62 

0.46 

] 

10 

4-33 

2.08 

0.11 

42 

2.87 

2.03 

0.44 

11 

1.60 

2.09 

0.32 

43 

2.42 

2.26 

0.24 

12 

2.82 

2.00 

0.34 

44 

4.53 

1-93 

-0.03 

13 

2.40 

2.00 

0.25 

45 

3-04 

1.78 

0.30 

14 

3.69 

2.09 

0.43 

46 

4.13 

2.00 

0.17 

15 

2.18 

2.19 

0.27 

47 

1.49 

2.01 

0.l4 

16 

2.09 

2.17 

0.25 

48 

5.37 

1.95 

0.10 

5 

17 

4.48 

1.91 

0.52 

49 

1.36 

2.11 

0.27 

. 4 

18 

4.95 

1.95 

0.18 

50 

2.31 

2.23 

0.18 

i 

19 

2.21 

1.97 

0.39 

1 

20 

3.41 

2.08 

0.40 

EX 

147.55 

96.24 

17.89 

, ‘ 2 

21 

4.82 

1.88 

0.25 

M 

2.951 

1.925 

0.358 

5 ? 

22 

1.76 

1.93 

0.23 

exx2 

503.777 9 285-5627 

51.1527 

v i 

23 

4.39 

1.74 

0.54 

EXEXg/N 

435.4200  284.0042 

52.7934 

| 

24 

3.23 

2.01 

0.51 

Exxg 

68.3579 

1.5585 

-1.6407 

' ] 

25 

R.S8 

2 04 

n Le; 

: j 

26 

1.64 

1.78 

0.63 

2XX3 

187.5912 

33.2598 

27 

4.58 

1.76 

0.45 

EXEXo/N 

185.2428 

34.4347 

) 

28 

3.26 

1.93 

0.59 

Txxi 

1.5585 

2.3484 

-1.1749 

29 

4.29 

1.81 

0.46 

30 

1.23 

1.89 

0.32 

8.1635 

i 

31 

3.44 

1.48 

O.96 

EXEXi/N 

6.4010 

32 

2.11 

1.97 

0.12 

Exxi 

x 

-1.6407 

-1.1749 

1.7625 

' = 

68.4  bg  + 

1.56  b3 

* -1.64 

(Eq.  23) 

l b2  = 

-.013 

1.56  b2  + 

2.35  b3 

= -1.17 

(Eq.  24) 

/ ‘3- 

-.49 

5 

a = O.36  + .013  (2.95)  + .49  (1.92)  = 1.34  (Eq.  28) 


= 1.34  - .OI3X2  - *49X2  (Regression  Equation,  eq.  22) 

R2  = Ijl01-3  .(ij-,,-.64 49 ■_!? ) „ >333  (Eqj  2gj 

*2  = 1 - (.662)  49/47  - .310  (Eq.  31)  R = .56 


EXHIBIT  22 


Standard  error  Standard  error 
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EXHIBIT  24 
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CIVIL  WORKS  INVESTIGATIONS  PROJECT  NO  151 

FLOOD  VOLUME  STUDIES-* 

WEST  COAST 

Sp 

FREQUENCY  INDEX  FOR 

PEAK  RUNOFF 

(See  par  7-1 

i) 

io  o io  2o  so  eo 
scale  tHHSffi=rrrs=i ts==~t  miles 

U.S  Army  Enjintor  Dietiict 

Socromento,  Californio 

L.H.  8.  1 Sept.  I960 

Filo  No.  SA-OT -26-68 

EXHIBIT  25 


;i.  ,lkuLL.  laul/riiWliiiE 


Central  pressure  index  in  inches  Central  pressure  index  in  inches. 


Occurrences  per  hundred  years 


Occurrences  per  hundred  years. 


O.l  0.5  I 2 5 10  20  40  60  BO  90 

Occurrences  per  hundred  years. 


NOTE: 

Drawings  prepared  for  the 
Corps  of  Engineers  by  the  U S. 
Weather  Bureau. 


ILLUSTRATIVE  EXAMPLE 

ACCUMULATED  FREQUENCY  OF  HURRICANE  CENTRAL  PRESSURES 
(PLOTTED  AS  FREQUENCY  PER  100  YEARS  BASED  ON  1900-1967) 


■ 


3 

| 


(See  par.  8-04) 


EXHIBIT  28 


Annual  sediment 


vx> 


•fitl.Br  »f*#t  l*»*l 


sett 


240 


*4o 


220, 


200 


EQTPt 

1.  Cirolad  point*  *rt  ki»4  a*  routiner*  ot 
indie*  Ud  hypothatia*l  flood*  with  indioalad 
initial  *ta#*a.  Ooaaareatloft  ralaaaaa  on.'ijr 
krt  mu  It  Itlfti  fcalov  200,000  fto-ft. 

2.  ^raqaasoT  ourr*  would  ha  W«*d  «o  aulti- 
purpoaa  routing*  ot  lowar  ilafti  and  on  4hU 
ooaputatiok  ot  tha  hijfear  atof»>  ahovn  Mreon. 


lottfla  of  f lood-qontrol  apao* 


Com  pat*  tic*  ot  ooopoaito  ourrei 


Max 

tttdaadano*  freauwnoy 

Stora** 

A 

8 

0 

D 

B* 

*» 

0.7 

0.1 

0.0 

0.0 

.16 

500 

2.0 

0.7 

0.2 

0.0 

.62 

280 

4.1 

2.0 

0.9 

0.5 

1.64 

260 

m 

4.2 

2.5 

0.9 

5.5 

240 

10.0 

7.0 

4.4 

2.1 

5.6 

220 

15.2 

10.0 

7.1 

5.« 

6.2 

200 

16.1 

15.5 

10.1 

5.4 

11.0 

•WiflhUd  avara## 


99.99 


1 


70  “SO  5o  40  50 


94.9  99.5  99  9« 95  90 


.5  .2  .1  .05 


Exoaadatioa  froquar.ct  par  hundred  mn 


Poroant  of  <•  4**  or*  axoaadad 


A - Initial  atonuf*,  200,000  ao-ft.  (20jC  of  tin* ) 

fi  - Initial  ■toraif*,  180,000  ao~ft.  (20jt  of  tin*) 

C - Initial  •torwuy* , 155,000  ao-ft.  (Hot  of  ti»a) 

0 - Initial  atora**,  106,000  ao-ft,  (50^  of  tlaa) 


ILLUSTRATIVE  EXAMPLE 

STORAGE  FREQUENCY  COMPUTATION  BASED  ON  COINCIDENT  FREQUENCIES 

(See  par.  8-06) 


J 


EXHIBIT  30 


ILLUSTRATIVE  example 


COMPUTATION  OP  SIMPLE  LItfEAR  CORHEUTICW 


(Sea  par.  9-04) 


Max . 12  hr  Mean  Annual 


7-0-12 

3.1 

.49 

10.0 

1.00 

7-0-309 

4.6 

.66 

19.2 

1.26 

7-0-15 

7.5 

.88 

30.8 

1.49 

7-P-20 

5.2 

.72 

13.3 

1.12 

7-0-22 

12.6 

1.10 

42.2 

1.63 

7-P-25 

3.1 

.49 

12.2 

1.09 

7-0-23 

10.4 

1.02 

36.2 

1.56 

7-P-61 

4.9 

.69 

12.2 

1.09 

7-0-36 

4.9 

.69 

15.4 

1.19 

8-0-8 

5-2 

.72 

25.2 

1.40 

7-0-39 

7.1 

.85 

20.0 

1.30 

8-0-18 

4.6 

.66 

18.5 

1.27 

7-0-43 

5-3 

.72 

14.2 

i.15 

8-0-29 

11.8 

1.07 

32.3 

1.51 

7-0-77 

9.7 

.99 

24.5 

1.39 

8-0-34 

12.4 

1.09 

33.8 

1.53 

7-0-84 

6.2 

.79 

19.3 

1.29 

8-0-35 

10.1 

1.00 

29.8 

1.47 

7-0-89 

10.0 

1.00 

18.8 

1.27 

8-0-45 

6.0 

.78 

19.7 

1.29 

7-0-93 

6.0 

.78 

16.2 

1.21 

8-0-60 

8.2 

.91 

22.5 

1.35 

7-0-95 

5.0 

.76 

18.2 

1.26 

8-0-67 

9.8 

.99 

28.7 

1.46 

7-0-99 

4.5 

.65 

15.8 

1.20 

8-0-75 

11.7 

1.07 

36.9 

1.57 

7-0-102 

5.3 

.72 

17.3 

1.24 

8-0-219 

10.0 

1.00 

27.8 

1.44 

7-0-110 

9.3 

.97 

21.3 

1.33 

7.0-434 

6.7 

.83 

23.1 

1.36 

7-0-114 

8.5 

.93 

26.8 

1.43 

7-0-120 

11.6 

1.06 

29.0 

1.46 

N 

42 

42 

7-0-122 

10.5 

1.02 

28.8 

1.46 

ac 

36.68 

56.81 

7-0-124 

11.6 

1.06 

31.1 

1.49 

M 

•873 

1.353 

7-0-125 

9-9 

1.00 

27.9 

1.45 

sxf 

7-0-127 

12.9 

1.11 

32.2 

1.51 

33.2134 

77.80J 

7-0-130 

13.4 

1.13 

44.6 

1.65 

(«);/& 

32.0339 

76.845 

7-0-133 

8.4 

.92 

23.8 

1.38 

Ex2 

1.1795 

.96] 

7-0-136 

7.9 

•90 

27.2 

1.43 

7-0-149 

8.0 

•90 

19.2 

1.28 

E(X,~  ) 

50 

.5601 

7-0-182 

5.2 

•72 

16.1 

1.21 

i ErExf/N 

49.6141 

7-0-190 

6.9 

.84 

20,9 

1-32 

E<V2) 

.9**60 

1 

b2  " 

.9460/. 9618 

a 0.< 

?8  (Eq.  20) 

a a 

•873 

- .98  l 

U.353) 

- -.453 

(Eq. 

21) 

log  ^ 


0.?86 


. 1.1795 

1 - (0.2l4)  4l/4o  = 0.78  (Eq.  Bl)  R 
- 0.45  (regression  eq.,  Eq.  22) 


.98  log  Xg  - 0.45  - log  (.36  X^’?8) 

• q8 

•36  Xg'^  (transformed  regression  Eq.) 
0.22  (l.l8)/4l  a .0063  (Bq.  35) 

.079  antilog  2 S a 1.44 


EXHIBIT  31 


30 


ILLUSTRATIVE  EXAMPLE 

COMPUTATION  OF  MULTIPLE  LINEAR  CORRELATION 


(See  par.  9-0 6) 


Water 

Log. 

Log. 

Log. 

Log. 

yea" 

Snow  Cover 

Ground  Water 

April  l'reoip. 

Runoff 

(1) 

(2) 

(3) 

(4) 

(5) 

1936 

.399 

.325 

.710 

• 939 

1937 

•3^3 

.385 

.634 

.945 

1938 

• 369 

.408 

.886 

1.052 

1939 

.246 

.428 

.581 

.744 

1940 

.181 

.316 

1.027 

.666 

1941 

.297 

.460 

1.315 

1.081 

1942 

.299 

.511 

1.097 

1.060 

1943 

.354 

.379 

•707 

.892 

1944 

.295 

.395 

1.240 

1.021 

1945 

.321 

• 376 

1.091 

.920 

1946 

.168 

.413 

I.038 

.755 

1947 

.280 

.410 

.979 

.960 

EX 

3.552 

4.806 

11.305 

11.035 

M 

.296 

.400 

.942 

.920 

DCEXf/ 

£(xxr,) 


e(xx3) 

£XEX,/n 

E(xxp 

S(XX4) 

EXEXl*/N 

£(xx4) 

£(xxL) 

EXEXj/N 

£(x*l) 


1.4197 

3.2898 

3.3365 

1.4226 

3-3463 

3.2664 

-.0029 

-.0565 

.0701 

-.0029 


-.0565 


1.9556 

1.9246 

.0310 


.0453 


4.5730 

4.5277 

.0453 


11.2796 

10.6502 

.6294 


10.5224 

10.3959 

.1265 


10.3468 

10.1476 

.1992 


.0545  b2  - .0029  b3  - .0565  b4  = .0701  (Eq.  25)  ] b2  = 1,623 
-.0029  bg  + .0310  b3  + .0453  b4  « .0392  (Eq.  26)  > b3  * 1.012 
-.0565  b2  + .0453  + .6294  b4  = .1265  (Eq.  27)  J bj,  = 0.274 

a = .920  - 1.623  (.296)  - 1.012  (.400)  - .274  (.942)  - -.223  (Eq.  28) 
X4  = 1.623  X2  + 1.012  Xj  + 0.274  X4  - 0.223  (Regression  equation) 


R2  s 1.623  (.0701)  + 1.012  (.0392)  4 0.274  (.1265)  „ <944  ^Eq<  ^ j 


“ .1992 

R2  a 1 - (.056)11/8  - 0.923 
s|  - .077  (.1992)/!!  = 0.00139 


R = O.96 


(Eq.  31) 
(Eq.  35) 


S*  = .037 


antilog  2 Se  = 1.18 


Exceedence  probability 


d.f , 

.99 

■ 95 

• 90 

.70 

-50 

.30 

.10 

.05 

.01 

i 

.0*157 

•00393 

.0158 

.148 

.455 

1.074 

2.706 

3.841 

6.635 

2 

.0201 

.103 

.211 

-713 

1.386 

2.408 

4.605 

5.991 

9.210 

3 

■ 115 

•352 

■ 584 

1.424 

2.366 

3.665 

6.251 

7.815 

11.345 

it 

.297 

• 711 

1.064 

2-195 

3.357 

4.878 

7.779 

9.488 

13.277 

5 

• 554 

1.145 

1.610 

3.000 

4.351 

6.064 

9.236 

11. 070 

15.086 

6 

.872 

1.635 

2.204 

3.828 

5.348 

7.231 

10.645 

12.592 

16.812 

7 

1.239 

2.167 

2.633 

4.671 

6-346 

8.383 

12.017 

14.067 

18.475 

8 

1.646 

2.733 

3.490 

5-527 

7.344 

9-524 

13.362 

15.507 

20.090 

9 

2.088 

3.325 

4.168 

6-393 

8.343 

IO.656 

14,684 

16.919 

21.666 

1° 

2.558 

3.940 

4.865 

7-267 

9.342 

11.781 

15.987 

18.307 

23.209 

11 

3.053 

4.575 

5.578 

8.148 

10.341 

12.899 

17.275 

19.675 

24.725 

12 

3-571 

5.226 

6.304 

9-034 

11.340 

14.011 

18.549 

21.026 

26.217 

13 

4.107 

5.892 

7-042 

9.926 

12.340 

15.H9 

19.812 

22.362 

27.688 

14 

4.660 

6.571 

7.790 

10-821 

13-339 

16.222 

21.064 

23.685 

29-141 

15 

5-229 

7.261 

8.547 

11.721 

14-339 

17-322 

22.307 

24.996 

30.578 

16 

5-812 

7-962 

9.312 

12.624 

15.338 

18.418 

23.542 

26.296 

32.000 

17 

6.408 

8.672 

10.085 

13.531 

16.338 

19.5U 

24.769 

27-587 

33-409 

18 

7-015 

9.390 

10.865 

14.440 

17.338 

20.601 

25.969 

28.869 

34.805 

.19 

7-633 

10.117 

U.651 

15-352 

18.338 

21.689 

27.204 

30.144 

36.191 

20 

8.260 

10.851 

12.443 

16.266 

19*337 

22.775 

26.412 

31.410 

37-566 

21 

8.897 

U.591 

13-240 

17.182 

20.337 

23.858 

29.615 

32.671 

38-932 

22 

9-542 

12.338 

l4,o4l 

18.101 

21.337 

24.939 

30.613 

33-924 

40.209 

23 

10.196 

13.091 

14.848 

19.021 

22.337 

26.018 

32.007 

35-172 

41.638 

24 

10.856 

13.848 

15-659 

19.943 

23-337 

27.096 

33.196 

36A15 

42.900 

25 

11.524 

14.611 

16.473 

20.867 

24.337 

28.172 

34.382 

37.652 

44.314 

26 

12.118 

15-379 

17-292 

21.792 

25.336 

29.246 

35.563 

38.805 

45-642 

27 

12.879 

16.151 

18.114 

22.719 

26.336 

30.319 

36.741 

40.113 

46.963 

28 

13.565 

16.928 

18.939 

23.647 

27-336 

31.391 

37.916 

4l . 337 

48.278 

29 

14.256 

17.700 

19.768 

24-577 

20.336 

32.461 

39.087 

42.557 

49.588 

30 

14.953 

18.493 

20.599 

25-508 

29.336 

33.530 

40.256 

43.773 

50.892 

For  higher  degrees  of  freedom,  values  of  -yfzX^  are  distributed  approximately  normally  about  • 


with  a standard  deviation  of  1. 


EXHIBIT  35 


■“ ■am.fcjL . . „ , ' ' ; , . .J! 
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TRANSFORMED  PLOTTING  POSITIONS  (l) 

k VALUES 

HAVING  ZERO  mm  AND  UNIT  STANDARD  DEVIATION 

(See 

par. 

4-05) 

It 

11 

12 

1) 

14 

15 

14 

1? 

18 

19 

20 

21 

2? 

2) 

24 

?) 

H 

27 

26 

29 

30 

31 

32 

3) 

W 

35 

taftt 

■*.  1 

1.71 

1.75 

1.79 

1.4? 

1.95 

1.88 

1.91 

1,9) 

1.93 

1,91 

1*99 

2.01 

2.0) 

2.05 

2.07 

2.09 

2.10 

2.11 

?.» 

».» 

2.16 

2.16 

2.17 

2.14 

2*19 

? 

1.11 

1,14 

un 

1.24 

1,30 

1.33 

1,34 

1.39 

1.4? 

1.45 

1.44 

1.50 

U5* 

i.V 

1.5* 

1.5* 

1*60 

1.62 

1.64 

1.6* 

1.61 

1*68 

1.70 

un 

1*7? 

) 

o.H 

0.8? 

0.0* 

0.93 

0,54 

1.0? 

1,04 

1.10 

1.1) 

1,1* 

1*19 

un 

1,24 

1.2* 

1.28 

1.30 

Ml 

1*  )4 

l.  36 

U38 

1.40 

1.4J 

1.44 

U4« 

U47 

k 

0.44 

<.54 

0.42 

0,48 

0.73 

0.78 

0,81 

0.87 

0.91 

0,9* 

0.97 

1.00 

1,0) 

1.06 

1,09 

l.U 

1.1) 
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0.4* 

0.49 

O.I) 

0.14 

0.14 

0.18 

0.19 

0.21 

0.2) 

0.24 

0,2* 

0.27 

0,29 

0.30 

0.32 

0.)) 

0.34 

0.3* 

0»)7 

0.38 

0.39 

0.40 

0.41 

0.4) 

0.44 

0.43 

0.44 

29 

0.04 

0.12 

0.14 

0.14 

0.17 

0.19 

0.20 

0.22 

0.24 

0.25 

0.2? 

0.28 

0.29 

0.51 

0.3? 

0.3) 

o.v 

0,5* 

0»)7 

O.M 

O.W 

0.40 

0.41 

0,4) 

>0 

0,04 

0.10 

0.12 

0.1) 

0.15 

0.17 

0.18 

0.20 

0.21 

0.?) 

0.24 

0.26 

0.27 

0.29 

O.W 

0.91 

O.M 

0.)) 

0.)5 

o.H 

0*)7 

0.  38 

0*>9 

si 

0,02 

0.04 

0.04 

0.08 

0.10 

0.11 

0.1) 

0.15 

0.1* 

0.14 

0.19 

0.21 

0.2? 

0*24 

0,25 

0.26 

0.28 

0.29 

0*)0 

O.M 

0*)) 

0.34 

0.35 

0.)* 

32 

0.02 

0.04 

0.04 

0.08 

0.09 

0.11 

0,1) 

O.I1 

0.16 

0.17 

0.19 

0*20 

0.22 

»•?) 

0.24 

0.26 

0.27 

0.28 

o.y> 

0*)1 

0.3? 

0.)) 

>3 

0 

0.02 

0.04 

0.0* 

0.0? 

0,09 

o.u 

0.12 

0.14 

0.15 

0.17 

0.16 

0.20 

0.21 

0.22 

0.24 

0-35 

o.?6 

0.28 

0.29 

0,)0 

54 

0 

0.0? 

0.04 

0,0* 

0,07 

0.09 

o.n 

0.12 

0.14 

O.I) 

0.1* 

0.18 

O.I) 

0.31 

0.23 

0.2) 

0*24 

0.?6 

0*27 

0 

0,0? 

0.04 

0.0) 

0,07 

0.09 

0.10 

0.1? 

0.1) 

0.15 

0.16 

0,17 

0.19 

0.20 

0.21 

0.2) 

0.24 

Si 

0 

0.02 

0.04 

0.05 

0,07 

0,08 

0.10 

0.11 

0.1) 

0,14 

O.H 

0.17 

0.18 

0.20 

0.21 

0 

0.0? 

0.0) 

0.05 

0*07 

0.08 

0.10 

0.11 

0.1? 

0*14 

0.15 

0.17 

0,18 

0 

0.02 

0*0) 

0.05 

0.06 

0,08 

0.09 

0,11 

0.12 

0.14 

0.15 

39 

0 

0,02 

0.0) 

0.05 

0.0* 

0.08 

0.09 

0.11 

0,12 

0 

0.02 

0.0) 

0,05 

Q*« 

0.08 

0.09 

0 

0.02 

0.0) 

0.05 

0,06 

**3 

3*i M f»*3itoU«  *•  »h**t  2 

0 

0*02 

0.0) 

0 

I 

•1 

■ 

-J 

3 

I 

i 
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i 

1 

* 

i 

i 
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TRANSFORMER  PLOTTING  POSITIONS  ( l) 
k VALUES  HAVING  mo  MEAN  AND  UNIT  STANDARD  DEVIATION 

(See  pur.  A-QS) _ 

n ««  *1  *****  * X *7  X 9*  *«>  »0  170  1*0  140  1*0  1*0  l?0  1W  1*0  TOO 


?,5* 

2.6* 

2.5* 

:.v* 

2/4 

7.55 

>.55 

7.55 

7.X 

If** 

».X 

?. *7 

1.57 

7.57 

J 

1.1 1 

2.1? 

7.1? 

7,n 

7.15 

2.1* 

>.l* 

7.14 

7.15 

7.15 

7.16 

7. 1C 

1.14 

?.l] 

\ 

1.64 

1.40 

1.90 

l.fl 

!.*1 

i.w 

t.*7 

»•♦: 

!♦*) 

1»*6 

1**4 

1.95 

1.95 

l.X 

4 

1.7* 

1,75 

1.76 

1.x 

l.tl 

1.11 

1.76 

1.7* 

1.7* 

1.7* 

1.00 

1,00 

1.90 

1.91 

5 

1.4* 

\M 

1*45 

1.64 

1.46 

1.6) 

1.66 

U44 

1.67 

1.47 

1.4* 

1.6# 

1.6* 

1»*9 

6 

1.6/ 

1.55 

1.5* 

1.64 

1.45 

1.55 

1.64 

1.56 

Ml 

l.X 

1**0 

1.5* 

1.5* 

1.60 

T 

1.44 

1.44 

1.4* 

1.45 

1.44 

1.4: 

l.*7 

1,40 

1.40 

’.4* 

l*4f 

1.50 

1.50 

1.51 

» 

t.X 

W9T 

1.67 

I.X 

I.X 

..X 

1.40 

1.40 

1.41 

1.41 

l.*f 

1.41 

1.4* 

1.4* 

6 

l •?* 

1.50 

1.50 

1.61 

UK 

1.57 

l.t* 

Ml 

»,*4 

l.X 

l.*9 

l.X 

l.X 

1**1 

10 

l.*7 

1.21 

1.7* 

1./6 

1.75 

1.76 

1.76 

1.77 

1,79 

1,70 

1.7* 

1.7* 

l.X 

l.X 

11 

l.l’ 

1.1) 

1.1" 

Ml 

1.1* 

1,70 

1.71 

1.71 

1.7? 

1.10 

1.7* 

l.?4 

1.04 

1.25 

\7 

1.11 

M' 

1.1/ 

1.1* 

l.l* 

1.14 

1.16 

1.16 

1.17 

1.11 

1.10 

1.10 

1.19 

1.19 

1* 

l.Oh 

UcVf 

l.OJ 

l.l  A 

1.0* 

1.0* 

1.10 

l.ll 

Ml 

1.17 

1.17 

1.1) 

1.16 

V,  14 

14 

l,11 1 

».a? 

1.0? 

1.0) 

l.C* 

* 

1.0* 

1.06 

1.07 

1,07 

1.00 

1.09 

1.0* 

1.10 

H 

®.v> 

0.9J 

0.9" 

0.7* 

0.9* 

1,00 

1.01 

1.0) 

1,02 

1.0* 

1.0) 

1.04 

1.05 

1.0* 

16 

9.f; 

0.*.' 

0.v\ 

0.94 

0,96 

0.96 

O.X 

0.9  / 

0.90 

Ot*0 

0.*6 

1.00 

1.00 

I.OI 

IT 

©.■7 

0.96 

H.W 

0.90 

0.91 

0.91 

0.97 

0.95 

».** 

0.94 

0.95 

0.*5 

o.X 

0.97 

u 

O.f) 

0.94 

0.6'- 

D.M. 

0.66 

o.»T 

0.0* 

0.6* 

0.9* 

O.VI 

0**1 

o.w 

0.97 

0.9* 

1* 

0.7* 

0.0© 

0,91 

0.97 

0.67 

0.6* 

0.94 

0.0* 

6.05 

O.U 

0.07 

0.09 

0.9« 

0.0* 

TO 

0.T6 

0.7* 

O.t? 

O.TV 

0.  76 

O.X 

0.90 

0.41 

0,97 

0.02 

0.6) 

0.04 

0.0* 

0.95 

71 

0.  + 3 

0.77 

U.7' 

0.74 

0,76 

0.T6 

0»7t 

0.77 

*.?* 

0.7* 

0,90 

0,91 

0.01 

0.9? 

.V 

o.« 

n,« 

0/9 

M,  W 

0.71 

0.7? 

0.7) 

0.74 

0.76 

0.75 

0,76 

0.77 

0.7« 

O.TU 

7* 

o.x 

0.65 

0.64 

0,47 

0.69 

0>9 

0.70 

0,70 

0.71 

*.12 

>4 

u.**o 

0,61 

0.4: 

0,6* 

0.54 

0.6* 

0.66 

0.61 

0.50 

0.60 

?4 

0.67 

0.60 

0.4* 

o.«> 

O.kl 

0.6? 

0.6) 

0.6) 

0.64 

0.65 

?i 

0.5) 

0.55 

O.X 

(1.57 

0.67 

O.X 

0.5* 

O.60 

O.M 

0.6? 

17 

t:?, 

8:2 

0.57 

0.5' 

O.X 

0.55 

o.X 

0.57 

O.X 

0.5* 

70 

0.49 

0.50 

o.*l 

0.5: 

0.5) 

0.54 

0.55 

o.X 

2* 

0.44 

0,4', 

0.46 

0.4? 

0.40 

0.49 

O.X 

0.31 

0.5? 

0.6) 

Y> 

0.41 

0.4? 

0.4) 

0.64 

0.45 

0.46 

0.47 

0.49 

0.4* 

0,60 

*1 

0,*t 

0.59 

0.40 

0.41 

0.4? 

0.4) 

0.46 

0.45 

0.46 

0.4f 

V 

0.*H 

0.  *5 

o.x 

o.x 

0.** 

0.40 

0.4) 

0.4? 

0.4) 

0.44 

*> 

0,  >1 

O.V 

0.** 

0.45 

O.X 

0.57 

O.X 

0.59 

0.40 

0,41 

X 

O.T0 

0.2  i 

rt.tO 

O.V 

0.4t 

O.X 

O.X 

O.X 

0.  >7 

O.V0 

74 

0.?6 

9.26 

u.?1 

0,29 

o.w 

O.M 

0.  6? 

0,1) 

O.X 

0, 13 

X 

0..*. 

o.?t 

0.24 

0.76 

0.21 

0.19 

1.?* 

0.  Vj 

0,  >1 

0.  *2 

S' 

A. 10 

o.?o 

0.77 

0.7) 

0.24 

0.76 

o.:f 

0.77 

0.29 

0.  W 

X 

0.16 

0,14 

0.1* 

0.70 

0.71 

0.77 

0,<) 

o.» 

C.26 

0.77 

>9 

0.1  J 

0.1‘ 

0.16 

0.1? 

0.14 

0.17 

0.21 

0.?? 

O.P) 

0.?4 

40 

0.10 

0.1.' 

0.1) 

0.14 

0.1* 

0.17 

0.10 

0.19 

C.?n 

o.ri 

VI 

0.01 

0.CT) 

0.10 

0.1s 

0.1; 

0.1* 

0.15 

O.t* 

0.16 

0.1* 

47 

0.04 

0.06 

0.0? 

0.09 

0.10 

O.U 

0.1? 

0.14 

0.15 

0.16 

4* 

0.02 

0.0) 

o.c* 

O.ffc 

0.07 

0.09 

0.10 

0.11 

0. 12 

0.1) 

4J. 

0 

*1.0 

0.0) 

0.04 

0.06 

0.07 

0.09 

0.10 

0.11 

4A 

0 

□ .0? 

0.0) 

0.04 

0.0. 

o.ol 

O.W 

66 

0 

0.0? 

0.0) 

0.04 

0.05 

0.70  0.71 

0.6?  0 M 


0.7*  0,76 
0.71  0,7 7 


0.5?  0/4 
o.*  0.55 

o.*i  o.** 
O.V*  0.4* 
0.46  0.46 
0.4*  0.44 
0.4©  0.41 
o. *7  0.5* 
O.X  o.J* 
0.57  0.5) 
0.7*  O.V) 
0.26  0.27 
0.?4  0.2* 
0.71  0.72 
0.10  0.1* 
0.1*  0,17 
0.1*  0.14 
0.10  0.17 
0,00  0.0* 
0.0*  0.06 
0.0*  0.04 
0 0.01 


0.*>*  0.60 
O.X  0.*? 
0.**  0.54 
0,50  0.51 
0.47  O.U 
0.4'  0.46 

0.42  0.4) 
D.J*  0.40 
0.56  ).»' 
0.54  0. 
0.5)  v <• 

0.2u 

0.J6 

0.2*  ©.  1 
0.20  . \ 
0.16  '.i 

0.1*  0 
0.1*  0.1* 
0.10  O.IC 

0.09  0.09 

0,05  0,06 
O.OJ  c.04 
0 0.01 


<■*  Llwv  lliltui  or.  Miiwl  probMbllltv  crt*. 

k ?•'-  low«r  ImIT  of  *•*■1*.  »r»  *nd  In  ***n!t\*l*  tc 

tfcn*.  **  k,  f (•«»■  n'*»*.*rv!  n r#*#r»n  irli-r  mrtir*  with  ■. 

' r'-t  . 


7,61  7,64 
1.71  7.74 

7.00  2.04 

no*  i.** 

».X  1,73 
1.64  1.6* 
1.56  1.61 
l.M  1,5* 
1.4?  1,47 
1 * *6  1 »4l 
1.51  1.56 
1.76  1.51 
l.?l  1.76 
1.16  1,21 
1.17  l.)7 
1.04  l.l) 
1.04  1.0* 

1.00  1,05 
0.*6  1.0? 
0.5)  0.59 
0.6*  0.*5 
0.66  0.*? 
0.65  0.6* 
0.7*  0,66 
0.76  0,6) 

8:W  8:1? 

0.67  0.74 
0.65  0.71 
0.67  0.6* 
0.5*  0.66 
0.56  0.64 
0.54  0.61 


0.46  O.54 
»J  0.51 
v 1 0.49 
0.47 
0.44 
' 0.W 

0.4<* 
•/  -»4  0.57 


0.22  0,51 
0.1*  0,2* 
0.17  0.26 
0.15  O.H 
0.1*  D.V2 
O.IC  0,20 
0.06  0.16 
0.06  0,16 
0.0)  0.14 
0.01  0.17 


1.6*  7.11 
7.50  7.*? 
1.10  7.1* 
1.94  1**9 

1.65  1.66 

1.76  M9 
I. 6«  1,71 
1.61  1,65 
1.55  1.5# 

i # ift 

1.'*  1.4) 

1.**  i.V) 

l.*0  l.*4 
S.?6  1**0 

1.77  1.76 
1.1*  1.7) 
1.15  1.1* 


1.05  1.10 

1.02  1,07 
0.**  1.04 
0.96  1.01 

0.**  0,*A 

8:21  8:21 
0.65  0.90 
0.6*  0,66 
0.60  0.15 
0.7*  0.6* 
0.76  0.60 
0,75  0.76 

0 .11  0.76 
0.69  0.74 
0,66  0.7? 
0.64  0.70 
0.67  0.67 

0.»0  0.65 

0.56  0.6* 
0.55  0,61 
0.5)  0.59 
0.51  0.57 
0.5*  O.55 
0.47  0.5* 
0.45  0.51 

0,4 6 O.50 
0.41  0.40 
0.50  0.46 
0.)8  0.44 
0. *6  0.62 
0.*4  0.40 
0.  *2  0.*9 
O.JO  0»)7 
0,26  0.)5 
0.26  0.5* 
0.24  0.)1 
0.2J  0.)Q 
0.21  0«?6 
0,19  0.76 
0.17  0,?4 
0.15  0.2* 
0*1)  0.21 
0.12  0.1* 
0.10  0.16 
0.06  0.16 
0.06  0.14 
0.05  0,1) 
0,0)  O.U 
0.01  0.0* 
0.07 


l.l?  7,79 
2. *9  7.41 

7.20  7,72 

7.07  7.0* 
1.96  1.96 
1.67  1.90 

1.60  1.67 
1.1*  1.1* 

1.61  1.70 

1.62  1.65 

1.57  1.60 
1.52  1.55 
1.46  1.91 

1*44  1.47 

1.40  1,4* 
1.)*  1.59 
1.*)  l,*6 

l.JO  l.** 

1.27  l.*0- 

1.74  1.77 

1.21  1.74 
1.16  1.7) 
1.15  1.16 

1.1?  1.16  : 

1.10  ),)* 

1.07  1.11 

1.0*  1.06 
1.09  1,06 

1.00  1.04 
0.96  1,01 
O.X  0.** 
0.**  0.97 
0,*1  0.9) 
0.6*  0.9)  I 
O.97  0.91 
0.65  0.9* 
0.9*  0.67 
0,61  0.65 
0.79  0.6)  ■ 
0.7/  0.61 
0.75  0.79 

0.74  0.76 

0.72  0.74 
0.10  U.74 
0.69  0.7) 
0.67  0,71 
0.65  0.6* 
0.4)  0,67 
0.61  0.66 
0.60  0.64 
0.56  0.6* 
0,56  0.61 
0.55  0,5* 
0.5*  0.56 
0.52  0.56 
0.50  0,55 
0,40  0.5* 
O.47  0.52 
O.45  0.50 

0.44  0.4* 
0.42  0.41 
0.41  0.44 
0.59  0.44 
0.)6  0.4) 
0,)6  0.41 
0.55  0.40 
0.*)  0.59 
0, )?  0.J7 
0.10  o.X 
0.?9  0,  V* 

O.V?  0.)) 
Q.?4  0.)2 
0.25  O.)0 
0.2*  0.79 
0.2?  0.27 
0.20  0.76 

0.1*  0.25 

0.17  0.?* 
0.16  0.2? 
-.15  0.21 
0.1*  0.1* 
0.1?  0.16 
0,10  0.17 
0.09  0,15 
O.OA  0,14 
0.06  0.1) 
0.05  0.11 
0.0'  o.i-> 

C.O?  O.OS 

0 0;  0,07 
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PLOTTING  POSITIONS  IN  PERCENT 

(EXCEEDENCE  FREQUENCY  IN  EVENTS  PER  HUNDRED  YEARS) 
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PLOTTING  POSITIONS  IN  PERCENT 
(EXCEEDENCE  FREQUENCY  IN  EVENTS  PER  HUNDRED  YEARS) 


1.1851  1.1712  1.1576  l.lW  1.1515  1.1186  1.1065  1.0942  1.0624  1.0709  1.0596  1.0485  1.0575  1.0269  1.0165  1.0061  .996? 

For  conputer  application,  the  foranla  (m— is  used. 
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For  use  with  samples  drawn  from  a normal  population 
(See  par.  4-03  d) 
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NOTE:  % values  above  are  usable  approximately  with  Pearson  Type  III 
distributions  having  small  skew  coefficients  * 


EXHIBIT  40 


0 l 2 3 4 5 6 7 8 


101000  OCL  00?  013  017  021  025  02?  O33  O3? 

11  041  C45  049  053  057  061  0£4  06*  072  0?o 

12  079  083  086  090  093  097  100  10li  107  111 

13  'ill;  117  121  121;  127  130  13?  lit 0 143 

lit  146  14$  152  15?  156  161  I6!t  167  170  173 

15  176  179  182  185  186  190  193  196  199  201 

16  20!;  20?  210  212  215  217  220  223  225  228 

17  230  '233  236  238  2ll  213  2lj6  248  250  253 

18  255  258  260  ri(P.  265  26?  270  272  27U  276 

19  279  281  283  286  288  290  292  291;  297  2 99 

20  301  303  305  308  310  312  314  316  318  320 

21  322  321;  326  328  330  332  33!;  336  338  340 

22  3l;2  3 lil;  316  318  350  352.351;  356  356  360 

23  362  361;  366  367  369  371  373  375  37?  378 

24  300  382  381;  306  387  389  3 91  393  394  396 

25  398  LOO  101  L03  L05  LO?  );03  L10  Ll2  L13 

26  L15  L17  118  L20  L22  L23  1*25  L27  L23  L30 

27  L31  1;33  L35  L36  L30  L39  U;i  ;;!;);  LL6 

28  447  449  L50  L52  );53  L55  L56  458  L59  L61 

,29  162  161;  L65  L67  L6S  L70  L71  Ji73  li7lt  4?6 

30  ,’(77  L79  LOO  L31  L83  464  1.86  L87  L89  L90 

31  L9l  L93  L9L  L96  L97  Ls8  500  501  502  SOL 

32  505  507  508  5oo  5ll  512  513  515  516  517 

*"  519  520  521  522  52L  525  526  ?2fi  529  530 

530.  533  531;  535  537  538  539  5L0  SL2  5L3 

544  5L5  5L7  5L8  549  550  551  553  55L  555 

556  558  559  560  561  562  563  565  <66  567 

568  569  571  572  573  57l;  575  576  <T(  579 

580  531  582  583  58L  585  58?  588  589  590 

591  592  593  594  596  597  598  5 99  600  601 

602  6O3  604  605  606  60?  6O9  610  611  612 

613  611;  615  616  617  618  619  620  621  622 

...  623  624  625  626  627  628  629  630  6.31  632 

633  63)1  635  636  637  638  639  640  641  642 

643  644  645  646  64?  648  61,9  650  651  652 

653  654  655  656  657  658  659  660  661  662 1 

663  661;  665  666  667  667  660  669  670  671 

672  673  674  675  676  677  673  679  679  680 

601  682  683  604  60S  606  667  688  68  8 609 

690  691  692  693  694  695  695  696  697  698 

699  700  701  702  702  703  704  705  706  707 

700  700  709  710  711  712  713  7V*  714  705 

716  717  718  719  719  720  721  722  723  723 

721;  725  726  727  728  728  729  730  731  732 

732  733  734  735  736  736  737  738  739  740 


55  740  741  742  743  744 

56  748  749  750  751  751 

57  756  757  757  758  759 
53  765  764  765  766  766 
59  771  772  772  773  774 

778  779  780  760  ?81 
785  766  767  787  788 
792  793  794  794  795 
799  800  801  eOi  802 
806  007  808  800  009 


813  014  814  815  816  I 

820  820  621  822  822  I 
826  82?  82?  828  829  i 
833  833  834  834  835  1 
839  839  840  841  841- 1 

845  846  846  847  840 
651  852  852  653  854 
85?  858  859  859  660 
863  864  865  865  666 
869  870  870  871  872 


875  876  876  877  077 
881  081  682  883  883 
886  887  880  888  889 
892  893  893  894  894 
■098  898  899  899  900 

903  904  904  905  905 
908  909  910  910  911 
914  914  915  915  916 
919  920  920  921  921 
924  925  925  926  926 

929  930  930  931  931 
935  935  936  936  9? 7 
9);0  940  941  941  942 
914  945  9l£  946  946 
91*9  950  950  953  951 

954  955  955  956  956 
S59  960  960  960  961 
964  964  965  965  966 
968  969  969  970  970 
973  974  974  975  975 

978  97e  979  979  900 
902  983  983  901;  984 
987  987  988  988  989 
991  992  992  993  993 
996  996  997  997  997 


HOlEi  Logarithm  of  a number  is  tho  algebraic 
sum  of  above  mantissas  (preceded  by 
a decimal  point)  and  tho  characteristic. 
The  characteristic  is  one  loss  than  tho 
number  of  digits  left  of . the  decimal  point 
if  number  is  one  or  greater.  If  number 
is  less  than  1,  characteristic  is  negative 
and  ecjual  to  one  more  than  the  number  of 
consecutive  ciphers  ismedia teiy  following 
the  decimal  point, 


Examples 


Number 

4380 

2.1 

0.21 

0.01; 


